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Overview 

Already in 19th century Euler and Lagrange established a mathematically satis- 
factory foundation of Newtonian mechanics. Hamilton developed analogous for- 
mulation of optics. Jacobi imported Hamilton's idea in mechanics, and eventually 
arrived at a new formulation (now referred to as the Hamilton -Jacobi formalism). 
The Hamilton-Jacobi formalism was a crucial step towards Liouville's classical 
definition of the notion of integrability. Moreover, Jacobi himself is also famous 
for having discovered (1839) that the geodesic motion on an ellipsoid is an inte- 
grable system which is solvable in terms of hyper-elliptic functions [7]. Jacobi's 
research on dynamical started in 1837 and deeply motivated by Hamilton's for- 
mulation of optics based on the least action. One may say that Hamilton's work 
played the same role as Abel's work in Jacobi's researches of elliptic functions. 

Liouville's definition of integrability [8] is based on the notion of first integrals, 
i.e. conservation laws. In his definition, a (Hamiltonian) system is said to be 
integrable if it has sufficiently many first integrals in involution. The same idea 
has ever been inherited in many variants of the notion of integrability. Liouville's 
definition of integrable Hamiltonian systems naturally covered many classical ex- 
amples. Among them are the Kepler motion solved by Newton, harmonic oscil- 
lators solvable by trigonometric functions, the rigid bodies ("spinning tops") of 
the Euler-Poisont type and the Lagrange type, and Jacobi's example of geodesic 
motion on an ellipsoid. The spinning tops and Jacobi's example were significant 
because they were known to be solvable by elliptic functions. Soon after the work 
of Liouville, C. Neumann discovered a new integrable Hamiltonian system, and 
pointed out that this system can be solved by hyper-elliptic functions. That was 
the beginning of subsequent discoveries of many integrable systems. Neumann's 
Hamiltonian system [9] was a prototype of integrable Hamiltonian systems dis- 
covered in the 19th century. They all are more or less connected with hyper- 
elliptic functions. 

Some generalizations of Langrange's rigid body problem (spinning top) has been 
done by G. Kirchhoff [10] who derived the equation of motion of a rigid body in 
an ideal fluid and by A. Clebsch [11] and V. Steklov [12] for some other cases. 
Weber [13] solved the motion of a rigid body in an ideal fluid in terms of genus- 
two hyper-elliptic functions. Soon afterwards Sophie Kowalevski [14] discovered 
her famous the third example of solvable rigid bodies problem and solved it in 
terms of genus-two hyper-elliptic functions. Triumphal procession of classical 
theory of integrable systems has been crowned by dissertation of P. Stackel [15] 
where a systematic classification of Hamiltoninan systems is presented that can 
be solved by separation of variables. 
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All the researches, mentioned above, belong to the main stream of studies on inte- 
grable systems in the 19th century. They are, however, not the direct origin of the 
breakthrough in the seventies of the 20th century. The breakthrough originates in 
quite different sources in 19th century — study of solitary waves (or "solitons" in 
the modem language) and surface geometry. Studies on solitons were initiated by 
J. Scott Russell when he observed a solitary wave in a canal in 1834. His report 
invoked a controversy. Sir G.B. Airy was suspicious, but Lord Rayleigh advocated 
Scott Russel's observation by calculation (1876). Finally, in 1895, D.J. Korteweg 
and G. de Vries proposed a nonlinear partial differential equation (the KdV equa- 
tion) that well fits in the result of Lord Rayleigh. Note that, unlike the equation 
of motion of point particles and rigid bodies, the KdV equation is a partial differ- 
ential equation. In other words, it has an infinite degrees of freedom. This is a 
common feature of many soliton equations discovered in the breakthrough in the 
seventies. Surface geometry is also a subject started and developed in the 19th 
century. The sine-Gordon equation, which is also an important soliton equation, 
is a partial differential equation that characterizes a special family of surfaces. 
Furthermore, members of this family are connected by the so called Backlund 
transformations. A good classical overview of this subject is provided in G. Dar- 
boux's book [16], in which one can even find the Toda lattice! 

The brilliant success of the search for integrable Hamiltonian systems was rapidly 
fading on the tum to the century. It was the beginning of a long blank that contin- 
ued until the sixties. Routes towards the breakthrough after that blank, however, 
were already prepared in these days. A route was opened by Paul Painleve and 
the contemporaries. Actually, the idea of Painleve was conceptually very close 
to Kovalevskaya's work on integrable rigid bodies. Kovalevskaya's method was 
to search for the cases where the solutions of the equation of motion, analytically 
continued to the complex plane, have no singularity other than poles. Painleve 
did the same for second order nonlinear ordinary differential equations, slightly 
relaxing the conditions on possible singularities of solutions. It should be noted 
that the celebrated six Painleve equations (and presumably their various gener- 
alizations) are NOT integrable in the classical sense. Rather, apart from special 
cases, they are not solvable by any abelian function, nor reducible to a linear 
ordinary differential functions; this property is known as the irreducibility of so- 
lutions of the Painleve equations. Nevertheless Gamier, one of the successors of 
Painleve's researches, pointed out a link with integrable systems [17]. In this pa- 
per Gamier discovered an integrable system solvable by hyper-elliptic functions 
as a byproduct of his researches on a generalization of Painleves equations. This 
system was derived by taking a special limit of his generalized Painleve equations. 
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Remarkably, his calculations implicitly uses a "Lax equation", which is a clue in 
the present approach to integrable systems. 

Another route to the revival of integrable systems was discovered in a quite differ- 
ent direction by Drach, Burchnall and Chaundy. Drach considered [18], [19] the 
second order linear ordinary equation 

^ = [Hx) + h]y (1) 

from a Galois-theoretic point of view, and discovered the case where the coeffi- 
cient (jc) and the solution ^(x) are both related to hyper-elliptic abelian functions. 
J.L. Burchnall and T.W. Chaundy considered [20] the case where the same linear 
ordinary operator L as Drach's commutes with another ordinary differential op- 
erator M, namely satisfies the commutation relation [L,M] = 0, and got results 
similar to Drach's work. In fact, Drach, Burchnall and Chaundy discovered or fi- 
nite gap operators, related to hyper-elliptic solutions of the KdV equations. Their 
work is thus a precursor of algebraic-geometric studies of soliton equations in the 
seventies. 

Theory of integrable systems resurrected as soliton theory in 1965 when M. Zabus- 
ky and M.D. Kruskal reported the celebrated numerical computation of solutions 
of the KdV equation. Their research was motivated by the work of E. Fermi, J. 
Pasta and S. Ulam (FPU) done in the fifties. Zabusky was re-examining this work, 
approximating the nonlinear lattice of FPU by a continuous system. The continu- 
ous system turned out to be described by the KdV equation. He and Kruskal thus 
attempted to solve the KdV equation numerically, and discovered numerical solu- 
tions in which many solitary waves co-existed. The numerical solutions revealed 
remarkable stability of the solitary waves, each of which behaved like a "particle". 
Because of this behavior, they called the solitary waves solitons. This observa- 
tion stimulated theoretical researches, and soon led to the discovery, by Gardner, 
Greene, Kruskal and Miura [21] of exact multi-soliton solutions and the inverse 
scattering method that produces those solutions. The clue of their discovery was 
a relation between the KdV equation and the stationary Schrodinger equation in 
one spatial dimension. Using this relation, they could derive the multi-soliton so- 
lutions, an finite number of conservation laws, etc. systematically. Peter Lax [22] 
soon proposed a more convenient and universal reformulation of the work [21] 
which is now called the Lax formalism. 

The Toda lattice was also discovered in these days. This discovery, however, was 
done independently, [88, 89]. 
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Being unaware of the work of FPU, Toda was looking for an exact model of heat 
conduction, and eventually arrived at his exponential lattice. His method was to 
go back the usual way from equations to solutions: starting from a nonlinear wave 
given by an elliptic function, he attempted to derive a nonlinear equation to be 
satisfied by the nonlinear wave, and eventually arrived to his exponential lattice. 
Within less than a decade, the method of [21] and [22] was extended to many other 
soliton equations, some of the most results in this direction are briefly presented 
below. 

Ablowitz et al., [75], extended the inverse scattering method to a matrix system. 
A wide range of new equations, such as the modified KdV equation, the nonlinear 
Schrodinger equation, and the classical sine-Gordon equations, thus turned out to 
be solvable by the inverse scattering method, [80, 80, 83] Flaschka and Manakov 
independently developed a Lax formalism of the Toda lattice. The inverse scatter- 
ing method was thus further extended to a spatially discrete system. 

Zakharov and Shabat, [92], extended the inverse scattering method to systems 
with two spatial-dimensions. A typical example was the Kadomtsev-Petviashvili 
(KP) equation. While the inverse scattering method was thus refined, several new 
techniques were also invented in the seventies. In particular, the direct method 
(bilinearization) of R. Hirota, the algebro-geometric method of B.A. Durovin, 
V.B. Matveev, S.P. Novikov and I.M. Krichever, and the group-theoretical (or Lie- 
algebraic) method due to M. Adler, B. Kostant, WW. Symes, A.G. Reyman and 
M.A. Semenov-Tian-Shansky emerged in the second half of the seventies, and 
grew up to the mainstream of the progress through the beginning of the eighties. 

Yet another route — Calogero systems The world of integrable Hamiltonian sys- 
tems of the 19th century, too, came back soon after the birth of the soliton theory, 
with a novel family of integrable systems — the Calogero systems. The Calogero 
systems are Hamiltonian systems of interacting particle on a line. Calogero dis- 
covered the simplest case of these systems as a quantum integrable system in [78]. 

This paper deals with the quantum integrability of the system with two-body po- 
tential l/x2. Sutherland soon proposed another version of Calogero's systems in 
[87] where he considered a variant of Calogero's system with two-body potential 

l/sin2x. 

This system is nowadays called the Sutherland system. Calogero conjectured the 
integrability of classical analogues of these systems, and presented a partial an- 
swer. Moser solved this probelm by constructing a Lax formalism of these sys- 
tems. 
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Calogero and Marchioro proved in [79] the classical integrability of the three- 
body case by explicit calculation. Moser, [84], constructed a Lax representation 
for both the Calogero and Sutherland systems, and proved the classical integrabil- 
ity (at least for the Calogero system). In fact, this paper is also known to be one 
of the papers in which Moser presented his method for sovling the non-periodic 
finite Toda (and Kac-van Moerbeke) system using Stiektjes' theory of continued 
fraction. Moser's work suggested that the Lax formalism, originally developed 
in soliton theory, would be also useful for finite-dimensional integrable systems. 
Moser's idea was further extended by Olshanetsky and Perelomov in [86]; this 
work is also one of the earliest attempts of the Lie algebraic methods for construct- 
ing (and solving) integrable systems. A similar work was done by Bogoyavlensky 
for the (finite) Toda lattice, [76]. 

After the aforementioned work, Moser turned to classical integrable systems of the 
19th century, i.e., Jacobi's geodesic flows on an ellipsoid, C. Neumann's system, 
etc., and demonstrated that these systems, too, can be treated in the Lax formalism, 
[85]. 



(...) 
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Chapter 1 



General notions and ideas 



1.1 Introduction 

In our first Chapter we are going to present some general notions and ideas of 
modem theory of integrability trying to outline its computability aspects. The 
reason of this approach is that though this theory was wide and deeply developed 
in the last few decades, its results are almost unusable for non-specialists in the 
area due to its complexity as well as due to some specific jargon unknown to math- 
ematicians working in other areas. On the other hand, a lot of known results are 
completely algorithmic and can be used as a base for developing some symbolic 
programm package dealing with the problems of integrability. Creation of such a 
package will be of a great help not only at the stage of formulating of some new 
hypothesis but also as a tool to get new systematization and classification results, 
for instance, to get complete lists of integrable equations with given properties as 
it was done already for PDEs with known symmetries [31]. 

Some results presented here are quite simple and can be obtained by any student 
acquainted with the basics of calculus (in these cases direct derivation is given) 
while some ideas and results demand deep knowing of a great deal of modem 
mathematics (in this cases only formulations and references are given). Our main 
idea is not to present here the simplest subjects of integrability theory but to give 
its general description in simplest possible form in order to give a reader a feeling 
what has been done already and what could/should be done further in this area. 
Most of the subjects mentioned here will be discussed in details in the next Chap- 
ters. 
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CHAPTER 1. GENERAL NOTIONS AND IDEAS 



We will use the word "integrability" as a generalization of the notion "exactly 
solvable" for differential equations. Possible definitions of differential operator 
will be discussed as well as some definitions of integrability itself. Numerous ex- 
amples presented here are to show in particular that it is reasonable not only to use 
different notions of integrability for different differential equations but sometimes 
it proves to be very useful to regard one equation using various definitions of in- 
tegrability, depending on what properties of the equation are under the study. Two 
classical approaches to classification of integrable equations - conservation laws 
and Lie symmetries - are also briefly presented and a few examples are given in 
order to demonstrate deep difference between these two notions, specially in case 
of PDEs. Two interesting semi-integrable systems are introduced showing one 
more aspect of integrability theory - some equations though not integrable in any 
strict sense, can be treated as "almost" integrable due to their intrinsic properties. 



1.2 Notion of differential operator 

There are many ways to define linear differential first-order operator D, starting 
with Leibnitz formula for product differentiating 

D{ab)^a'b + ab'. (1.1) 

This definition leads to 

Da = a +aD, D'^a = aD^ + 2a D + a\...., (1.2) 

• « = L ( ^ ) D^i^W^ (1-3) 

where 

( n\ n{n-l)...{n-k+\) .^^ (n\ , 
[k)= 1.2... k [o)=' 

are binomial coefficients. Thus, one linear differential first-order operator and its 
powers are defined. Trying to define a composition of two linear operators, we get 
already quite cumbersome formula 



DiD2a^Di{D2{a)+aD2)^Di{a)D2 + DiD2{a) +D2{a)Di + aDiD2, (1.5) 



which can be regarded as a definition of a factorizable linear differential second- 
order operator. Notice that though each of Di and D2 satisfies Leibnitz rule, their 
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composition D1D2 does not! An important notion of commutator of two oper- 
ators [D1D2 — D2D1] plays a role of special multiplication (see Ex.1). In case 
of non-factorizable second-order operator an attempt to generalize Leibnitz rule 
leads to a very complicated Bourbaki-like constructions which we are not going 
to present here. All these problems appear due to coordinateness of this approach. 

On the other hand, a linear differential operator being written in coordinate form 
as 

^ = IM- (1-6) 

j 

leads to 

L= fad", a« = a«i---a««, \a\ = ai + ... + an (1.7) 

\a\<m 

which is definition of LPDO of order m with n independent variables. Composi- 
tion of two operators L and M is defined as 

LoM = Y^fad^Y^gpdl^ = Y^hjdr (1.8) 

and coefficients hy are to be found from formula (1.3). 

Now, notion of linear differential equation (LDE) can be introduced in terms of 
the kernel of differential operator, i.e. 

Ker(L):={<p| /«a«<p = 0}. (1.9) 

\a\<m 

Let us regard as illustrating example second-order LODO with one independent 
variable x: 

L = fo+fid+f2d\ d:^^. (1.10) 

ax 

Notice that if two functions 91 and ^ belong to its kernel Ker(L), then 

ci9i+C2^ e Ker(L), (1.11) 
i.e. Ker(L) is a linear vector space over constants' field (normally, it is M or C). 

Main theorem about ODEs states the existence and uniqueness of the Cauchy 
problem for any ODE, i.e. one-to-one correspondence between elements of the 
kernel and initial data. In our case Cauchy data 

(p\x=xo = <P°, d(p\x=xo = (P^ (1.12) 
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form a two-dimensional vector space and, correspondingly, dimension of kernel 
is equal 2, J/m(Ker(L)) = 2. Any two functions 9i , ^ G Ker(L) form its basis if 
Wronskian < ^i, ^ > is non- vanishing: 



< (Pi,9z>-= 



7^0 



(1.13) 



while an arbitrary function ^^ e Ker(L) has to satisfy following condition: 



< (pu(p2,\\f >:= 



(pi (p2 ¥^ 

II II II 



= 0. 



(1.14) 



Now we can construct immediately differential operator with a given kemel as 

L(va) =< <pi,<j?2,V^> . (1.15) 

For instance, if we are looking for an LODE with solutions sinjc and \fx, then 
corresponding LODE has form 



In '13 1^ 

iJ/ (1 tanx)+tan^w iJ/ -im- = 0. 

2 Lx Ax^ 



(1.16) 



Coming back to LODO of order m, we re-write formula for the kemel in the form 



L(VA) = 



< <pi,...,<j[^,y > 
< (pi,...,(pn> 



(1.17) 



which provides that high-order coefficient = 1. It is done just for our conve- 
nience and we will use this form further. 



Now, some constructive definition of linear differential operator was given and 
importance of its kemel was demonstrated. Corresponding differential equation 
was defined in terms of this kemel and constmction of operator with a given kemel 
was described. All this is not possible for nonlinear operator because in this case 
manifold of solutions has much more complicated structure then just linear vector 
space - simply speaking, the reason of it is that in this case linear combination of 
solutions is not a solution anymore. Due to this reason some other notions are to 
be used to study properties of nonlinear operators - symmetries, conservation laws 
and, of course, as the very first step - change of variables transforming a nonlinear 
operator into a linear one. We will discuss all this in the next sections. 
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1.3 Notion of integrability 

• 1.3. Solution in elementary functions: 

y"+y = 0. (1.18) 

General solution of this equation belongs to the class of trigonometrical 
functions, y — asm{x + b), with arbitrary const a^b. In order to find this 
solution one has to notice that this equation is LODE with constant coeffi- 
cients which possess fundamental system of solutions, all of the form e^^ 
where A is a root of characteristic polynomial. 

• 1.2. Solution modulo class of functions: 

1.2.1. Integrability in quadratures 

y" = f{y)- (1.19) 

In order to integrate this equation let us notice that y"y' — f{y)y' which 
leads to 

= J fiy)dy + const = Fiy) (1.20) 

and finally 

dx= -., (1.21) 

which describes differential equation with separable variables. In 

case when F(y) is a polynomial of third or fourth degree, this is defini- 
tion of elliptic integral and therefore the initial nonlinear ODE is inte- 
grable in elliptic functions. Particular case when polynomial has 
multiple roots might leads to a particular solution in elementary func- 
tions. Let us regard, for instance, equation /' = 2y^ and put const = 0, 
then = y'^ and y = l/x, i.e. we have ONE solution in the class of 
rational functions. General solution is written out in terms of elliptic 
functions. Conclusion: equation is integrable in the class of elliptic 
functions and not integrable in the class of rational functions. 



The case when F{y) is a not-specified smooth function, is called in- 
tegrability in quadratures. This is collective name was introduced 
in classic literature of 18th century and used in a sense of "closed- 
form" or "explicit form" of integrability while analytic theory of spe- 
cial functions was developed much later. 
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Notice that as a first step in finding of solution, the order n of initial 
ODE was diminished to n — 1, in our case n = 2. Of course, this is 
not possible for any arbitrary differential equation. This new ODE 
is called first integral or conservation law due its physical mean- 
ing in applications, for instance, our example can be reformulated as 
Newton second law of mechanics and its first integral corresponds to 
energy conservation law. 

1.2.2. Painleve transcendent PI 



This equation defines first Painleve transcendent. About this equa- 
tion it was proven that it has no solutions in classes of elementary or 
special functions. On the other hand, it is also proven that Painleve 
transcendent is a meromorphic function with known special qualita- 
tive properties ([25]). 

This example demonstrates us the intrinsical difficulties when defining 
the notion of integrability. Scientific community has no general opin- 
ion about integrability of Painleve transcendent. Those standing on 
the classical positions think about it as about non-integrable equation. 
Those who are working on different applicative problems of theoreti- 
cal physics involving the use of Painleve transcendent look at it as at 
some new special function (see also Appendix 1). 

• 1.3. Solution modulo inexplicit function: 



This equation describes so called shock waves and its solutions are ex- 
pressed in terms of inexplicit function. Indeed, let us rewrite this equation 
in the new independent variables i = t, x = u and dependent one u=x, i.e. 
now X = G{t,u) is a function on t,u. Then 



y =y +x. 



(1.22) 



Ut — luUx- 



(1.23) 



dt = dt, dx= Ofdt + Oudu, Ut\ dx=Q = 



. 1 

Ux\dt=0 = TT 



'u 



(1.24) 



and 




6 = 2ut-(p{u) 



(1.25) 
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and finally x+ Itu = q){u) where q){u) is arbitrary function on u. Now, we 
have finite answer but no explicit form of dependence u = u{xj). Has the 
general solution been found? The answer is that given some initial condi- 
tions, i.e. t — 0,we may define solution as w = (p~^{x) where is cp^^ denotes 
inverse function for (p. 

1.4. Solution modulo change of variables (C-Integrability): 

\l/xy + a\l/x + ^\l/y + \l/xWy = 0- (1-26) 

This equation is called Thomas equation and it could be made linear with 
a change of variables. Indeed, let ij/ = logd for some positively defined 
function 6: 

Vx^{lo8d)x=j, xify = {logd)y = ^, (1.27) 

¥xy-[^>y- ^2 ' y^xWy--^ (i-^o) 

and substituting this into Thomas equation we get finally linear PDE 

exy+aex+pey = o. (1.29) 

Suppose for simplicity that j8 = and make once more change of variables: 
d = (j)e^'y, then 

0x = 0xe^'^ exy = (l>xy^'^ + ki<l>xe^^y and 0xy+ (^i + a)0x = 0, 

(1.30) 

and finally 

dx{^y + {ki + a)<^)^Q, (1.31) 

which yields to 

0y-^20 = /(y) with k2 = -{ki + a) (1.32) 

and arbitrary function f{y). Now general solution can be obtained by the 
method of variation of a constant. As a first step let us solve homogeneous 
part of this equation, i.e. 0^ — ^^2^ = and ^{x,y) = g{x)e^^y with arbitrary 
g{x). As a second step, suppose that g{x) is function on jc,_y, i.e. g{x,y), 
then initial equation takes form 

{g{x,y)e^^y)y-k2g{x,y)e'^^y = f{y), (1.33) 
g{x,y)ye''^y + k2g{x,y)e''^y -k2g{x,y)e'^y = f{y), (1.34) 

g{x,y)ye''y = f{y), g{x,y) = J f{y)e-''ydy + h{x), (1.35) 
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and finally the general solution of Thomas equation with a = can be writ- 
ten out as 

Hx,y) ^g{x,y)e'^y = e'^y{ J f{y)e-''^ydy + h{x)) ^ f{y)+e'^yh{x) 

(1.36) 

with two arbitrary functions f{y) and h{x). 
• 1.5. Solution modulo Fourier transform (F-Integrability): 

ut = luux + euxx (1-37) 

where e is a constant. This equation is called Burgers equation and it 
differs from Thomas equation studied above where change of variables was 
local in a sense that solution in each point does not depend on the solution 
in some other points of definition domain, i.e. local in (jc,_y) -space. This 
equation can be transformed into 

Ut^2uux + Uxx (1.38) 

by the change of variables x = ex, u = eu and t = s^t. We will use this form 
of Burgers equation skipping tildes in order to simplify the calculations be- 
low. To find solution of Burgers equation one has to use Fourier transform 
which obviously is nonlocal, i.e. here solution is local only in ^-space. In 
order to demonstrate it let us integrate it using notation 



udx — v, (1-39) 
then after integration 

vt = vl + vxx = e-'ie')xx (1.40) 

and change of variables w — e^. Burgers equation is reduced to the heat 
equation 

Wt = Wxx (1.41) 

which is linear. Therefore, solutions of Burgers equation could be obtained 
from solutions of heat equation by the change of variables 

Wr 

u^Vx = — (1.42) 
w 



as 



/oo 
Qii^{ikx)w{k)dk (1.43) 
-oo 

, , J exp{ikx — k^t)w{k)ikdk (144) 
' Jexp{ikx — k'^t)w{k)dk ' 
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where w(x,0) is initial data and w{k) is a function called its Fourier trans- 
form and it can be computed as 

1 f°° 

w{k) = — / exp{—ikx)w{x,0)dx. (1-45) 
271 J -co 

In fact, it is well-known that any linear PDE with constant coefficients on an 
infinite line can be solved using as standard basis {e^'^\k e m} because they 
are eigenfunctions of these operators. Thus, Thomas equation with a = 
where the general solution was found explicitly, is an exception while heat 
equation demonstrates the general situation. 

• 1.6. Solution modulo 1ST (S-Integrability): 

ut = 6uux + Uxxx ■ (1-46) 

This equation is called Korteveg-de Vries (KdV) equation and it is non- 
linear PDE with nonconstant coefficients. In this case, choosing set of func- 
tions {e'^} as a basis is not helpful anymore: Fourier transform does not 
simplify the initial equation and only generates an infinite system of ODEs 
on Fourier coefficients. 

On the other hand, some new basis can be found which allows to reduce 
KdV with rapidly decreasing initial data, m ^ 0, jc ^ ±oo, to the linear equa- 
tion and to solve it. This new basis can be constructed using solutions of 
linear Schrodinger equation 

\l/xx + k^¥ = W (1.47) 

where function u is called potential due to its origin in quantum mechan- 
ics. Solutions of linear Schrodinger equation are called Jost functions, 
\l/'^{t,x,k), with asymptotic boundary conditions: 

y/^{t,x,k;u{x,t))e-^'^^'+''''^ ^1, jc^ ±oc. (1.48) 

Jost function (^{x,k) — \{f~^{t,x,k;u{x,t))e~'^'^'^'^*'> is defined by the integral 
equation 

(l){x,k) = 1+ ri:i^^^P|fc:^„(x')0(x',;fc)jx' (1.49) 
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with t playing role of a parameter. Second Jost function is defined analo- 
gously with integration over [— . Notice that asymptotically for x — )■ ±00 
linear Schrodinger equation 

Yxx + k^V=W is reduced to \i/xx + k^\i/ = (1.50) 

as in case of Fourier basis {e'^^\k E r}. It means that asymptotically their 
solutions do coincide and, for instance, any solution of linear Schrodinger 
equation 

V^~cie'(^+^')+C2e-'(^+^'^ for x^oo. (1.51) 

It turns out that solutions of KdV can be regarded as potentials of linear 
Schrodinger equation, i.e. following system of equations 

Ut = 6uUx + Uxxx, 

is consistent and any solution of KdV can be written out as an expansion of 
Jost functions which in a sense are playing role of exponents e'*^ in Fourier 
basis, [35]. 



On the other hand, there exists a major difference between these two basis: 
Fourier basis is written out in explicit form via one function while Jost basis 
is written out in inexplicit form via two functions with different asymptotic 
properties on the different ends of a line. The crucial fact here is that two 
Jost functions are connected by simple algebraic equation: 

Y~{x,k,t) =a{k)^f'^{x,-k,t) + b{k)e'^^'^f^{x,k,t) (1.53) 

while it allows us to construct rational approximation of Jost functions for 
given a{k) and h{k) and, correspondingly, general solution of KdV. The 
problem of reconstruction of function u according to a{k),b{k) is called in- 
verse scattering problem and this method, correspondingly, inverse scat- 
tering transform (1ST). 

• 1.7. Solution modulo Dressing method (D-Integrability): 

iut^Uxx±\u\^u. (1-54) 

This equation is called nonlinear Schrodinger equation (NLS) and it is 
very important in many physical applications, for instance, in nonlinear op- 
tics. Dressing method is generalization of 1ST and in this case role of auxil- 
iary linear equation (it was linear Schrodinger equation, second order ODE, 
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in the previous case) plays a system of two linear first order ODEs [62]: 

■i>^''";r^'i (1.55) 

where v = ±m. It turns out that system of equations 

is consistent and is equivalent to Riemann-Hilbert problem. Solutions of 
this last system are called matrix Jost functions and any solution of NLS 
can be written out as an expansion of matrix Jost functions which in a sense 
are playing role of exponents e''^ in Fourier basis, [34]. 



1.4 Approach to classification 

Our list of definitions is neither full nor exhaustive, moreover one equation can 
be regarded as integrable due to a few different definitions of integrability. For 
instance, equation for shock waves from § 1.3, Ut = 2uUx, is a particular form of 
Burgers equation from § 1.5, Ut = luu^ + eu^^ with e = 0, and it can be linearized 
as above, i.e. it is not only integrable in terms of inexplicit function but also 
C-integrable and F-integrable, with general solution 

, , J exp{ikx — k^t)i2{ek)ikdk (157) 
J exp{ikx — k^t)u{£k)dk 

What form of integrability is chosen for some specific equation depends on what 
properties of it we are interested in. For instance, the answer in the form of in- 
explicit function shows immediately dependence of solution form on initial con- 
ditions - graphically presentation of inverse function u — (p^^{x) can be obtained 
as mirrored image of x = q>{u). To get the same information from the formula 
above is a very nontrivial task. On the other hand, the general formula is the only 
known tool to study solutions with singularities. This shows that definitions of 
integrability do not suit to serve as a basis for classification of integrable systems 
and some more intrinsic ways should be used to classify and solve them. Below 
we present briefly two possible classification bases - conservation laws and Lie 
symmetries. 
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1.4.1 Conservation laws 

Some strict and reasonable definition of a conservation law (also called first in- 
tegral for ODEs) is not easy to give, even in case of ODEs. As most general 
definitions one might regard 

j/iy)-0, (1.58) 

^y = /(j), y = {yi,-,yn), /(j) = (/i, -,/«), 

d f 

— J G{u,Ux,Uy,Uxx,Uxy,Uyy...)dxdy... = 0, (1.59) 



for ODE 



and 



forPDE 

= S{^1^XlUy,Uxx-,Uxy,Uyy, ...). 

Obviously, without putting some restrictions on function F or G these definitions 

are too general and do not even point out some specific class of differential equa- 
tions. For instance, let us take any second order ODE, due to well-known theorem 
on ODEs solutions we can write its general solution in a form 

F{t,y,a,b) = Q (1.60) 

where a^b are two independent parameters (defined by initial conditions). Theo- 
rem on inexplicit function gives immediately 

a = Fi{t,y,b) V& and b = F2{t,y,a) Va, (1.61) 

i.e. any second order ODE has 2 independent conservation laws and obviously, 
by the same way n independent conservation laws can be constructed for ODE 
of order n. For instance, in the simplest case of second order ODE with constant 
coefficients, general solution and its first derivative have form 

'y = cie^i^ + C2e^2x 
I -1 -1 (1-62) 

and multiplying y by Ai and A2 we get equations on ci and C2 correspondingly: 

A2);-y = ci(A2-Ai)e^i^ ^ (^2X + C2 = log{y'-?iiy) 

hy-y = c2(Ai ->l2)e^2^ [Aijc + ci = log{y' -?^2y) 

and two conservation laws are written out explicitly: 

C2 = log (y -Xiy)- X2X, ci = log {y - X2y) - Aix. ( 1 .64) 
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To find these conservation laws without knowing of solution is more complicated 
task then to solve equation itself and therefore they give no additional informa- 
tion about equation. This is the reason why often only polynomial or rational 
conservation laws are regarded - they are easier to find and mostly they describe 
qualitative properties of the equation which are very important for applications 
(conservation of energy, momentum, etc.) 

On the other hand, conservation laws, when known, are used for construction of 
ODEs solutions. Indeed, let us rewrite Eq.(1.58) as 



-F{y) = ifidi +f2d2 + ... +fndn)F = ^{F) = 0, 



(1.65) 



i.e. as an equation in partial derivatives ^{F) = 0. Such an equation has {n — 
1) independent particular solutions (^i, 9„_i) if its Jacobian matrix has 
maximal rank 



a((Pi,(P2,...,(p„_i) _ 



with notation 



d{y\,y2,-,yn-\) 

d{(pi,(p2,...,CPn-l) 

d{yi,y2,-,yn-i) 



(1.66) 



/ dfi d(pi ... <?<Pi \ 
dy\ '3y2 ' ' ' '3y^ 



. d(pn . . . . 

\ dyi dy2 dyn / 



(1.67) 



Now we can write out the general solution as <j!)„_i) with arbitrary 

function F and initial ODE can be reduced to 

dz 



= dt 



(1.68) 



and solved explicitly in quadratures (see § 1.2). 



1.4.2 Symmetry properties 

In order to give definition of canonical form for n-order ODE 

3;W=F(x,3;,y,. ..,/"-!)) 

let us first introduce vector 



y^ 



X 

y 
y 



(1.69) 



(1.70) 



\y 



(n-l) 
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which is called vector of dynamical variables with all its coordinates regarded 
as independent, and its first derivative 



dy 
dt 



' 1 " 




T 


y 




y 


n 




II 


y 




y 


j(«)_ 




F 



(1.71) 



with respect to some new independent variable t such that dt = dx, then the equa- 
tion 

1 



dy 
dt 



(1.72) 



is called canonical form of an ODE. This canonical form is also called dynam- 
ical system, important fact is that dimension of dynamical system is (n+l) for 
n-order ODE. 



Definition 1.1. Dynamical system 



dy 
dT 



= 8iy) 



is called a symmetry of another dynamical system 

dy 



dt 



if 



d ,dy. d ,dy . d , „ d , 

d-.^i^ = dt^i^ ^ d-.^f^^ = dt^'^^^ 



(1.73) 



(1.74) 



holds. Symmetry g of dynamical system / is called trivial if g = const •/. 



Obviously, Eq.(1.74) gives necessary condition of compatibility of this two dy- 
namical systems. It can be proven that this condition is also sufficient. There- 
fore, construction of each Lie symmetry with group parameter T is equivalent to 
a construction of some ODE which have T as independent variable and is consis- 
tent with a given ODE. System of ODEs obtained this way, when being written 
in canonical form is called dynamical system connected to the element of Lie 
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symmetry algebra. 

We will discuss it in all details later (also see [59]), now just pointing out the fact 
that Lie algebra can be generated not only by Lie transformation group (normally 
used for finding of solutions) but also by the set of dynamical systems (1.72) (nor- 
mally used for classification purposes). 

Two following elementary theorems show interesting interconnections between 
conservation laws and symmetries in case of ODEs. 

Theorem 1.1 Let dynamical system 




is a symmetry of another dynamical system 
and 

jP{y) = {fidx +f2d2 + ... +fndn)F = ^{F) = (1.75) 
is a conservation law. Then Fg is symmetry as well. 

► Indeed, let us introduce 

^ ^{gidi+gidi + .-.+gndn), (1.76) 
then consistency condition Eq.(1.74) is written out as 

and after substituting FM instead of M we get on the left hand of (1.77) 

^{FJi)^^{F)Ji^F^Ji^F^Ji^FJi££, (1.78) 
and on the right hand of 

{FJ^)££ = FM££. (1.79) 
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Corollary 1.2 Let n-th order ODE of the form (1.69) has a symmetry 



dt 



go 
gi 

g2 
gn 



(1.80) 



then go is a conservation law and consequently without loss of generality we may 
putgo= l>ifg0 7^0. 



Theorem 1.3 An ODE of arbitrary order n having of {n— 1) independent con- 
servation laws (i.e. complete set), has no nontrivial symmetries consistent with 
conservation laws. 

► Indeed, using all conservation laws we can reduce original n-order ODE into 
first order ODE 

da , s 

^ = /(a) (1.81) 
and look for symmetries in the form 

^ = g(a). (1.82) 



Then 



and finally 



f ^ ^ ln(l) = const, (1.84) 



/(«) g{a) 
i.e. functions f{a) and g{a) are proportional 



Simply speaking, Eq.(1.73) defines a one-parameter transformation group with 
parameter t which conserves the form of original equation. The very important 
achievement of Lie was his first theorem giving constructive procedure for ob- 
taining such a group. It allowed him to classify integrable differential equations 
and to solve them. The simplest example of such a classification for second order 
ODEs with two symmetries is following: each of them can be transformed into 
one of the four types 



(I) / = h{y) (II) / = h{x), (III) / = (IV) / = h{x)y. 
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where h denotes an arbitrary smooth function while explicit form of correspond- 
ing transformations was also written out by Lie. 

Symmetry approach can also be also used for PDEs. For instance, first-order 
PDE for shock waves Ut = luu^ (§ 1.3) has following Lie symmetry algebra 
Ux = (p{u)ux where any smooth function (p = (p{u) defines one-parameter Lie sym- 
metry group with corresponding choice of parameter T. Indeed, direct check gives 
immediately 

(mj)t= {2q)Ux)ux + 2u{q>Uxx + ^ ul) (1.85) 

and 

(mt)j = 9 Ux2uux + q>{2ux + 2uuxx) (1.86) 

while relation {ut)^ = (mx)? yields to the final answer (here notation (p' = (pu was 
used). In order to construct dynamic system for this equation, let us introduce 
dynamical independent variables as 

M, U\=Ux, U2 = Uxx, ■■■ (1.87) 

and dynamical system as 

d 
dt 

This system can be transformed into finite-dimensional system using characteris- 
tics method, [66]. 

For PDE of order n> I analogous dynamic system turns out to be always infi- 
nite and only particular solutions are to be constructed but no general solutions. 

Infinite-dimensional dynamical systems of this sort are not an easy treat and also 
choice of dynamical variables presents sometimes a special problem to be solved, 
therefore even in such an exhaustive textbook as Olver's [59] these systems are 
not even discussed. On the other hand, practically all known results on classifica- 
tion of integrable nonlinear PDEs of two variables have been obtained using this 
approach (in this context the notion of F-lntegrability is used as it was done to 
integrate Burgers equation, § 1.5.) For instance, in [38] for a PDE of the form 

uxy = f{x,y) (1.89) 

with arbitrary smooth function / on the right hand it was proven that this PDE is 
integrable and has symmetries iff right part has one of the following forms: e". 



u 




2mm 1 


Ml 




2mm2 + 2m J 


M2 




2mM3 +6miM2 



(1.88) 



18 



CHAPTER 1. GENERAL NOTIONS AND IDEAS 



sinw or cie" + C2e ^" with arbitrary constants ci,C2. Another interesting result 
was presented in [66] where all PDEs of the form 



(1.90) 



have been classified. Namely, a PDE of this form is integrable and has symmetries 
iff if it can be linearized by some special class of transformation. General form of 
transformation is written out explicitly. 

1.4.3 Examples 

Few examples presented here demonstrate different constellations of symmetries 
(SYM), conservation laws (CL) and solutions (SOL) for a given equation(s). 

1.4.3.1: SYM +, CL +, SOL +. 



Let us regard a very simple equation 

II 

y =1, 

then its dynamical system can be written out as 



(1.91) 



dt 



"1" 




"1" 


y 




y 


F 




1 



with dt — dx 



(1.92) 



and its general solution is )^ = ^x^ + c\x + C2 with two constants of integration. In 
order to construct conservation laws, we need to resolve formula for solution with 
respect to the constants c\,C2'. 

1 



c\=y -X, C2 =y+-x^' 

Now, we look for solutions F(y) of the equation 
d 



-xy 



dt 



F{y) = {dx + y dy + di)F = ^{F) = 



(1.93) 



(1.94) 



with ^ = dx + y dy + d^i . Direct check shows that F = y —x and F = y+ jx^ — 
xy' are functionally independent solutions of this equation. Moreover, general 
solution is an arbitrary function of two variables 



F = F{y -x,y+ ^^-xy'), 



(1.95) 
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for example, 

F = (y -xf - l{y+ -xy') =y - 2y. (1.96) 

On the other hand, if there are no restriction on the function F, the conservation 
laws may take some quite complicated form, for instance, 

F=Arcsin{y' -x)/{y+^x^ -xy'f'"^^. (1.97) 

Now, that dynamical system, conservation laws and general solution of the origi- 
nal equation have been constructed, let us look for its symmetry: 

giy)- % = siy)^ g{y)^{gi.g2.g2)- (1-98) 

Demand of compatibility 

i,(f(y) - Jt^siy)) is equivalent to ^(gi) = ^(^3) = 0, ^{g2)=g3, 

(1.99) 

and it can be proven that any linear combination of two vectors (1,0,0) and 
{0,Xj 1) with (some) scalar coefficients provides solution of compatibility prob- 
lem. Thus, Lie symmetry group corresponding to the vector (1,0,0) is shift in x 
while the second vector (0,x, 1) corresponds to summing y with particular solu- 
tions of homogeneous equation. 



1.4.3.2: SYM - CL +, SOL +. 



{n2 — n'i)a2a2, (1.100) 
{nj,-n\)a\a2, (1.101) 

(«i-«2)«i«2 (1.102) 

with variables a; and constants This system is well-known in physical ap- 
plications - it describes three-wave interactions of atmospheric planetary waves, 
dynamics of elastic pendulum or swinging string, etc. Conservation laws for this 
system are: 

• Energy conservation law is obtained by multiplying the /-th equation by ai, 
i= 1,2,3 and adding all three of them: 



niai+n2a2 + n2a2 = const. (1.103) 



Let us regard as a system of ODEs 

da I 

"1^ = 
da2 

"^^ = 
da^ 
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• Enstrophy conservation law is obtained by multiplying the i-th equation by 
nitti, i= 1,2,3 and adding all three of them: 

n\a\ + n\a\ + n\a^ — consl. (1.104) 

Using these two conservation laws one can easily obtain expressions for and 
a3 in terms of ai. Substitution of these expressions into the first equation of 
Sys.(l.lOO) gives us differential equation on ai: 

{yf = f{.y) with y^a, (1.105) 

whose explicit solution is one of Jacobian elliptic functions while a2 and a^, are 
two other Jacobian elliptic functions (e.g. [82]). In fact, Sys. (1.100) is often re- 
garded as one of possible definitions of Jacobian elliptic functions. 

Notice that from Theorem 1.3 one can conclude immediately that all symmetries 
consistent with conservation laws, are trivial. 



1.4.3.3: SYM +, CL ±, SOL +. 



Let us regard heat equation 

dtU = dxxU, (1.106) 

which generates solutions of Burgers equation (see § 1.5). Direct check shows that 
Eq.( 1 . 1 06) is invariant due to transformations x — xT and t = with any constant 
T, i.e. dilation transformations constitute Lie symmetry group for Eq.( 1.106). 
Moreover, Eq.(1.106) is integrable and its solution 

u= [ Qxp(ikx — k^t)u(k)dk with u(k) = -;^ f Qxp(—ikx)u(x,0)dx (1.107) 

J LTt J-oo 

is obtained by Fourier transformation. 

This example demonstrates also some very peculiar property - heat equation (as 
well as Burgers equation) has only one conservation law: 

/ udx = 0. (1.108) 

dt j-oc 

Nonexistence of any other conservation laws is proven, for instance, in [66]. 
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1.5 Semi-integrability 
1.5.1 Elements of integrability 

Dispersive evolution PDEs on compacts is our subject in this subsection. In con- 
trast to standard mathematical classification of LPDO into hyperbolic, parabolic 
and elliptic operators there exists some other classification - into dispersive and 
non-dispersive operators - which is successfully used in theoretical physics and is 
not complementary to mathematical one (for details see [82]). Let regard LPDE 
with constant coefficients in a form 

P(|4)=0 (U09) 

where t is time variable and x is space variable, and suppose that a linear wave 

Y{x) =Aexpi{kx- at) (1.110) 

with constant amplitude A, wave number k and frequency CO is its solution. After 
substituting a linear wave into initial LPDE we get P{—io), ik) = 0, which means 
that k and ft) are connected in some way: there exist some function / such, that 
/(«,^) = 0. 

This connection is called dispersion relation and solution of the dispersion rela- 
tion is called dispersion function, (O = (o{k). If condition 

^^0 (1.111) 

holds, then initial LPDE is called evolution dispersive equation and its solutions 
obviously are completely defined by the form of dispersive function. All these 
definitions can be easily reformulated for a case of more space variables, namely 
x\,X2, ...,Xn. In this case linear wave takes form 

= Aexp/(H- (0?) (1.112) 

with wave vector ^ = (A:i , . . . .^„) and space-like variable x= (xi , . . . , x„) . Then 

'(T,'-k-7;r) = ^, (1.113) 

dispersion function can be computed from P{—i(0,ik\,...,ikn) =0 and the condi- 
tion of non-zero second derivative of the dispersion function takes a matrix form: 
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Notice now that solutions of linear evolution dispersive PDE are known by defi- 
nition and the reasonable question here is: what can be found about solutions of 
nonlinear PDE 



with dispersive linear part ^{^f) and some nonlinearity ^(i/a)? 

Nonlinear PDEs of this form play major role in the theory of wave turbulence and 
in general there is no final answer to this question. Case of weak turbulence, i.e. 
when nonlinearity ^( V^) is regarded small in a sense that wave amplitudes A are 
small enough (smallness of an amplitude can be strictly defined), is investigated 
in much more details. Two qualitatively different cases have to be regarded: 

1. coordinates of wave vector are real numbers, {k = {ki, ....kn)\ki e M} 
(corresponds to infinite space domain); 

2. coordinates of wave vector are integer numbers, (k = {ki, ....kn)\ki e Z} 
(corresponds to compact space domain). 

In the first case method of wave kinetic equation has been developed in 60-th 
(see, for instance, [27]) and applied for many different types evolution PDEs. Ki- 
netic equation is approximately equivalent to initial nonlinear PDE but has more 
simple form allowing direct numerical computations of each wave amplitudes in 
a given domain of wave spectrum. Wave kinetic equation is an averaged equation 
imposed on a certain set of correlation functions and it is in fact one limiting case 
of the quantum Bose-Einstein equation while the Boltzman kinetic equation is its 
other limit. Some statistical assumptions have been used in order to obtain kinetic 
equations and limit of its applicability then is a very complicated problem which 
should be solved separately for each specific equation, [36]. 

In the second case, it is proven, [82], that solving of initial nonlinear PDE can be 
reduced to solving a few small connected systems of ODEs of the form 



i^(V/)=^(V/) 



(1.115) 




(1.116) 



in case of quadratic nonlinearity. 



< 



'Ai = TAIA3A4, 
A2 = TAIA3A4, 
A3 = -TA1A2AI 
A4 = -TA1A2AI 
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in case of cubic nonlinearity and so on. Some of the resulting systems are proven 
to be integrable, e.g. [77, 90, 91]. 

Notice, that in contrast to a linear wave with a constant amplitude A ^ A{t,x), 
waves in nonlinear PDE have amplitudes A, depending on time. It means that so- 
lutions of initial nonlinear PDE have characteristic wave form as in linear case 
but wave amplitudes can be expressed in Jacobian elliptic functions on time, 
cn{T),dn{T) and sn(T). Notice that Sys.(1.116) has been studied in § 1.3.2 and 
its conservation laws were found. Exact solutions of Sys.(1.116) can be written 
out explicitly as functions of initial conditions (see [77] for details). 



1.5.2 Levels of integrability 

Let us formulate classical three-body problem whose integrability attracted at- 
tention of many investigators beginning with Lagrange. Computing the mutual 
gravitational interaction of three masses is surprisingly difficult to solve and only 
two integrable cases were found. For simplicity we regard three-body problem 
with all masses equal, then equations of motion take form 

^^=212/12 + 213/13 (1-117) 
-^ = 221/12+223/23 (1.118) 

-^ = 231/13+232/23 (1-119) 

where zj is a complex number, zj — xj + iyj, describing coordinates of j-th mass 
on a plane, fj^ is a given function depending on the distance between j-th and 
^-th masses (physically it is attraction force) while following notations are used: 

Zjk = 2; - Zk and fjk = f{\zjk\^). 

This system admits following conservation laws: 

• Velocity of center of masses is constant 
Summing up all three equations, we get 

^(2l+22 + 23) = 0. (1.120) 

This equality allows us to choose the origin of coordinate system in such a 
way that 

21 +22 + 23 = 



24 



CHAPTER 1. GENERAL NOTIONS AND IDEAS 



which simplifies all further calculations significantly. That is the reason 
why till the end of this section this coordinate system is used. Physically it 
means that coordinate system is connected with masses center 

• Conservation of energy 

Multiplying j-th equation by Zj, summing up all three equations and adding 
complex conjugate, we obtain on the left 

3 J 3 

i.e. left hand describes derivative of kinetic energy. 
On the right we have derivative of potential energy : 

j^'^ ^j^iFi\zn\^)+F{\zi3\^)+F{\z23\^)) with notation 
and finally energy conservation law takes form 

3 

ZjZj = F{\zi2\^) +F(|zi3|2) +F(|z23p) + const 

• Conservation of angular momentum 

By differentiating of angular momentum 

3 

with respect to f, we get 

Im + /12Z12Z12 + /13Z13Z13 +/23Z23Z23 

Im (^/I2|zi2p+/l3kl3p + /23k23|^) = 0, 

while force / is some real- valued function. 

In general case there are no other conservation laws and the problem is not in- 
tegrable. On the other hand, one may look for some periodical solutions of 
Sys.(1.117) and try to deduce the necessary conditions of periodicity. Importance 



F' = f 
(1.122) 

(1.123) 



(1.124) 



(1.125) 
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of the existence of periodical solutions was pointed out already by Poincare and 
is sometimes even regarded as a definition of integrability - just as opposite case 
for a chaos. 

Theorem 1.4. If fij > 0, (so-called repulsive case), then Sys.(1.117) has 

no periodical solutions. 

► Indeed, in case of periodical solution, magnitude of inertia momentum 

^:=|zi2p+|zi3p+k23p (1.126) 

should have minimums and maximums as sum of distances between three masses. 
On the other hand, 

1 n |2 I |2 I |2\ 
2j^(l^l2| +1^131 +1^231 ) = 

Izi2p+|zi3|^ + |z23|^+/l2|zi2|^+/l3kl3|^+/23|z23|^>0, (1.127) 

which contradicts to the fact that function ^ has to have different signs in the 
points of minimum and maximum. I 

One interesting case - Poincare case - though does not lead to integrable reduc- 
tion of Sys.(1.117), give quite enlightening results and allows to regard this case 
as "almost" integrable. In this case there exists one more conservation law - con- 
servation of inertia momentum 

ki2|^ + ki3p + k23|^ = const (1.128) 

and it is possible due to a special choice of function fj^ = l/|z;fe|^ which allows 
us to reduce initial system to the ODE of the form 

^" = a{^')^ + b{^'f + c^' +d. (1.129) 

in new polar coordinates. This equation describes a geometrical place of points 
on the plane, i.e. some plane curve providing solutions of initial system. The 
curve ^ is of figure-eight form and can not be described by any known algebraic 
curve. On the other hand, it can be approximated with desirable accuracy, for 
instance, by lemniscate 



(1.130) 
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Very comprehensive collection of results and graphics one can find in [26] 

There exists hypothesis that the only periodical solution of Eq.( 1.129) is this 
eight-form curve (not proven). Existence theorem for non-equidistant periodical 
solutions is proven for a wide class of functions / (in variational setting). 

The simplest possible case of periodical solution can be obtained if one of Zi is 
equal to zero (obviously, the problem is reduced to two-body problem, see Ex.6). 
Two more complicated classical integrable cases with periodical solutions for par- 
ticular choice of z-s are known: 

- Lagrange case: |zi2| = jzol = |223|- 

It means that distances between three masses are equal as well as all corre- 
sponding attraction forces, and the masses are moving along a circle. La- 
grange case is also called equidistance case. In this case Sys.(l.l 17) can be 
reduced to ODE y = f{y) and solved in quadratures. 

- Calogero case: all zj are real and fjk= 1 / kj/tl^ for j,k= l,...n 

It means that all masses are moving along a line (in fact, along a real axes) 
and this is generalization of Euler case of three-body problem which after 
appropriate change of variables Sys.(1.117) takes form 

dt^ {xi —X2)^ ~'~ {xi — X3)^ ^ ^ 

d^X2 1 1 



dt'^ {xi-X2)^ {X2-X3)'^ 

dh3 1 1 



+ 7- -TT (1-132) 

(1.133) 



dt^ (xi-X3)3 (X2-X3)3 

Euler 's system was generalized in [23], [24] to 
d^x- 1 

-^ = ^ Witt, ^ = Lj^2 and J=l,...,n. (1.134) 

All n independent conservation laws were found and it was proven that the 
system is integrable. 



1.6 Summary 

In our first Chapter we introduced a notion of differential operator and gave few 
different definitions of integrable differential equation. It was shown that some 
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of them can be equivalent for a given equation and it is reasonable to choose an 
appropriate one depending on what properties of the equation are under the study. 
Some interesting and physically important examples of "almost" integrable sys- 
tems were described. Very intrinsic question on interconnections of conservation 
laws and symmetries was also discussed. 



Simplest possible example of nonlinear equation - famous Riccati equation which 
is ODE of first order with quadratic nonlinearity. We will use this equation in the 
next Chapter for demonstrate many properties of differential equations described 
above. Riccati equation will also be very useful for introduction of some new 
notions like singularities of solutions, integrability tests, etc. 



1.7 Exercises for Chapter 1 

1. Using (1.5) prove that = D2D1 — D1D2 satisfies Leibnitz rule. 

2. Prove that operator 

can be transformed into an operator with constant coefficients by the following 
change of variables: x = e'. (Euler) 

3. Transform equation ut = (p{u)ux into Vt = vvx by appropriate change 
of variables. 

4. Prove that LODE with constant coefficients has conservation laws of the form 

^^^e-^3. = eonst. (1.135) 

5. For potential energy ^ from § 5.2 prove Lagrange- Jacobi identity: 

d% 

— ^y\^ — ^-^2^ — yyi-^ — ^-^3^ — ^>'33— = 

ax\ ay\ 0x2 oy2 0x3 oy3 

2(/l2kl2p + /l3kl3|^+/23k23|^). (1.136) 



6. For particular case zi = and our choice of coordinate system solve Sys.(l . 117) 
explicitly. 



CHAPTER 1. GENERAL NOTIONS AND IDEAS 



Chapter 2 
Riccati equation 



2.1 Introduction 

Riccati equation (RE) 

^^ = a{x)^^ + b{x)^+c{x) (2.1) 

is one of the most simple nonlinear differential equations because it is of first or- 
der and with quadratic nonlinearity. Obviously, this was the reason that as soon 
as Newton invented differential equations, RE was the first one to be investigated 
extensively since the end of the 17th century [39]. In 1726 Riccati considered the 
first order ODE 

Wx = w^ + u{x) (2.2) 

with polynomial in x function u{x). Evidently, the cases degM = 1,2 correspond 
to the Airy and Hermite transcendent functions, respectively. Below we show that 
Hermite transcendent is integrable in quadratures. As to Airy transcendent, it is 
only F-integrable (See Ex.3) though the corresponding equation itself is at the first 
glance a simpler one. 

Thus, new transcendents were introduced as solutions of the first order ODE with 
the quadratic nonlinearity, i.e. as solutions of REs. Some classes of REs are 
known to have general solutions, for instance: 

y+ay'^ = bx'' (2.3) 

where all a,b, a are constant in respect to x. D. Bernoulli discovered (1724-25) 
that this RE is integrable in elementary functions if a = —2 or a = —4k{2k — 
1), k= 1,2,3, Below some general results about RE are presented which 
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make it widely usable for numerous applications in different branches of physics 
and mathematics. 



2.2 General solution of RE 

In order to show how to solve (2. 1) in general form, let us regard two cases. 
2.2.1 a{x) = 

In case a{x) = 0, RE takes particular form 

(l)x^b{x)(j)+c{x), (2.4) 

i.e. it is a first-order LODE and its general solution can be expressed in quadra- 
tures. As a first step, one has to find a solution z{x) of its homogeneous part (see 
Ex.1), i.e. 

z{x): zx^b{x)z. (2.5) 

In order to find general solution of Eq.(2.4) let us introduce new variable 0(x) = 
(^{x)/z{x), i.e. z{x)^{x) — (^{x). Then 

(z(x)0(x))x = &(x)z(x)0(x)+c(x), i.e. z{x)^{x)x = c{x), (2.6) 

and it gives us general solution of Eq.(2.4) in quadratures 

(x) = z(x) (x) = z(x) ( / 44^-^ + const) . (2.7) 

J z[x) 

This method is called method of variation of constants and can be easily gener- 
alized for a system of first-order LODEs 

y ^A{x)y + f{x). (2.8) 

Naturally, for the system of n equations we need to know n particular solutions 
of the corresponding homogeneous system in order to use method of variation of 
constants. And this is exactly the bottle-neck of the procedure - in distinction with 
first-order LODEs which are all integrable in quadratures, already second-order 
LODEs are not. 
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2.2.2 a{x)^0 

In this case one known particular solution of a RE allows to construct its 
general solution. Indeed, suppose that <pi is a particular solution of Eq.(2. 1), then 

c=(pi^x-a(pi-b(pi (2.9) 

and substitution ^ —y + (pi annihilates free term c yielding to an equation 

yx = a/ + by (2.10) 

with b — b + 2a(pi. After re-writing Eq.(2. 10) as 

yx .b 

-iy=a + - (2.11) 

y y 

and making an obvious change of variables 0i = 1 /y, we get a particular case of 
RE 

0l,x + ^0l-« = O (2.12) 
and its general solution is written out explicitly in the previous subsection. 



Example 2.1 

As an important illustrative example leading to many applications in mathematical 
physics, let us regard a particular RE in a form 

yx+y^^x^ + a. (2.13) 

For a = 1, particular solution can be taken asy = x and general solution obtained 
as above yields to 

J e~^^ dx + const 

i.e. in case (2.13) is integrable in quadratures. Indefinite integral / e~^dx though 
not expressed in elementary functions, plays important role in many areas from 
probability theory till quantum mechanics. 



For arbitrary a, Eq.(2.13) possess remarkable property, namely, after an elemen- 
tary fraction-rational transformation 



y = x-\ 

y+x 



(2.15) 
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it takes form 

y^+f = x^ + a, a = a + 2, (2.16) 

i.e. form of original Eq.(2.13) did not change while its rhs increased by 2. In 
particular, after this transformation Eq.(2.13) with a = 1 takes form 

y, + f=x^ + 3 (2.17) 

and since y = x is a particular solution of (2.13), then y = x+l/xisa particular 
solution of the last equation. It means that for any 

a = 2A:+l, A: = 0,l,2,... (2.18) 

general solution of Eq.(2.13) can be found in quadratures as it was done for the 
case a = 1. 



In fact, it means that Eq. (2. 1 3) is form-invariant under the transformations (2.15). 
Further we are going to show that general RE possess similar property as well. 



2.2.3 Transformation group 

Let us check that general fraction-rational change of variables 

rW0 + 5(x) ^""'^^ 

transforms one Riccati equation into the another one similar to Example 2. 1. No- 
tice that (2.19) constitutes group of transformations generated by 

^, a{x)(j>, + j8(x), (2.20) 

thus only actions of generators have to be checked: 

• = l/(^ transforms (2. 1) into 

^x + c{x)^^ + b{x)^+a{x) =0, (2.21) 

• = a(jc)0 transforms (2.1) into 

$x - - [b{x) + (loga(^))j0 - a{x)c{x) = 0, (2.22) 
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• = + j8 (x) transfonns (2.1) into 

4 - a{x)f + [2/3 {x)a{x) -b{x)]^-c = 0, (2.23) 

where 

c = a{x)P^{x)-b{x)P{x)+c{x) + I5{x)x. (2.24) 

Thus, having one solution of a some Riccati equation we can get immediately 
general solutions of the whole family of REs obtained from the original one under 
the action of transformation group (2.19). 



It is interesting to notice that for Riccati equation knowing any three solutions 
01) 02, 03 we can construct all other solutions using a very simple formula called 
cross-ratio: 



0-02 03-02 

with an arbitrary constant A, where choice of A defines a solution, 
verify this formula let us notice that system of equations 

'0 = a{x)^^ + b{x)^ +c{x) 

01 = a{x)<j)f + b{x)<j)i+c{x) 

02 = a(x)02 + &(x)02 + c{x) 

^ 03 = a(x)0| + b{x)<l)3 + c{x) 



(2.25) 
In order to 



(2.26) 



= 



is consistent if 

02 ^ 

01 02 01 

02 02 02 
.03 03^ 03 

and direct calculation shows that this condition is equivalent to 

J_.0-01 03-01- 

dx 



(2.27) 



-)=0. 



(2.28) 



■0-02 03-02' 

As it was shown, REs are not invariant under the action of (2. 1 9) while (2. 19) con- 
serves the form of equations but not form of the coefficients. On the other hand, 
it is possible to construct new differential equations related to a given RE which 
will be invariant with respect to transformation group (2.19) (see next section). 



At the end of this section we consider a very interesting example [45] show- 
ing connection of Eq.(2.28) with first integrals for generalization of one of Ko- 
valevskii problems [48]. 
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Adler's example. System of equations 

n 

yj,x + ^yj = syj, s^Y^yj^ 7=l,2,...n (2.29) 

was studied by Kovalevskii in case n = 3 and it was shown that there exist two 
quadratic first integrals 

Fi = (yi - y2)y3, Pi = (yi - y3)yi (2.30) 

and therefore Kovalevskii problem is integrable in quadratures. 

In case n>4 the use of (2.28) gives us immediately following some first integrals 

yi-yi yk-yi 



yt-yjyk-yj 

i.e. Sys.(2.29) has nontrivial first integrals for arbitrary n. 



(2.31) 



It is interesting that for this example solution of Sys.(2.29) is easier to construct 
without using its first integrals. Indeed, each equation of this system is a Riccati 
equation if a is regarded as given, substitution = ^i^/2^i gives 

0/,xc = 50!,xc, 0!,xc = a{x)+Ci, S = Qxx/ % (2.32) 

and equation for a has form 

^^^(^ + ... + ^). (2.33) 

Ux 2 a — c\ a — Cfi 

After integration 

(2^ = const((3 — ci)...((2 — c„), (2.34) 

i.e. problem is integrable in quadratures (more precisely, in hyper-elliptic func- 
tions). 

In fact, one more generalization of Kovalevskii problem can be treated along the 
same lines - case when function s is not sum of yj but some arbitrary function 
s = s(xi, ...,x„). Then equation on a takes form 

— =axs{ ) (2.35) 

Ux a — c\ a — Cti 

which concludes Adler's example. 
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2.2.4 Singularities of solutions 

All the properties of Riccati equations which have been studied till now, are in the 
frame of local theory of differential equations. We just ignored possible existence 
of singularities of solutions regarding all its properties locally, in a neighborhood 
of a point. On the other hand, in order to study analytical properties of solutions, 
one needs to know character of singularities, behavior of solutions at infinity, etc. 

One can distinguish between two main types of singularities - singularities, not 
depending on initial conditions (they are called fixed) and depending on initial 
conditions (they are called movable). Simplest possible singularity is a pole, and 
that was the reason why first attempt of classification of the ordinary nonlinear 
differential equations of the first and second order, suggested by Painleve, used 
this type of singularities as criterium. Namely, list of all equations was written out, 
having only poles as movable singularities (see example of PI in Chapter 1), and 
nice analytic properties of their solutions have been found. It turned out that, in 
particular, Painleve equations describe self-similar solutions of solitonic equations 
(i.e. equations in partial derivatives): P2 corresponds to KdV (Korteweg-de Vries 
equation), P4 corresponds to NLS (nonlinear Schrodinger equation) and so on. 

Using cross-ratio formula (2.25), it is easy to demonstrate for a Riccati equation 
that all singularities of the solution <j), with an exception of singularities of par- 
ticular solutions 01, 02, 03, are movable poles described as following: 

03 = YZI^*^2-A0i) (2.36) 

where A is a parameter defining the solution 0. Let us construct a solution with 
poles for Eq.(2.13) from Example 2. 1. We take a solution in a form 

1 9 

+ao + ai{x + £)+a2{x + £Y + ... (2.37) 

x+e 

with indefinite coefficients a;, substitute it into (2.13) and make equal terms cor- 
responding to the same power of (x + e). The final system of equations takes 
form 

ao = 0, 

3ai — a — — 0, 
4a2 + 2e = 0, 
< 5a3-l+flf = 0, (2.38) 
6a4 + 2aia2 = 0, 
las + 2flifl3 +0-2 = 
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and in particular for a = 3, e = the coefficients are 

ai — I, a2 — ut, — ... = (2.39) 
which corresponds to the solution 

y = x+- (2.40) 

X 

which was found already in Example 2.1. 

This way we have also learned that each pole of solutions have order 1 . In general 
case, it is possible to prove that series (4.7) converges for arbitrary pair (e, a) 
using the connection of RE with the theory of linear equations (see next section). 
In particular for fixed complex a, it means that for any point xq = — e there exist 
the only solution of (2.13) with a pole in this point. 

As to nonlinear first order differential equations (with non -quadratic nonlinearity), 
they have more complicated singularities. For instance, in a simple example 

yx = y^ + i (2.41) 

if looking for a solution of the form y = cu^ + .... one gets immediately 

afcc*^-^ + ... = aV^+... ^ k-\ = ?>k ^2k=-l (2.42) 

which implies that singularity here is a branch point, not a pole (also see [40]). It 
make RE also very important while studying degenerations of Painleve transcen- 
dents. For instance, (2.13) describes particular solutions of P4 (for more details 
see Appendix). 



2.3 Differential equations related to RE 
2.3.1 Linear equations of second order 

One of the most spectacular properties of RE is that its theory is in fact equivalent 
to the theory of second order homogeneous LODEs 

¥xx = b{x) ¥x + cix) Xjf (2.43) 

because it can easily be shown that these equations can be transformed into Ric- 
cati form and viceversa. Of course, this statement is only valid if Eq.(2.1) has 
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non-zero coefficient a{x), a{x) ^ 0. 

► Indeed, let us regard second-order homogenous LODE (2.43) and make change 
of variables ^ 

= ^, then 0^=^^, (2.44) 

which implies 

^ = 0.+ ^ = 0. + 0' (2.45) 
and after substituting the results above into initial LODE, it takes form 

(l)x = (l)^ + b{x)(l)+c{x). (2.46) 
which is particular case of RE. I 

< On the other hand, let us regard general RE 

(l>x = a{x)(j>^ + b{x)(l>+c{x) (2.47) 

and suppose that a{x) is not = while condition of a{x) =0 transforms RE into 
first order linear ODE which can be solved in quadratures analogously to Thomas 
equation (see Chapter 1). Now, following change of variables 

= -^ (2.48) 

transforms RE into 

-J^ + ^(^)\^^=a(x)(^Y-^-^^+c(x) (2 49) 
a(x)\j/ a(x)\\j// a(x)2 Xj/ \a{x)\i/J a{x) Xjf 

and it can finally be reduced to 

a{x)\i/xx- [a{x)x + a{x)b{x)]\i/x + c{x)a{xf\i/ = (2-50) 
which is second order homogeneous LODE. I 

Now, analog of the result of Section 2.2 for second order equations can be proved. 

Proposition 2.1 Using one solution of a second order homogeneous LODE, we 
can construct general solution as well. 
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► First of all, let us prove that without loss of generality we can put b{x) =0 in 
¥xx + b{x) Yx + c{x) i/A^ = 0. Indeed, change of variables 



and finally 



¥x^{¥x-\byif)e~y'^^^^' (2.51) 



1 1 

il/xx + cil/ = 0, c = c--b^--bx. 



(2.52) 



Now, if we know one particular solution Yi of Eq.(2.52), then it follows from the 
considerations above that RE 



(2.53) 



has a solution 0i = ij/i^/ ^i. The change of variables = — 0i annihilates the 
coefficient c{x): 



(0 + 0i)'+(0 + 0i)^) + c(x)=O 0;, + 02 + 20i0 = O 



(2.54) 
(2.55) 



i.e. we reduced our RE to the particular case Eq.(2.4) which is integrable in 
quadratures. Particular solution z = 1/0 of homogeneous part of Eq.(2.55) can 
be found from 



Zx = 2z-r^ as z = V^T 



and Eq.(2.7) yields to 



¥2{x) = m 



dx 
xfffix) 



(2.56) 



(2.57) 



Obviously, two solutions and i//2 are linearly independent since Wronskian 
< \j/i,\j/2> is non- vanishing (see Ex.2): 



< ¥h¥2 >■= 



¥i ¥2 
¥i ¥2 



= V'iV^2-V^2V^i = 1^0. 



(2.58) 



Thus their linear combination gives general solution of Eq.(2.43).l 



Proposition 2.2 Wronskian < i/^i , 1/^2 > is constant iff xj/i and 1/^2 are solutions of 

¥xx = c{x)\i/. (2.59) 



► Indeed, if Yi and Y2 are solutions, then 

(V^i v4 - ¥2¥'i)' = ¥1 ¥2 -¥2¥i= cix) (yi ¥2 -¥i¥2) = 
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■< if Wronskian of two functions Xj/i and is a constant, 

V^iv4-V^2V^l = const ^ V^iv4'-V^2V'i =0 

" // 

Y2 ¥1 ' 



(2.61) 



Conservation of the Wronskian is one of the most important characteristics of 
second order differential equations and will be used further for construction of 
modified Schwarzian equation. 

To illustrate procedure described in Proposition 2. 1, let us take Hermite equation 



cOxx - ^xcOx + 2A(0 = 0. 
Change of variables z = (£>xl^ yields to 

= — -2^5 Zjc + z^-2xz + 2/l = 

(£> 

and with 3; = z — x we get finally 

Jx+J^ =x2-2A - 1, 



(2.62) 



(2.63) 



(2.64) 



i.e. we got the equation studied in Example 2.1 with a = —2X — 1. It means that 
all solutions of Hermite equation with positive integer A , A = n, n EN can easily 
be found while for negative integer X one needs change of variables inverse to 
(2.15): 



y= -x+ ■ 



r 



y-x 

It gives us Hermite polynomials 



{y-x){y^-x) = y, 7=d-l, a = a + 2. 



(2.65) 



X 


= 0, 


y = 


-X, (0=1 


X 


= 1, 


y = 


— x+i, (0 = 2x 


A 

< 


= 2, 


y = 





(2.66) 



A = „, 3; = -x+f, a) = //„(x) = (-l)V'^(e--0 
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as solutions. 



Notice that the same change of variables 

= ^ (2.67) 

which linearized original RE, was also used for linearization of Thomas equa- 
tion and Burgers equation in Chapter 1. This change of variables is called log- 
derivative of function i// or Dlog(i//) and plays important role in many different 
aspects of integrability theory, for instance, when solving factorization problem. 



Theorem 2.3. Linear ordinary differential operator L of order n could be factor- 
ized with factor of first order , i.e. L = Mo[dx — a) for some operator M, iff 

a=^, where v^eKer(L). (2.68) 
>■ L = Mo{dx — a), a= Vj^/ l/A implies {dx — a)\\f = Q,i.Q. veKer(L). I 



< Suppose that = 1 is an element of the Ker(L), i.e. y/i G Ker(L). It leads to 
a = and operator L has zero free term and is therefore divisible by dx- 



If constant function i^i = 1 does not belong to the kernel of initial operator, fol- 
lowing change of variables 

W^— (2.69) 



lead us to a new operator 
which has a constant as a particular solution y/i for f = Yi 



L = f-'Lof (2.70) 



Remarli. Operators L and L given by (4. 1 1), are called equivalent operators and 
their properties will be studied in detailed in the next Chapter. 



Notice that Theorem 2.3 is analogous to the Bezout's theorem on divisibility cri- 

terium of a polynomial: A polynomial P{z) = is divisible on the linear factor, 
P{z) = P\ (z) {z — a), iff (2 is a root of a given polynomial, i.e. P{a) = 0. Thus, in 
fact this classical theorem constructs one to one correspondence between factor- 
izability and solvability of L{^) — 0. 
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The factorization of differential operators is in itself a very interesting problem 
which we are going to discuss in details in Chapter 3. Here we will only regard 
one very simple example - LODO with constant coefficients 

d"\i/ d"^^\i/ 

In this case each root A, of a characteristic polynomial 

X" + aiX"'^ + ... + an = (2.72) 
generates a corresponding first order factor with 

?ii = — (2.73) 

ant it yields to 

\l/x = ^iW W = Cie^'' (2.74) 



and finally 



This formula allows us to construct general solution for L{\j/) = 0, i.e. for Xj/ e 
Ker(L), of the form 

V/-£c,e^'^ (2.76) 
in the case of all distinct roots of characteristic polynomial. 

In case of double roots X,k with multiplicity rrik it can be shown that 

llf = Y,Pk{x)e^'' (2.77) 

where degree of a polynomial i\(jc) depends on the multiplicity of a root, degPi^{x) < 
rrik - 1 (cf. Ex.4) 

2.3.2 Schwarzian type equations 



Let us regard again second-order LODE 

\ifxx + b{x)\if:c + c{x)\if=Q 



(2.78) 
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and suppose we have two solutions Xj/i , \j/2 of (2.78). Let us introduce new function 
9 = V^i/ 1/^2, then 

<Px = 2 ' 

^xc = b[x) 2 ^2 3 V^2jc, (2.80) 

which yields 

^ = -b(x)-2^ (2.81) 
and substituting = ^ = (log 1/^2 )x into (2.1) related to (2.78) we get finally 



3 f (pxx\^ 1 (Pxxx 



c{x). (2.82) 



Left hand of (2.82) is called Schwarz derivative or just Schwarzian and is invari- 
ant in respect to transformation group (2. 19) with constant coefficients a, j8 , 7, 5: 



It is sufficient to check only two cases: 



(2.83) 



^ = - and 9 = «()!) + j8. (2.84) 

which can be done directly. 

Remark. The schwarzian equation (2.82) with 

^W = Er^+Ir^ (2.85) 

plays major role in the theory of conformal transformations of polygons [46]. In 
this case (p{x) maps polygon with vertexes ai , . . . fl„ on the complex x— plane into 
unit circle. 

Notice now that for linear second order ODE in the canonical form 

¥xx = v{x)\j/ (2.86) 

we have 

WuWi-WixWi const 

^>X— 9 — T- 

¥2 V2 
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and it gives rise the question what should be equation for the A = i/^iv^2 in- 
stead (Px = Wi^- This question leads us to equation which generalize the classical 
Schwarzian equation (2.82). 

Theorem 2.4 Let '^f\ , i/^2 are two linear independent solutions of Eq.(2.86) Then 
functions 

Wh ¥i¥2 (2.88) 
constitute a basis in the solution space of the following third order equation: 

= 4v{x)Ax + 2vxix)A. (2.89) 

Moreover function cp — Yi/Wi satisfies Schwarzian equation (2.82) as well as 
A = 1/^1 1/^2 is a solution of the equation 

4v{x)A^ +aI - lAAxx = (2.90) 
where w = \if\xW2 -W2xWi- 



► Using notations 



M = V/l^ A2 = V/|, A3 = v/1 VA2, (2.91) 



we can compute Wronskian W of these three functions 

W =< Ai,A2,A3 >= {\l/i\l/2,x-¥2Vi,xf =< ¥h¥2 >^ (2.92) 

and use Proposition 3.2 to demonstrate that 

W = const ^ 0, 

i.e. functions A, are linearly independent. 
After introducing notations 



it is easy to obtain 
which yields to 



w =< v^i, V^2 > and fj = ^ (2.93) 

Wj 



W At. r 

l-=/2-/i, ^=/2+/i (2.94) 

A3 A3 

At, x — w Al V + w 



Substitution of these fj into 

fj,x + fj=v{x) (2.96) 
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gives rise Eq.(2.90) with A = A3 and it's differentiation of with respect to x gives 
Eq.(2.89) and it easy to see that equations (2.90) and (2.89) are equivalent. 

Analogous reasoning shows that Ai, A2 are also solutions of Eq.(2.89) but in this 
case w = in Eq.(2.90). 

At last, rewrite Eq.(2.90) as 

, , A? 2Axx 
and introduce notation H =1 /A, then 

and compare this equation with Schwarzian Eq.(2.82) one can see immediately 
that it corresponds to Eq.(2.98) with w — and H — (p^.t 



Substitution 9 = V^i/ (^2 allowed us to get invariant form of the initial Eq.(2.78), 
and this is the reason why we call it modified Schwarzian equation. This form 
of Schwarzian equation turns out to be useful for a construction of approximate 
solutions of Riccati equations with parameter (see next section), is called modi- 
fied Scliwarzian equation. 



Notice that after the substitution a — e , rhs of modified Schwarzian equation, 
i.e. modified Schwarzian derivative, Dmod, takes a very simple form 

Dmod{a) := \^-\—= b„ + bl (2.99) 
which is in a sense similar to Dlog. Indeed, for l/A = e^, 

Dlog{xif) = ^ = (p, = e-^-^e^, (2.100) 
X]/ ax 

while 

Dmod{e^f) = e^'^-^e'P. (2.101) 

At the end of this section let us stress the following basic fact: we have shown 
that from some very logical point of view, first order nonlinear Riccati equation, 
second order linear equation and third order nonlinear Schwarzian equation are 
equivalent! It gives us freedom to choose the form of equation which is most 
adequate for specific problem to be solved. 
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2.4 Series in spectral parameter 

In our previous sections we have studied Riccati equation and its modifications 
as classical ordinary differential equations, with one independent variable. But 
many important applications of second order differential equations consist some 
additional parameter A, for instance one of the most significant equations of one- 
dimensional quantum mechanics takes one of two forms 

= {X+u)\j/ (2.102) 
Vxx = (AVmiA + M2)v^ (2.103) 

where Eq.(2.102) is called Schrodinger equation and Eq.(2.103) can be consid- 
ered as modified Dirac equation in quantum mechanics while in applications to 
solitonic hierarchies it is called Zakharov-Shabat equation. 

Notice that Schrodinger equation, yet not (2.103), can be readily rewritten as the 
equation for eigenfunctions of the operator L = D^ — u: (That fact will be used in 
next chapter). 

\l/xx= {X + u)\i/ <^ L\i/ = Xxj/ (2.104) 

Correspondingly, we call L = — u as Schrodinger operator and A is called 
spectral parameter. In order to unify terminology it is convenient to name the 
whole coefficient before y/ as generalized potential allowing it sometimes to be 
a polynomial in A of any finite degree. Thus, coming back to Eq.(2.86), the gen- 
eralized potential here is just the function c{x) but in Example 2.1 the parametric 
dependence has been important and as matter of fact this example corresponds to 
Schrodinger equation (2.102) with the harmonic potential u{x) = x^. 

Now, with the equation having a parameter, problem of its integrability became, 
of course, more complicated and different approaches can be used to solve it. If 
we are interested in a solution for all possible values of a parameter A, asymptotic 
solution presented by a formal series can always be obtained (section 4.1) while 
for some specific exact solutions can be constructed (sections 4.2, 4.3) in a case of 
truncated series. It becomes possible while existence of a parameter gives us one 
more degree of freedom to play with. Cf. with Example 4.2 where exact solution 
has been obtained also as a series and its convergence resulted from the main 
theorem of the theory of differential equations on solvability of Cauchy problem. 
On the other hand, this solution is valid only for some restricted set of parameter' 
values, namely for integer odd a. 



46 



CHAPTER 2. Rice ATI EQUATION 



2.4.1 RE with a parameter A 

Let us show first that the RE with a parameter A corresponding to Eq.(2.103), 
namely 

/x + / = A2 + miA+M2, with f^D,\og{\j/), (2.105) 
has a solution being represented as a formal series. 

Lemma 2.4 Eq.(2.105) has a solution 

/ = A+/o + ^ + ... (2.106) 

where coefficients fj are differential polynomials in ui and uj. 

► After direct substituting the series (2.106) into the equation for / and making 
equal corresponding coefficients in front of the same powers of A, we get 

2/o = Ml 

2/l+/o,x + /o =M2 

< 2/2 + /i,x + 2/o/i=0 (2.107) 
2/3+/2,x + 2/o/2+/f = 



and therefore, coefficients of (2.106) are differential polynomials of potentials ui 
and M2- ■ 

Notice that taking a series 

g=-A+go + |- + || + || - (2.108) 

as a form of solution , we will get a different system of equations on its coefficients 
8i- 

-Igo = Ml 

-2gl +gQ,x + gl = M2 

< -2g2+gi,. + 2gogi=0 (2.109) 

-2^3 + g2,x + 2gog2 +gj = 
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Solution of Eq.(2.105) constructed in Lemma 4.1 yields to the solution of original 
Zakharov-Shabat equation (2.103) of the form 

Y,ix,X) = .//(-'^)^- = e'^iUx) + ^ + ^ + + ....) (2.110) 
and analogously, the second solution is 

V.2(x, A) = ./^M)'^- = e-'^mx) + ^ + ^ + ^ + ....) (2.111) 

In fact, it can be proven that Wronskian < i/Aj , i//2 > is a power series on X (see 
Ex.8) with constant coefficients. Notice that existence of these two solutions is 
not enough to construct general solution of initial Eq.(2.103) because linear com- 
bination of these formal series is not defined, also convergence problem has to 
be considered. On the other hand, existence of Wronskian in a convenient form 
allows us to construct family of potentials giving convergent series for (2.110) 
and (2.1 1 1). We demonstrate it at the more simple example, namely Schrodinger 
equation (2.102). 

Let us regard Schrodinger equation (2. 102), its RE has form 

fx + f = ?i + u, with f = Dlogixif), (2.112) 
and it can be regarded as particular case of (2.103), i.e. the series for / yields to 

f = k + fo + ^ + ..., A=^, (2.113) 
k 

and g{x,k) = f{x, —k). We see that in case of (2.102) there exists a simple way to 
calculate function g knowing function / and it allows us to construct two solutions 
of Schrodinger equation (2. 102): 

Xlfi{x,k) = .//(^'^)^- = e'\l + iiM + ^ + ^ + ....) (2.114) 

and 

\l/2{x,k) = \i/i{x,-k). (2.115) 
Substitution of say xj/i into (2.102) gives a recurrent relation between coefficients 

Cj+i,x = ^{uCj - Q,xx), Co = 1- (2.116) 

In particular, 

w = 2Ci,x (2.117) 
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which means that in order to compute potential u it is enough to know only one 
coefficient of the formal series ! Below we demonstrate how this recurrent re- 
lation helps us to define potentials corresponding to a given solution. 



Example 2.5. Let us regard truncated series corresponding to the solutions of 
(2.102) 

y,l=/-(l + ^), y,2 = e-^"(l-|), (2.118) 

then due to (2.116) 

M = 2Ci,x, Ci,xc = 2Ci,xCi (2.119) 
and Wronskian W of these two functions has form 

W^<yjfi, W2 >= -2^+ ^(Cf - Ci,x)- (2.120) 

Notice that 

(Ci'-Ci,x)x = (2.121) 

and it means that W does not depend onx,W = W{k). Introducing notation ki for 
a zero of the Wronskian, W{ki) = 0, it is easy to see that 

Ci-Ci,x = kj (2.122) 

which implies that Xj/l ¥2 are solutions of (2.102) with 

2ki 

^i=ki- ^^^_2fc^(,_,„) = -^itanh^i(x-xo) (2.123) 

and potential 

u = -2 = ^ (2 124) 

(e*:i(x-^o) + e-fei(^-^o))2 cosh^(A:i(x-xo))' 

where xq is a constant of integration. 



It is important to understand here that general solution of Schrodinger equation 
(2.102) can be now found as a linear combination of xj/i xj/i for all values of 
a parameter k =k^ with exception of two special cases: k = and k = ki which 
implies functions Xj/i Xj/i are linearly dependent in these points. 
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At the Fig. 1 graph of potential u is shown and it is easy to see that magnitude 
of the potential in the point of extremum is defined by zeros of the Wronskian 
W. At the end of this Chapter it will be shown that this potential represents a 
solitonic solution of stationary KdV equation, i.e. solution of a Riccati equation 
generates solitons! 

2.4.2 Soliton-like potentials 

It this section we regard only Schrodinger equation (2.102) and demonstrate that 

generalization of the Example 4.2 allows us to describe a very important special 
class of potentials having solutions in a form of truncated series. 

Definition 2.5 Smooth real- valued function u{x) such that 

u{x)^0 for x^±oo, (2.125) 

is called transparent potential if there exist solutions of Schrodinger equation 
(2. 102) in a form of truncated series (2. 1 14) 

Another name for a transparent potential is soliton-like or solitonic potential due 
to many reasons. The simplest of them is just its form which is a bell-like one 
and "wave" of this form was called a soliton by [49] and this notion became one 
of the most important in the modem nonlinear physics, in particular while many 
nonlinear equations have solitonic solutions. 

Notice that truncated series for and \]/2 can be regarded as polynomials in k 

of some degree A'^ multiplied by exponent e^^^ (in Example 4.2 we had A'^ — 1). 
In particular, it means that Wronskian W =< i//2 > is odd function, W{—k) = 
—W{k), vanishing at A; = and also it is a polynomial in k of degree 2N+ 1: 

N 

W{k)^-2kY[{k^-kj). (2.126) 
1 

As in Example 2.4., functions and 1/^2 are linearly dependent at the points kj, 
i.e. 

xif2{x,kj) = {-iy+%xifi{x,kj), 7=l,2,...,Ar (2.127) 
with some constant proportionality coefficients which we write down as ( — 1 ) ^+ ^ Ay . 



In order to find the coefficients of the polynomial 

e~'^\l/l{x,k) =J^ + ail^'~'^ + ...+aN 



(2.128) 
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it is enough to write out explicitly and \j/2 in roots k = kj of polynomial (2.126) 
and substitute these in equations (2.127). Thus we find that equations (2.127) 
are equivalent to the following system of N linear equations on A'^ unknowns 
ai, . . . ,a;v : 



+ Ei=0 



- k2 



+ ^+£2^ + ---+l=0 



l+^3| + | + -+^3 = 



(2.129) 



where following notations has been used: 



e } —e 1 



: = tanh Xj , Ty = kjx + j8y , Aj = e^^j . 



(2.130) 



e'^j+e '^i 

Lemma 2.6. Let numbers kj are ordered in following way 

(2.131) 

and the values of variables Ej in coefficients of system (2.129) obey conditions 

-1<£; <1, V (2.132) 

Then the determinant A(£'i, . . . ,£'7v) of system (2.129) of linear equations on 
unknowns ai, . . . ,aiv does not vanish and satisfies inequality as follows 



A{Ei,...,En) > 



1 



n 



ki ■ ■ -kf^ \ki 



1 1 



(2.133) 



► Obviously determinant A(£'i , . . . , Ej^) is linear functions in Ej and as such takes 
minimal value on the boundary of the set (2.132) at some sequence Ej — Ej, Ej = 
il.Thus 



/ ^1 ^ £1^1 



A>det 



E2k 



V 



-1 
\ 

^3 



,-2 



£3^; 



1 

-3 
^3 



... \ 



E2k 



def . 
= Ag. 



(2.134) 



Since = 1 latter determinant can be computed in the closed form. Indeed 



£264 ■ 



ki- ■ - kisf 



-det 



/ 1 (£i^i)-i (eifci)-2 ... \ 

1 (£2^2)"^ (£2^2)"^ ... 
1 (£3/^3)-! (£3^:3)-^ ... 



V 



(2.135) 
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Thus, use formula for Vandermond determinant we obtain 




(2.136) 
(2.137) 



Theorem 2.7. Suppose we have two sets of real numbers 

{kj}, {/3,}, j = l,2,...,yv, kj,l5jeR (2.138) 

and numbers kj are ordered as it shown in (2.131). Then conditions (2.132) are 
fulfilled and defined by the linear system (2.129) the functions \i/i{x,k), Xj/iix^k) 

\l/i{x,k) = e'^{k^ + aik^'^ + ...+aN), VA2(^,A:) = (-l)Vi(^, -^) (2.139) 
are such that 

i) Wronskian W{k) v/i, 1/^2 >= -2kUiik^ - kj) 

ii) Schrodinger equation (2.102) is satisfied with transparent potential u such that 

-2kl <u<0. (2.140) 



Proof of i). 

In order to compute the Wronskian W of functions (2.139), notice first that W is a 
polynomial with leading term —2^^^+^. Condition of proportionality (2.127) for 
functions Xj/i {x, kj) , Y2{x, kj) provides that kj are zeros of the Wronskian and that 
W is an odd function on k, i.e. (2.126) is proven. 

Formal proof of ii) is a bit technical and we postpone it up to next chapter. In any 
case in order to prove that functions (2.139) satisfy Schrodinger equation 

^l^^k^ + u, u = 2ai^ (2.141) 
it is sufficient to verify equations (2. 116) 



= -^^iuaj-Uj^xx), «o = 1 



(2.142) 
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for solution ai,...,aN of the linear system (2.129). In order to illustrate how to 
use this theorem for construction of exact solutions with transparent potentials let 
us address two cases: A^^ = 1 and N = 2. 

Example 2.8. In case N = I explicit form of functions 

Xl/i = e'°'ik + ai), \if2 = e-'°'{k-ax) (2.143) 

allows us to find a\ immediately: 

a\ = -kiEi = -ki tanhyi = -ki tanh(A:ix + j8i) (2.144) 

which coincides with formula for a solution of the same equation obtained in Ex- 
ample 2.5. Notice that using this approach we have found solution of Schrodinger 
equation by pure algebraic means while in Example 4.2 we had to solve Riccati 
equation in order to compute coefficients of the corresponding truncated series. 

The system (2.129) for case N = 2 takes form 

kiai+Eia2 + klEi^0 {kl-kj)Ei 

^ ai = ; — (z.i4D; 

£■2^:20 1 +a2 + k^ = ki- k2E\ E2 

which yields 

a\ = Dx log {{k2 - ki ) cosh(Ti + T2) + (^2 + ki ) cosh(Ti - T2) ) (2. 146) 
and corresponding potential u = 2ai^x has explicit form 

u = 2Dl\og {{k2 - cosh(Ti + T2) + (/:2 + cosh(Ti - T2)) (2.147) 

where 

x^±oo =^ ai^±{ki+k2), (2.148) 
i.e. M is a smooth function such that 

M^O for x^±oo. (2.149) 

Definition 2.9. The point A G C is called eigenvalue of Schrodinger operator 
L = D^ — u if there exist solution (p eigenfunction of equation (2.102) such that 



q){x)^0 for x^±oo. 



(2.150) 
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Remark 2.10 In spectral theory (see classical textbook [41]) it is proven that for 
smooth real-valued potentials u{x) vanishing rapidly as jc — )■ ±oo the set of eigen- 
values of Schrodinger operator L = u is finite and they are all real positive 
numbers X =^j, j = l-,...-,N. Moreover, in the case of standard numeration 

L(pj = Xj(pj, Ai > A2, . . . , Aat > (2. 15 1) 

the eigenfunction has 7 — 1 change of sign on the x— axis. 

In application to transparent potentials from Theorem 4.6 we find that 

±00 ^ Ej^ ±1 (2.152) 

and, therefore, in these limits 

\\fi{k,x)e-'''' ^A^{k)ioix^±oo (2.153) 

where A* [k) are polynomials in k of A^— th degree. Come back to Lemma 4.5 and 
the linear system of equations (2.129) for oi, . . . ,fl;v one finds that 

A' N 

A-{k) = \[{k + kj) 2ind.A-^{k)=Y[{k-kj) (2.154) 
1 1 

It is easy to see now that for the our case the eigenvalues Xj = kj correspond to 
zeros of Wronskian (2.126) and eigenfunctions cpj are just \ifi{kj,x) with positive 
kj numbered as in Eq.(2.131). Indeed, for any positive k the solutions v^i(^,jc) 
and 1/^2 (fc,x) vanish when x — > ^00, respectively, and are proportional each other 
atA; = A;^byEq.(2.127). 

To proceed further and apply more deep results of the spectral theory it is crucial 
to introduce roots variables i.e. zeros k = rj{x), 7 = 1, ... of the polynomial 

e~'^\l/i{k,x): 

N 

\l/iik,x) = e'^Yl(k-rj(x)). (2.155) 
1 

Consider the function \]/i{k,x) upon the interval —00 < x < T one should notice 
that the positive roots rj{T) are eigenvalues of boundary value problem for Eq. 
(2.102) on interval —00 < x < T with zero boundary conditions. Similarly, the 
negative roots rj{T) are eigenvalues of boundary value problem on interval T < 
X < 0° < X < T for the same equation. Apply the spectral theory one can prove 
that roots variables rj{T) increase monotone (Cf Exercise 8) along with T and 



T^±oo ^ rj^±kj. 



(2.156) 



54 



CHAPTER 2. Rice ATI EQUATION 



Furthermore, the 2N—th degree polynomial in k 

N 

^=xi/l{k,x)xi/2{k,x)=ll{k^-rj) (2.157) 

1 

due Eq.(2.127) change the sign on ends of intervals 

[kj,kj+i], k\>k2> ... > > (2.158) 

and it imply that 

k]>r]>kl>...>kl>rl>0 (2.159) 

Substitute the expansion of the polynomial O in A = /c^ into Eq.(2.90) from The- 
orem 3.4 and equate coefficients on one can get the important formula 

A' 

u{x) = 2Y^{r]-l^) (2.160) 
1 

of transparent potentials in terms of roots variables. Thus reality conditions and 
inequalities (2.159) yields the estimate (2.140) formulated in Theorem 4.5. 

Summing up we notice that the fraction 

is totally defined by it values at the zeros k = knOf the Wronskian and due (2.127) 
we have 

login^^l =2(^„x+j8„) (2.162) 
1 '^n rj 

It is solution Eq.(2.82) and as can be shown Ansatz j8„ = klt + 5n gives rise 

ut + 6uux + Uxxx = (2. 1 63) 

which is important model equation in the theory of surface waves. Formulae 
(2.146) and (2.147) have been generalized for the case of arbitrary by Hirota 
whose work gave a rise to a huge amount of papers dealing with construction of 
soliton-like solutions for nonlinear differential equations because some simple 
trick allows to add new variables in these formulae (see [47] and bibliography 
herein). 
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2.4.3 Finite-gap potentials 

In previous section it was shown how to construct integrable cases of Schrodinger 
equation with soliton-like potentials vanishing at infinity. Obvious - but not at 
all a trivial - next step is to generalize these results for construction of integrable 
cases for Schrodinger equation with periodic potentials. In the pioneering work 
[50] the finite-gap potentials were introduced and described in terms of their spec- 
tral properties but deep discussion of spectral theory lays beyond the scope of this 
book (for exhaustive review see, for instance, [51]). The bottleneck of present 
theory of finite-gap potentials is following: spectral properties formulated by 
Novikov's school provide only almost periodic potentials but do not guarantee 
periodic ones in all the cases. 

We are going to present here some simple introductory results about finite-gap po- 
tentials and discuss a couple of examples. For this purpose, most of the technique 
demonstrated in the previous section can be used though as an auxiliary equation 
we will use not Riccati equation but its equivalent form, modified Schwarzian 
(2.98). 

In particular, it allows us to generalize Lemma 2.1 for the case of arbitrary poly- 
nomial U{x,X) in X. Namely we have 

Lemma 2.11. Equation for modified Schwarzian 

l^-l-^ + m^ = U{x,X):=?i"' + uiX"'-^ + ... + u,n (2.164) 

with any polynomial generalized potential i7(x, A) has unique asymptotic solution 
represented by formal Laurant series such that: 

oo 

h{x,X) = 1 + £ X-%{x) (2.165) 
where coefficients h are differential polynomials in all M]^ , . . . , Uf^i. 

► The proof can be carried out directly along the same lines as for Lemma 2.1. 1 

In particular, for m = 1 which corresponds to Schrodinger equation (2.102) with 
potential U = A -|- m we find 

hi = ^u, 2h2=^hi^xx-hi,... (2.166) 

We remind that for the function 4> = l//i we have (see Theorem 2.4) 
4U{x,X)^^ + ^l-2^^^ = X'^ ^^ = 4U{x,X)^x + 2Ux{x,X)^. (2.167) 
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Definition 2.12. Generalized potential 

U{x,X)^X'" + uiX"'-^ +... + Um (2.168) 

of an equation 

xif,,^U{x,X)xi/ (2.169) 
is called A^-piiase potential if Eq. (2.89), 

(Pxxx = 4U{x,?i)(p, + 2U:c{x,?i)(p, (2.170) 
has a solution ^> which is a polynomial in A of degree A^: 

A' 

4>(x, A) = + (pi(x)A^-i + ... + (Pn{x) = - W)- (2.171) 

7=1 

Roots Pj{x) of the solution <^{x,A.) are called root variables. 

In particular case of Schrodinger equation this potential is also called finite-gap 
potential. As it follows from [51], original "spectral" definition of a finite-gap 
potential is equivalent to our Def. 2.12 which is more convenient due to its ap- 
plicability not only for Schrodinger equation but also for arbitrary equation of the 
second order. 

Example 2.13. Let a solution 4>(x, A) = A — p(x) is a polynomial of first degree 
and potential is also linear, i.e. m = 1, = 1. Then after integrating the equation 
from definition above, we get 

4(?i + u)^'^ + ^l-2^^^ = a(X). (2.172) 

Obviously the "constant of integration" in left part of (2.172) is a polyno- 
mial in A of degree 3 

i^(A) =4A^ + ciA2 + c2A+c3 =4(A-Ai)(A-A2)(A-A3), (2.173) 

where A,- are all roots of the polynomial fl(A). Eq.(2.172) is identity on A and 
therefore without loss of generality we write further C(A) for both sides of it. 
This identity has to keep true for all values of A, in particular, also for A = y{x) 
which gives 

^ = c{r) = 4{r-Xi){r-h){r-h). (2.174) 

Now instead of solving Eq.(2.172), we have to solve Eq.(2.174) which is inte- 
grable in elliptic functions. 

If we are interested in real solutions without singularities, we have to think about 
initial data for Eq.(2.174). For instance, supposing that all Ay are real, without 
loss of generality 

Ai > A2 > A3, (2.175) 
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and for initial data (xq, Yq) satisfying 

V(;co,ro): h<ro<h. (2.176) 

Eq.(2.174) has real smooth periodic solution expressed in elliptic functions 

M = 27-^1-^2 -A3 (2.177) 

It is our finite-gap potential (1 -phase potential) and its period can be computed 
explicitly as 

T= / , (2.178) 

We have regarded in Example 4.8 particular case niN = 1. Notice that in general 
case of niN > 1 following the same reasoning, after integration we get polynomial 
C(A) of degree 2A^+m 

4U{x,X)cp^ + (pi - 2cpcp^ = C{X):= 4X^^+"' + ... (2.179) 

and correspondingly a system of 2N + m — l equations on functions, i.e. the 
system will be over-determined. On the other hand, choice of A = Yj makes it 
possible to get a closed subsystem of N equations for N functions as above: 

yj,x = ciYj)/U(rj-yk?- (2.180) 

Following Lemma shows that this over-determined system of equations has unique 
solution which is defined by Sys.(2.180). 

Dubrovin's Lemma. Let system of differential equations (2.180) on root vari- 
ables Pj defined by (2.171) with 

= 4(^2^+" + CiA2^+'"-1 + • • ■+C2N+m =^ «^ (2.181) 
and ^> = Y[{X — Pk)- Then following keeps true: 

1. a{X)\x^p. = (pl{x,x)\x=p., j=i,...,N, 

2. expression 

4)-i(24)^+ ^^^^~^' ) (2.182) 

is a polynomial in A of degree m and leading coefficient 1. 
< Differentiate the product we get 

-^x = PuYl(^-pk)+ P2,x Yl(^-Pk) + --- (2.183) 
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which imphes by pjk Pj — Pk 

^x\x=p, = -pj,.UPjk ^ ^i^)\x=pj = ^l{xA)\x=pj, (2.184) 

i.e. first statement of the lemma is proven. 

The second differentiation yields 

^xx I X=pj = -P j,xx n pjk + ^Pj,x E n Pi^^Pjk (2- 1 85) 

On the other hand differentiate Ct)(py) we find 



Pj,x-TYh=pj = Pmx llpjk+ (2. 186) 



It remains only to apply I'Hopital's rule to the fraction from (ii). I 
Let a function V{x, A) is defined by Eq.(2.179) as 

4V{x,X) = (p-i(2<p« + ^^^i:^) (2.187) 

(....) Below it will be shown that also transparent potentials themselves can be 

computed algebraically. 

The idea of constructing also transparent potentials by pure algebraic means is 
very appealing because (....) In order to show how to do it, we concentrate further 
not on the coefficients cij of the polynomials and \j/2 but on their zeros. More 
precisely, let us study properties a polynomial 

P=Vi¥2 (2.188) 

and its roots. (....) 



2.5 Summary 

In this Chapter, using Riccati equation as our main example, we tried to demon- 
strate at least some of the ideas and notions introduced in Chapter 1. Regard- 
ing transformation group and singularities of solutions for RE we constructed 
some equivalent forms of Riccati equation. We also compared three different 
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approaches to the solutions of Riccati equation and its equivalent forms. The clas- 
sical form of RE allowed us to construct easily asymptotic solutions represented 
by formal series. Linear equation of the second order turned out to be more conve- 
nient to describe finite-gap potentials for exact solitonic solutions which would be 
a much more complicated task for a RE itself while generalization of soliton-like 
potentials to finite-gap potentials demanded modified Schwarzian equation. 

In our next Chapter we will show that modified Schwarzian equation also plays 
important role in the construction of a differential operator commuting with a 
given one while existence of commuting operators allows us to obtain examples 
of hierarchies for solitonic equations using Lemma 4.6. In particular, for m = 
1 coefficients hk{x) of Eq.(2.165) describe a set of conservation laws for KdV 
equations. 

2.6 Exercises for Chapter 2 

1. Prove that general solution of z — a{x)z has a form 

z(x) = ei'«W^^ (2.189) 

2. Deduce formula (2.57) regarding < y/i , v>2 >= 1 as a linear first order equation 
on ^^2- 

3. Check that Airy function ... is a solution of Eq.(??). 

4. Prove that for L = ^ its kernel consists of all polynomials of degree < m — 1 . 

5. Let functions Ai and A2 are two solutions of (2.89). Prove that the Wronskian 
< Ai , A2 >= A1A2 —A2A[ is solution as well. 

6. Let the function A satisfies (2.90). Prove that the functions 

/± = ^DlogA±^ (2.190) 

satisfies the Riccati equations (2.1). 

7. Proof that 

^4-^—=^' ^ fl=(eie^ + e2e-*^)-2. (2.191) 

8. (Duke Math.J, v.l04(2),2000) , (Nuovo Cimento, v.43,1978) Prove for roots 
variables 

rj^ + 2rj^rj= £ j = l,...,N. (2.192) 

l<i<N,i^j O' '^i 
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9. Proof that Wronskian of two functions (2.1 10) and (2.1 1 1) is well-defined and 
has a form 

(2.193) 



Chapter 3 



Factorization of linear operators 



3.1 Introduction 

Let us notice first tliat different definitions of integrability, as a rule, use lineariza- 
tion of initial equation and/or expansion on some basic functions which are them- 
selves solutions of some linear differential equation. Important fact here is that 
linearization of some differential equation is its simplification but not solving yet. 
For instance, in case of linear Schrodinger equation, Yjcx + = mi/'^, we are not 
able to find its solutions explicitly but only to name them lost functions and to 
exploit their useful properties (see previous Chapters). 

On the other hand, well-known fact is that for LODE with constant coefficients 
operator itself can always be factorized into first-order factors and thus the prob- 
lem is reduced to the solving of a few first-order LODEs: 

which are solvable in quadratures. 

In case of differential operators with variable coefficients factorization is not al- 
ways possible but for the great number of operators BK-factorization gives fac- 
torization conditions explicitly which we are going to demonstrate in the next 
Section. 
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3.2 BK-factorization 

Speaking generally, BK-factorization produces following result: in case of LPDO 
with characteristic polynomial having at least one distinct root, factorization is 
constructed algebraically for an operator of arbitrary order n while in case of some 
multiple roots of characteristic polynomial of LPDO, factorization is formulated 
in terms of Riccati equation(s). Factorization of LODO is always equivalent to 
solving some Riccati equation(s). Below explicit procedure for order 2 and 3 is 
briefly described. 



3.2.1 LPDO of order 2 

Let us outline here BK-factorization procedure [67] for the simplest case of bi- 
variate LPDO of second order. Consider an operator 

A2= ajkd^dy = a2od^ + andxdy + ao2dy+aiodx + aoidy + aoo. (3.1) 

j+k<2 

with smooth coefficients and seek for factorization 

M = {pidx + pidy + P3){P4dx + psdy + pe)- 

Let us write down the equations on /?, explicitly, keeping in mind the rule of left 
composition, i.e. that dx{OLdy) — dx{OL)dy-{- adxy- 

Then in all cases 



«20 = PlPA 

an = P2P4+P1P5 

«02 = P2P5 

aio = J^{p4)+p3P4 + pipe 

aoi = ^{P5)+P3P5+P2P6 

laoo = ^{P6)+P2,P6 



{ISysP) 



where we use the notation ^ = pid^-i- pidy. In generic case we assume that (after 
a linear change of variables if necessary) 



020 7^ and p\ = \. 
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Then first three equations of 2SysP, describing the highest order terms are equa- 
tions in the variables P2,P4, Ps and to find them we have to find roots of a quadratic 
polynomial 

^2{(0) :=a2o(«)^ + aii(«)+ao2 = 



and this leads to a linear system for p4, p^ with co as parameter: 



" 1 0" 




P4 




«20 




P4 




"1 0" 




«20 


-co 1 




.P5. 




an 




.P5. 




(0 1 




an_ 



Thus 




{2Pol) 



and choice of a root (0 generates different possible factorizations of operator A2. 



Having computed P2-,P4-,P5 one can plug them into two next equations of 2SysP 

fllO = -Sf(p4)+P3P4+PlP6 
flOl = ■^iP5)+P3P5+P2P6 

and get a linear system of equations in two variables P3,p6 which can easily be 
solved 

©aio + aoi - (0=^020 --^(«20<o + «ii) 

= ^ TT, ' 

2a2oCO + aii 

{a2oCO + an){aio-^a2o)-a2o{cioi--^{a2oCO + an)) 
2a2o(0 + an 

(3.2) 



P6 



if ^2{(o) = 2a2oQ) + an 7^ 0, i.e. co is a simple root. At this point all coefficients 
pi,P2,---,P6 have been computed and condition of factorization 

aoo = ^(pe) + P3P6 (3.3) 

takes form 



_ Off f O}a\Q+am-.iC{2a20(O+ail) \ , QJaio+aoi-^(2a2oQ)+aii) 
"00 ^ 2a2i)(0+au )^ 2a2oCO+au 

aioiaoi -^(a20(<J+aii))+(a2oQJ+aii)(aio--^a2o) 



(3.4) 



^ 

2a2oG)-\~an 
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3.2.2 LPDOoforder3 

Now we consider an operator 

A3 = ^ Ujkd^dy = azQdl + fl2l d^dy -\-andxdy + fl03<5/ (3.5) 

j+k<-i 

+a2odx + ^ndxdy + ao2^y + awdx + aoidy + aoo- 
with smooth coefficients and look for a factorization 

A3 = {pidx + P2dy + P3){p4dx + P5dxdy + pedy + P7dx + Psdy + pg). 

The conditions of factorization are described by the following system: 

«21 
«12 
«03 
^ «20 
flu 

«02 

aio 

«01 
^«00 

with ^ = pidx + pidy 

Once again we may assume without loss of generality that the coefficient of the 
term of highest order in dx does not vanish, and that the linear factor is normalized: 

«30 7^0, Pl = l. 

The first four equations of SSysVar describing the highest order terms are equa- 
tions in the variables P2-iP4-,P5-iP6- Solving these equations requires the choice 
of a root —p2 of a certain polynomial of third degree. Once this choice has been 
made, the remaining top order coefficients PA^Ps^Pe easily found. The top 
order coefficients can now be plugged into the next four equations of 3SysP. The 
first three of these four equations will now be a linear system of equations in the 
variables P3,pi,p^ which is easily solved. The next equation is now a linear equa- 
tion on variable pg which means that all variables pi,i= 1 , • • • , 9 have been found. 



P1P4- 

P2P4+PIP5 
P2P5+PIP6 

P2P6 

^{P4)+P3P4+P\P7 

'^{P5)+P3P5+P2P7+P1PS 

■^(Pe) +P3P6+P2PS 

■^{Pl)+P3Pl+P\P9 

'^{Ps)+P3PS+P2P9 

■^iP9)+P3P9 



(SSysP) 



3.2. BK-FACTORIZATION 



65 



The last two equations of 3SysP will give us then the conditions of factorization. 



Namely, at the first step from 



f«30 = P4 
«21 =P2P4+ P5 

an = P2P5 + P6 

it follows that 

^3{-P2) ■= a3o{-p2f + a2i{-p2)^ + an{-p2) + ao3 =0. 

As for the case of second order, taking p2 = —(O, where © is a root of the charac- 
teristic polynomial ^3 we get a linear system in paiP5iP6 with (o as parameter. 
Then again 

P2 = -G), ^ = dx-G)dy, 

which leads to 



«30 = Pa 

a\2 = -0)p5+P6 



1 0" 




PA 




«30 




PA 




-CD 1 




P5 






1 


P5 




-a 1 




P6 




an 




P6 





(0 





1 

(0 1 





^30 




«21 




_ai2_ 



Thus 



Pi 


= 1 


P2 


= -a 


PA 


= «30 


P5 


= (23o£0 + a21 


^P6 


= a'iQOp- + a2\a) + an 



{3Pol) 



At the second step, from 
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j!f{p4)+P3P4 + PlP7 
■^(Ps) +P3P5+P2P7+P1PS 
■^{P6)+P3P6+P2PS 



and (2Pol) we get 



020 -^«30 = P3«30 + P7 

ail -jSf(a30<o + a2i) = p3(a3oCO + a2i) - copj + ps (3Lm*) 

ao2 -^(«30<»^ + «21<O + «12) = /73(«30<»^ + «21<O + «12) - <0P8- 
As a linear system for P3, p-j, ps this has determinant 

3a3o(0^ + 2fl2i<o + ai2 = =^3(0), 
so if (0 is a simple root the system has unique solution 

^^^(^^20 ~-^^3o) + tt>(aii — ^ (a^oO) + a2i)) +ao2 — ^(asott)^ + a2itt) + 012) 

P3 ^ 5 

3fl3ott)^ + 2a2i W + ai2 

020 — ^«30 «30 

Pi ~ * P3 5 

3a3oCO^ + 2a2i(0 + ai2 3a3o<o^ + 2a2i<o + «i2 

(O(fl20--^«30)+«ll -■^(«30<O + «2l) fl30<O + a2l 
PS ~ * P3 • 

3a3oCt)2 + 2a2i © + 012 3fl3o<o2 + 2a2i co + ai2 

In order to find the last coefficient pg we use the next equation of (SSysVar), 
namely: 

aio = ^{pj)+P3P'j+pipg, (3Lm**). 

At this point all coefficients pi,i— 1, ...9 have been computed, under the assump- 
tion that (0 is a simple root. 

At the third step from 

flOl = (P8 ) + P3PS + P2P9 ^2 6) 

fl00= ■^{P9)+P3P9 



all the necessary conditions for factorization can be written out. We do not do so 
here because the formulas are tedious and do not add anything to understanding 
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the main idea. If the conditions are satisfied, the expUcit factorization formulae 
could be written out as for the second-order operator. The difference is that in this 
case the polynomial defined by the highest order terms is of degree 3 and we have 
not one but two conditions of factorization. 



3.2.3 Constant coefficients 

• Obviously, in case of constant coefficients atj all formulae above can be 
simplified considerably and used for classical factorization problem of a 
polynomial. For instance, a bivariate second order polynomial 

W = X^-Y^ + aioX + aoiF + aoo 

can be factorized into two linear polynomials, 

X^-Y^ + aioX + aoii" + «oo =ipiX + piY + ^3) {PA^ + P5Y + P6): 

iff 



'^'01 ^'10 



4 

In each case coefficients pi can be written out explicitly, for instance if 



2 2 

«oo = 



^01 ^10 



then 

X -Y + aioA + aoir + aoo = (a + F H 7, )(^-FH ). 



As in case of order two, constant coefficients aij simplify all the formulae 
and reduce the problem under consideration to the classical factorization of 
a polynomial. For instance, a bivariate third order polynomial 

W = X^Y+ XY^ + a2QX^ + ai iXF + 002^^ + aioX + ooiF + ooo 

can be factorized into the product of one linear and one second order poly- 
nomials, 

X^Y +XY^ + a2o^^ + anXY + uqiY^ + uxqX + aoiF + aoo = 
= {piX + P2Y + p3) {P4X^ + PsXY + peY^ + pjX + p^Y + pg) 
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for instance, if 



«00 = («11 — «20 — «02)[«10 + «20(«11 — «20 — «02)]5 



and the result of factorization then is (with notation y = an — a2o — ao2): 
X^Y +XY^ + a2oX^ + anXY + ao2Y^ + uiqX + aoiY + am = 
= {X + Y + y)(XY- a2oX + (fl20 - an + r)Y + «io + «2o7) • 

3.3 Laplace transformation 

3.3.1 Main notions 

The most important question now is - what to do when conditions of factorization 

are violated? Do we still have a way to solve an equation ^{^1^) = correspond- 
ing to the initial operator? In order to answer these questions let us re-write results 
of BK-factorization for generic case of second order hyperbolic operator as 



and corresponding LPDE as {dxdy + adx + bdy + c)^f\=Q and introduce new func- 
tion = {dy + a) v^i . Our main goal now is to construct some new LPDE having 
1/^2 as a solution and to check its factorization property. If this new LPDE is fac- 
torizable, then its solution is written out explicitly and due to the invertibility of 
a transformation y/i \j/2 the formula for solution of initial LPDE can also be 
obtained immediately. 

Let us first introduce some definitions. 
Definition 3.1 Two operators of order n 



are called equivalent operators if there exists some function / = f{x,y) such that 




(3.7) 



=^1= E ajkdjdy' and =^2 = £ Mi^y 




/jSfl=^2o/. 
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Definition 3.2 Expressions 

d = ab + Ux — c and b — ab + by — c 
are called Laplace invariants. 

Lemma 3.3 Two hyperbolic operators and J^2 of the form (3.7) are equiva- 
lent iff their Laplace invariants coincide pairwise. 

► Indeed, 

f^i{Y)=fdxdy + faidx + fbidy + fci 

and 

-^lifV) = dxdy{f)V+dyfdx¥+dxdy{\j/)f + dy\l/dxf+ 
aidx if) V + a2dx{\lf)f + &2^y (/) + ^2<3y ( V)f + f¥C2 , 
i.e./J5fi = J5f2o/iff 

fli = fl2 + ^ = fl2 + (log/)y; 
&l = &2 + y^ = ^2 + (log/).; 

CI = C2 + ^ + a2^+b2^ =C2 + (l0g/)xy + (log/);,(log/),+ 

+«2(l0g/);c + &2(l0g/)y. 

Direct substitution of these expressions into formulae for Laplace invariants gives 
(we use notation 9 = log/): 

di=aibi+ai^x-ci = ia2 + qfy)ib2 + q>x) + (02 + 9y)x- 

C2-(Pxy- (Px(Py - ai^Px - ^(Py = «2^2 + «2,;c - C2 = a2- 

Analogously one can obtain bi = S2 and it means that for two equivalent hyper- 
bolic operators their Laplace invariants do coincide. I 



► First of all, let us notice that two operators 

^i = (dx + bi)(dy + ai)+Ai 
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and 

^2={dx + b2){dy + a2)+A2 
can be transformed into some equivalent form 

and 

^2 = {d, + b2)dy+A2 
(perhaps) by different functions /i , /2 such that 

and as was proven above, Laplace invariants of the initial operators coincide with 
those of the equivalent ones. Operator of the form 

has following Laplace invariants: A and —by +A, i.e. Ai = A2 and biy = b2y for 
operators JSfi and ^2 with the same Laplace invariants. It yields to 

bx=b2 + (p{x), ^2^{dx + bx + (p{x))dy+Ax 

with some arbitrary smooth function <p(x). Now operator ^2 differs from 
only by term ^{x) which can be "killed" by one more equivalent transformation, 
namely, for some function /3 = /3 {x) 

A = f^^^2 o h = {d, +b, + (p{x) + (log/3 + Ai 

and choice (p{x) f completes the proof. I 

Now let us rewrite initial operator (3.7) as 

^1 : dxdy + aidx ^-b\dy^- ci 

and notice that 

{dx + bi)\i/2 = di\i/i 

leads to 

r (^x+^i)v^2 = aiv^i ^ f i^{dx+bi)\ir2 = m 

\ (<9y + ai)v/i = v/2 1 {dy + ai)y/i^\}/2 

and standard formula for log-derivative: 

e^dye~^ = dy — (Py with 9 = logo, 
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gives finally a new operator ^2 with corresponding LPDE 



V^2 = 0. (3.8) 



In order to check whether these two operators ^1 and ^2 are different, let us 
compute Laplace invariants of the new operator 

a\ =aibi+ai^x-ci j a2 = fli (logai)xy + ai,;c 

b\=a\b\+b\^y-c\ \b2 = a\ 

Now one can see that operators J^i and J^2 are not equivalent and operator ^2 is 
factorizable if 02 = or &2 = (see example below). 



Definition 3.4 Transformation ^1 =^ J^2, i-e. 

bi =aibi+bij-ci l^2 = «i 

is called Laplace transformation. 

If first new operator is also not factorizable, the procedure can be carried out for as 
many steps as necessary in order to get some factorizable operator. At the step 
when the first factorizable operator is found, algorithm stops because the division 
on corresponding = is not possible any more. In fact, it is possible to write 
out formulae for Laplace transformation in terms of Laplace invariants only. 



Tlieorem 3.5 Let Un is one of Laplace invariants an,bn obtained at the step n. 
Then 



Un+l=2Un+{\ogUn)xy-Un-l. (3.9) 

► Indeed, due to Lemma 3.3 it is enough to regards sequence of operators of the 
form 

■^n '■ dxdy + andx + bndy + Cn with a„ = 
because it is easy to find some function fa„ (for instance, /a„ = /«n^y) such that 

K,Q=f-„^.^nOfa„ = dxdy + bndy + Cn. (3.10) 
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From now on 

jSf„ : dxdy + bndy + Cn 

and tilde-s are omitted for simplicity of notations. Now formulae for Laplace 
transformation take form 

¥n,y = -Cn¥n+h Yn = [dx + bn + {^OgCn)x\Vn+\ (3.11) 

and Eq.(3.8) can be rewritten for the function Yn+i as 

[dx + bn+ {logC„)x]dy\l/n+l + [Cn + bn,y + (logC„)xy] Wn+1 = 0, 

i.e. 

Cn+l -C„ = bn,y + (logC„)xy, &„+l - = {\ogCn)x, 
Cn+1 - Cn+l = bn+l,y + (logC„+i)xy, &„+2 - &n+l = (logC„+i)x 

and finally 

Cn+2 = 2c„+i + (logc„+i)xy-c„. (3.12) 
Notice that in this case, first Laplace invariant is 

and obviously satisfies to Eq.(3.12), i.e. for the first invariant the statement of 
the theorem is proven. In order to prove it for the second invariant bn one has to 
choose another sequence of operators with = generated by some function fi,^ 
such that 

^n,0 = fb!^nOfb, = dxdy + andx + Cn. (3.13) 

I 

Notice that in order to obtain the recurrent formula for Laplace invariants, we used 
separation of variables x and y given by (3.11). Moreover, introduction of a new 
discrete variable n allows us to regards these equations as difference-differential 
ones. In order to deal with this sort of equations one needs a couple of definitions. 

Definition 3.6 An operator T acting on the infinite sequences of functions 

(..., v^_2, V^-i, V^o, Wh V2, Wn, •••) 

as 
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is called shift operator. For convenience of notation sometimes infinite vector- 
function is introduced 

V/oo = (..., Y-2, ¥-1, ¥o, ¥h ¥2, ¥n, •••) 
and matrix of operator T then has the following form: 



T = 



\ 






1 



















1 



















1 





















1 




















(3.14) 



i.e. it is infinite matrix with all zero elements but the elements over main diagonal 
- they are equal to 1. 



Definition 3.7 Commutator ^ =[A,B] of two operators A and B is defined as 

'^^AB-BA. 



Obviously, following properties hold true: 

• [dx, dy] = (cross-derivative rule), 

. [dx,T] = [dy,T]=0, 

• [T ,a] — T o a — aoT — [T (a) — a) oT (Leibnitz rule analog). 

Let us now regard two operators corresponding Laplace transformations from 
Theorem 3.5 rewriting slightly formulae (3.11) in terms of shift operator: 

¥n,y^-CnT¥n, ¥n,x + bn¥n ^ T^Wn- (3.15) 

The use of shift operator makes it possible to present their commutator 

'^=[dy + CnT, dx + b„-T-^] = [dy + cT, dx + b-T-^] 
omitting low index n. 
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Lemma 3.8 Commutator 

'rf^ldy + cT, d^ + b-T-^] 

is equal to zero iff 

Cx = c{T{b)-b), by = c-T-\c). (3.16) 
► Indeed, by definition 

<i^ = ljy-CxT + cTb-bcT + T-^cT-c=-c+by+T-^{c)-{cx-cT{b)+bc)T. 

Now, if ^ = 0, then -c + by + T-\c) = and C;, - (Z?) + Zjc = 0, i.e. (3.16) 
holds. 

If (3.16) holds, then coefficients of are equal to zero, i.e. ^ — O.t 

At the end of this Section let us notice that in the original Eq.(3.7) two variables x 
and J played symmetrical role which can also be observed in commutation relation 
of Lemma 3.8 after appropriate gauge transformation: 

e^" {dy + CnT)e-''" = dy-an + T, e"" {dj, + bn + 1-^)6-"" = + Cn-lT'^ 

with 

~ ^"^ bfi = (Jrixj ~ 9ny 

It can be shown that qn satisfy the following equation 

?n,xy = e^"+^"^"-e^''"^"-'. (3.17) 

The equation (3.17) is usually called two-dimensional Toda lattice and plays 
fundamental role in the theory of Laplace transformations. 

3.3.2 Truncation condition 

Definition 3.9 Truncation condition for Eq.(3.9), namely 

(logM„)xy = M„+i -2m„ + m„_i, n = l,...,N+ 1 
is defined by Dirichlet boundary conditions, i.e. uq = Uf^+\ — 0. 
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For example, case N =1 gives us Liouville equation^ 

(10gMl)xy + 2Ml =0 



while case N = 2 yields to the system 



(logMl)xy = M2-2M1 + M0 ^ \ {logUl)xy = -2ui + U2 
(logM2)xy = M3-2M2 + Ml \(loga2)xy= Ml - 2m2 

These both cases are known to be integrable in quadratures. 

Case of arbitrary given N corresponds to the system of equations with following 
matrix 



An 



/ -2 1 
12 10 
12 1 



\ 














1 





-2 1 
1 -2 



(3.18) 



on the right hand. This matrix is some Cartan matrix and another choice of 
boundary conditions leads to Cartan matrices of other form. It is interesting to 
notice that all Cartan matrices can be constructed in this way. Moreover, it is 
proven that system of equations corresponding to each Cartan matrix is integrable 
in quadratures ([69]). Most well-known source of Cartan matrices is semi-simple 
classification of algebras Lie - there is exists one-to-one correspondence between 
these matrices and semi- simple algebras Lie ([69]). 

Explicit expression for any Laplace invariant m„ after n Laplace transformations is 
given by following lemma ([69]) which for simplicity is formulated for the special 
case Toda chain. 



Lemma 3.10 Let us regard an infinite sequence of functions 

{dn}, n = 0,1,2,... 

such that 



dxdylogdn = 



J2 



Vn 



'See Ex.1 
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dn = dQt{d^d^w) , ij^l,...,n. 



Corollary 3.11 Sequence of functions 




is solution of Eq.(3.9) while sequence of functions 



dn-l 



is solution of Eq. (3.17). 

Theorem 3.5 describes an infinite chain of equations corresponding to Laplace 
transformations and to start with this chain, we need nothing more then two in- 
variants. On the other hand, many applications of this theorem are connected with 
some special problems in which different sort of finite chains are considered. In 
the next Section we will discuss two most usable ways to construct some finite 
chain of invariants and we close this Section with an example of 2-steps chain. 

Example 3.12 Let us regard operator 
then its Laplace transformation gives 




and operator ^2 is factorizable: 



It is a simple task to write out explicitly solution Y2 of LPDE 



(3.19) 
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and afterwards solution i/^i of 

^i(V^i) := {d,dy+xd, + 2){Wi)^Q 
can be computed by formula 

Wi = -dx¥2- (3.20) 

Indeed, introducing in Eq.(3. 19) notation (p = ((9y + ^) we find that q) = Y{y) is 
arbitrary function of one variable (...) 



V2 



3.3.3 Periodic closure 

Definition 3.13 Classical periodic closure for equation on Laplace invariants 

{logUn)xy = Un+l-2Un + Un-l, n= l,...,N 

is defined by periodic boundary conditions, i.e. Un+N = 

In this case Cartan matrix is replaced by matrix An and for > 3 its form is 



An 



(-11 
1-210 
12 1 





VI 



1 





... 1 

... 

... 

-2 1 

1 -2 



(3.21) 



Notice^ that matrix Aa^ is degenerated. It will be shown below that for N = I 
initial equation can be solved explicitly, while for N — 2,3 initial system of equa- 
tions allows some reduction to one scalar equation. 



Let us regard first case |A^= 1 1 , it yields to (logMi)xy = M2 — 2mi + mq = and 
obviously ui = gi{x)g2{y) with arbitrary smooth functions gi{x),g2{y)- 



2Ex.2 
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Case I = 2 1 is more interesting due to the huge amount of appUcations (surfaces 
with constant curvature, relativity theory, etc. ) and gives rise to the system of 
equations 

(log"iW = 2(Mi-M2), ^2 22) 

(logM2)xy = 2(M2-Ml) 

with many important properties: it has conservation laws, symmetries, soliton- 
type particular solutions, etc. In particular, the reduction of this system can easily 
be constructed to the one scalar equation: 

(logMi)xy + (logM2)xy = =^ (logMiM2)xy = =^ U1U2 = X {x)Y {y) 

with two smooth arbitrary functions X{x),Y{y). Suppose now u\U2 = 1, then set- 
ting Ml = , U2 = we find solutions of Sys.(3.22) from solutions of equation 

0xy + sinh0 = O. (3.23) 
This equation is called sinh-Gordon equation [70]. 



Case N =3 corresponds to the system 



(l0gMl);,y 


= U2 + UT, 


— 2mi 


(l0gM2)xy 


= Ml +M3 


-2m2 


(10gM3)xy 


= Ml -|-M2 


-2m3 



which allows reduction mi = M3 

(log Ml )^ = M2 + M1 -2mi 
(log M2)xy = Ml + Ml - 2m2 
{l0gUl)xy = M1-I-M2-2M1 



(logMi);,^ = M2-M1 
(logM2)xy = 2(Mi-M2) 



and it can be treated analogously with the case above: 

M1M2M3 = 1, U2 = e , ui = e , M3 = e =^ 

0^ = e-^^-e^. (3.24) 

Eq.(3.24) is called Tzitzeica equation and its solutions give solution of initial 
system of equations. Tzitzeica equation is also very important for various appli- 
cations [71]. 



3.3. LAPLACE TRANSFORMATION 



79 



For both Eq.(3.23) and Eq.(3.24) their general solutions are not available as well 

as for the case of general N. Attempts to solve appearing systems of equations 
directly demand some tedious technique of inverse scattering and produce partial 
solutions with singularities. Method to reduce the initial system to one equation 
simplifies drastically construction of smooth solutions. 

At the end of this section let us notice that the same three equations which were 
obtained while studying truncated and periodical cases - namely, Liouville, sinh- 
Gordon and Tzizeica equations - do appear together in some other context [72]: 

Theorem. Nonlinear PDE of the form Uxy = f{u) has higher symmetries iff one 
of three cases take place: 

/ = e", / = e" + e-", / = e" + e-2«. (3.25) 



As we know already (see Chapter 1) that integrability of a differential equation is 
intrinsically related to its symmetry properties. Of course, possession of a symme- 
try does not mean that equation is integrable in some sense but this fact gives us 
a good hint on what equations might be integrable. Moreover, for some classes 
of differential equations it is proven [73] that integrability is equivalent to some 
well-defined symmetry properties. 

From this point of view the theorem above justifies hypothesis that two of these 
three equations (integrability in quadratures of Liouville equation we have demon- 
strated already) have good integrability properties. As we will show in the next 
section, integrability of these PDEs can be reduced to the integrability of well- 
know ODEs. 

Thus, in contrast to the truncation condition which leads to integrability in quadra- 
tures for arbitrary order A'^, periodical closure generates a more complicated situa- 
tion. Here we are not able to get the answer directly in terms of invariants but we 
need first to study the properties of solutions at the all N steps simultaneously. 

3.3.4 Separation of variables 

Obviously, in case of periodical closure with period A'^ functions ^f\ and ^^N+\ 
satisfy the same equation which was presented in Theorem 3.6: 
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Wn^ + bnYn,y + CnYn = 0, V n = 0, ±1, ±2, ... (3.26) 

or, in matrix form, 

idxdy + {b}dy+{c})% = 

with diagonal matrices {b} and {c}. This fact allows us to regard finite vector- 
function Y, 

which defines completely all the properties of the initial infinite system of equa- 
tions. Notice that the fact of coincidence for coefficients of two LPDEs does not 
mean that their solutions Yn+N and !//'„ also do coincide - they might differ, for 
instance, by a constant multiplier. Therefore we need now some notion of periodic 
closure for solutions. 



Definition 3.14 Blocli periodic closure is defined for the components of func- 
tion \ff as follows: 

Notice once more that truncation and classical periodical closure are defined for 
Laplace invariants, i.e. for coefficients of Eqs.(3.26), while Bloch closure deals 
with solutions of the same equations. 



Corollary 3.15 Shift matrix in case of Bloch periodical closure has form 



Tn = 



\ 

where e C is a free parameter. 






1 



















1 



















1 





















1 



















(3.27) 



Corollary 3.16 Basic Lemma 3.8 holds true also for the case of periodic clo- 
sure, i.e. commutator of Laplace transformations 

iff 

Cn,x = Cn{bn+1 - bn) , &n,y = C„ " C„_l , n G In, 

which can be checked directly. 
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Example 3.17 As it was shown above, in case of N = 2 which corresponds 
to the classical periodic closure of chain of invariants degenerates into sinh- 
Gordon equation and this closure has the form (3.22). Corollary 3.14 allows us 
to construct a connection between closure for invariants and closure for solutions 
of(3.11) : 



















-c{ 






T2 


X 


1 -b2 




y2_ 




T2 


y 


_-k^C2 




T2. 



Then 

Wix + bi\i/i= k-^Y2, Wi =b2W2 + ¥2x, 

and excluding i/Ai or 1/^2 we get an equation of the second order on one scalar 
function, for instance 

¥2,xx + ih + b2)\\f2x + {b2x + bib2- k-^)\\f2 = 0. 

This equation is obviously equivalent to linear Shrodinger equation 

Wxx={^ + u)\\f, X=k-^, 

i.e. Sinh-Gordon equation is S-integrable. This important fact plays role, for in- 
stance, while constructing surfaces of constant curvatures (see very exhaustive 
review [74]). 



3.4 General invariants and semi-invariants 

Before discussing the notion of general invariant, let us notice that arbitrary 
LPDO of second order, A2, can be represented in the form of factorization with 
reminder 

where reminder I2 is defined by (3.3): 



r ©aio + aoi - J^(2a20<o + «ii) 1 
y 2a2o(0 + an ) 



©aio + aoi - J^(2a2o<o + «ii) 
2a2o(0 + an 



a2Q{aoi-^{a2oCO + an)) + {a2oCO + an){aio--^a2o) 
2a2oCO + au 
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Similar to the case of order two, arbitrary LPDO of third order A3 can be repre- 
sented in the following form 

^3 = {Pidx+P2dy+p3){p4dx+P5dxdy+P(,dy+pTdx + psdy+p9)-hdy-hl- 

In contrast to the second order LPDO, in this case factorization with reminder 
gives us not a function but a linear first order operator and it is convenient for our 
further investigations to regard in this case two "reminders" and /31 which are 
defined by Sys.(3.6): 

^3 = «0l - {Pidx + Pldy + P3)Pi - P2P9 
hi =cioo-{pidx + P2dy+P3)P9- 

In the next Section it will be shown that "reminders" I2, I3, I31 are invariants of 
the corresponding operators uner the equivalence transformations. 

3.4.1 Construction of invariants 

Let us first recollect definition of two equivalent operators. 

Definition 4.1 Two operators of order n 

=^1= E ajkdidy' and £ Mx^' 

j+k<n j+k<n 

are called equivalent operators if there exists some function / = f{x,y) such that 

The definition is given for an operator jSf of arbitrary order n and obviously any 
factorization 

can be written out in equivalent form 

r 1^ o / = /-1^1^2 o / = (/-l=2f2 o /), 

as well as sum of operators ^ = ^2'- 

f-'^of = + jSf2) 0/ = (/-l^i of) + (/- 1^20/). 
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Below we will take function / in a form f = for convenience. In order to 
formulate theorem on invariants we need following notations: 

= a2od^ + andxdy + ao2dy+aiodx + aoidy + aoo, 
Mp = {pidx + P2dy + P3){p4dx + P5dy + pe) -h, 

Ma = a3odx + a2id^ dy + andxdy +<3o3<?j +«20<?x + 
+andxdy + ao2<?y + aiodx + aoidy + aoo, 

= {pidx+P2dy + P3){P4dx+P5dxdy+P6dy + 

+Pldx + p^dy + P9)-l3dy-hh 
Pi = r^Pi ° f, at J = f'^aij o f, 

Ai = f^^AiO i — 2a,2p,3a,3p. 

Above A2p = A2a, i-C- A2p and A2a are different forms of the same operator - its 
initial form and its form after the factorization with reminder. The same keeps 
true for and As^, i.e. Asp = As^. 



Theorem 4.2 For an operator of order 2, its reminder I2 is its invariant under 
the equivalence transformation, i.e. 

I2 = h. 

For an operator of order 3, its reminder I3 is its invariant, i.e. 

while reminder /31 changes its form as follows: 

hi = hi + hfy- 



► Indeed, for operator of order 2 

Ma = Mp-h, 
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I.e. 



Ma = r'A2a O / = r\A2p - h) O f 
= A2p-r\h)of = A2p-h. 



For operator of order 3 



A3a=A'ip-hdy 



I.e. 



A3a = r^A^a o / = r\A^p -l^dy-hl)of = 

hp - (hi + hdy) o / = A3p - hdy -h<py- h\ .1 



Corollary 4.3: If /s = 0, then /31 becomes invariant. 
That is the reason why we call /31 further semi-invariant. 

Corollary 4.4: If ^ 0, it is always possible to choose some function / : 

/3I = /3</>y + /31 =0. 

Notice that for second order operator, if its invariant /2 = then operator is fac- 
torizable while for third order operator two its invariants have to be equal to zero, 
^3 = ^31 = 0- On the other hand, if operator of third order is not factorizable we 
can always regard it as an operator with only one non-zero invariant. Of course, 
all this is true for each distinct root of characteristic polynomial, so that one ex- 
pression, say, for l^, will generate three invariants in case of three distinct roots 
of corresponding polynomial. Expressions for invariants I2 and and also for 
semi-invariant /31 can be easily written out explicitly using formulae given by 
BK-factorization (Section 2.1 and 2.2). 

As it was show already, for an important particular case - hyperbolic operator of 
second order in the form 



- there exist two Laplace invariants which coincide pairwise for equivalent opera- 
tors (Lemma 3.3). After rewriting hyperbolic operator in the form 



by appropriate change of variables, we can construct Laplace invariants as a sim- 
ple particular case from the formulae for general invariants (see next Section). 



dxdy + adx -\-bdy-\-c 



(3.29) 



d^-dy+ adx + bdy + c 



(3.30) 
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Corollary 4.5: Two hyperbolic second order operators having the same normal 
form, (3.29) or (3.30), are equivalent iff their general invariants coincide. 

3.4.2 Hierarchy of invariants 

As it was shown above, every general invariant is a function of a distinct root O) 
of the characteristic polynomial and each distinct root provides one invariant. It 
means that for operator of order n we can get no more than n different invariants. 
Recollecting that BK-factorization in this case gives us one first order operator 
and one operator of order n — 1, let us put now following question: are general 
invariants of operator of order n—l also invariants of corresponding operator of 
order n? 

Let regard, for instance, operator of order 3: 

^3a = =AiA2a -hdy-l3i=Ai (Ajp - h) -hdy-h\= A\A2p -hA\- hdy - /31 
and obviously 

= AiA2p - hA\ -hdy- /31 = A\A2p - hM - hdy - h\ , 

i.e. I2 is also invariant of operator A-^a- Let us notice that general invariant /s = 
l-i{(0^^^) is a function of a distinct root (O^^^ of the polynomial 

=^3 {z) = a^.QZ' + a2\7^ + a\2Z + flo3 

while general invariant h = hico^^^) is a function of a distinct root ©^^^ of the 
polynomial 

^2{z) = P4Z^ + P5Z + P6 

with p4, p5, p(, given by {3Pol) for (O = co^^\ In case of all distinct roots of 
both polynomials ^3(2) and M2{z), one will get maximal number of invariants, 
namely 6 general invariants. Repeating the procedure for an operator of order n, 
we get maximally nl general invariants. In this way for operator of arbitrary order 
n we can construct the hierarchy of its general invariants 

■) In— li •••lil 

and their explicit form is given by BK-factorization. 

For instance, let us regard a third order hyperbolic operator in the form 
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C = a-iod^ + a2\ d^dy + a^dxdy + ao3<9y^ + terms of lower order (3.31) 

with constant high order coefficients, i.e. aij = const \/ i + j = 3 and all roots of 
characteristic polynomial 

are distinct and real. Then we can construct three simple independent general 
invariants in following way. Notice first that in this case high terms of (3.31) can 
be written in the form 

{aidx + fi\dy){a2dx + fi2dy) {a-idx + fi^dy) 

for all non-proportional aj,/3, and after appropriate change of variables this ex- 
pression can easily be reduced to 

dxdy{dx-\-dy). 

Let us introduce notations 

d\ = d2 = dy, d3 = di + d2 = <9f , 
then all terms of the third and second order can be written out as 

Qjk = {di + ai){dj + aj) (4 + = didjdk + a^didj + a jdidk + aidjdk+ 
+{dj + Uj) {ukjdi + {di + Qi) {ak)dj + {di + Qi) {aj)dk + {di + ai){dj + a j) {aj,) 

with 

a2o = a2, ao2 = ci\, an = a\+a2 + aj, 

and Cijk = C — Cij^ is an operator of the first order which can be written out ex- 
plicitly. As it was shown above, coefficients of Cij^ in front of first derivatives 
are invariants and therefore, any linear combination of invariants is an invariant 
itself. These invariants have the form: 

for C123 we have 

fl2«3 + «l«2 + dy{a3) + dx{a2) - aio, axa^, + a\a2 + <?x(«3) + <?x:(«2) - «oi; 
for C312 we have 

fl2«3 + «l«2 + <5f(«2) + <?«:(«2) — «105 «l«3 + a\a2-\- dt{ai) -\- dx{a2) — aoi; 
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for C231 we have 

a2a2 + aia2 + dy{a^) + dy{ai) — aio, «i«3 -\-a1a2 + («i) + dy{a\) — aoi- 

Direct calculation gives us three simplest general invariants of the initial operator 
C: 

Proposition 4.6 General invariants /21, /32, /31 are all equal to zero iff operator 
C is equivalent to an operator 

L = did2d3 + bidi+b2d2 + c, (3.32) 

i.e. 3 function / : f~^Cof = L. 

► Obviously 

r\d,dydt) Of={d,+ (log/),) {dy + (log/),) {dt + (log/),) 

for any smooth function /. Notice that it is the form of an operator and 
introduce a function / such that 

«1 = (l0g/)x, «2 = (l0g/)y, 03 = (log/);. 

This system of equations on / is over-determined and it has solution /o iff 

«2,x - = 0, a^j - fl2,f = 0, fl3,;c = 0, 

i.e. /21 = ^32 = ^31 = 0. 1 

■4 Indeed, if C is equivalent to (3.32), then a2o — «02 = an = and obviously 
hi = hi = hi = 0. 1 

At the end of this section let us notice that the reasoning above can be carried 
out for the hyperbolic operator of order n with constant leading coefficients and 
analog of Proposition 4.6 keeps true with n order terms as 

did2d3....dn 

with 

di = dx, d2 = dy, d3 = dx + dy, 84 = a^dx + ^Ady, ... = Ondx + ^ndy 
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and terms of order n and n—l can be written as 

Cii,i2,-,in = {di+ai){d2 + a2)...{dn + an) 
while corresponding linear general invariants will take form 

lij = di{aj)-dj{ai). 

3.4.3 Examples 

Operator of order 2 

Let us regard first for simplicity LPDO of the second order with constant leading 
coefficients, i.e. 

aij = const V(? + j) = 2 

and all roots of characteristic polynomial are distinct. Then obviously any root CO 
also does not depend on x^y and expressions for pij can be simplified substantially. 
Let us introduce notations 

aoo<o°= =^o(<o), 
aioo^+aoi = ^i{(0), 
a2oG)^ + aiiCO + aQ2 = ^2{(X>) 
and notice that now CO and ^2(0)) 7^ are constants. 

Using formulae from Section 2. 1 we get 



which yields to 



P6 = 



Pl = l 

P2 = -CO 

_ ^ijco) 

PA = «20 
P5 =«20<O + «ll 
(a20<O + «ll)«10 — «20«01 



a2QC0 + ai\ a2o ^, . 
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and invariant h takes form 

h = -^{pe) -P3P6 + aoo = ^^^t^^(flio) - ^^=Sf(aoi)+ 

+^ ^) 



Let us regard hyperbolic operator in the form 

+ aiodx + aoidy + aoo, (3.34) 

i.e. a2o — = 0,ao2 = —1 and CO — ±1, £^ — d^ — (ody. Then /2 takes 
form 

7 r>9/®^io~^OK ©aio — aoi ct)<3io + «oi 

^^ = "«°-^^^^) 2^^ 2^ 

which yields, for instance for the root 0) = 1, to 

/2 = «00--^ I j ^ClW)-[.Ox-Oy)K ) 

and after obvious change of variables in (3.34) we get finally first Laplace 
invariant d 

I2 — C — djcU — ab = d, 

where 

«io— «oi , «10 + «01 
a= , b= , c = floo- 

Choice of the second root, ft) = — 1, gives us the second Laplace invariant 
b, i.e. Laplace invariants are particular cases of the general invariant so that 
each Laplace invariant corresponds to a special choice of ft). 

Let us proceed analogously with an elliptic operator 

dxx + dyy + aiQdx + aoi dy + aoo, (3.35) 
then ft) = ±i, = dx — (ody and 

l2 = «00 + [Ox T lOy){ 2 4 

where choice of upper signs corresponds to the choice of the root ft) = / and 
choice of lower signs corresponds to ft) = — /. 
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Operator of order 3 

Now let us regard LPDO of the third order with constant leading coefficients, i.e. 

Qij = const V(/ + 7) = 3 
with at least one root distinct of characteristic polynomial 

030©^ + a2i (O^ + ai2<0 + ao3 = ^3(0)) 
and notice that now (O and (ft)) 7^ are constants. 



Using formulae from Section 2.2 we get 



Pl = l 

P2 = -03 



Pl = 



PS = <33oa» + a2i 

P(y = ^30 + (321 w + ai2 
«20 «30 <^2(«) 



P8 = 



0)a2o + ciii coa^o + a2i ^2{(o) 

""^^ ^^(co) =^^((0) =^^((0) ^'^{(oy ^'^{(0) ^'^{(o) ^'-^{(oy 

and /3,/3i are 

-{p\dx + P2dy + P^)pi-piP9 + ao\ = I3, 

-iPidx + P2dy + P3)p9 + aoo = hi- 

The formulae are still complicated and in order to show the use of them let us 
regard here one simple example of an operator 



B = dxdy + dxdy +an dxdy + aiodx + aoi dy + aoo , 
with 030 = ao3 = 020 = «02 = 0, 021 = an = 1. Then its invariant 

h = dxan -«oi 



(3.36) 



3.5. SUMMARY 
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and semi-invariant 

/31 = dxaiQ — aoo 

have very simple forms and gives us immediately a lot of information about the 
properties of operators of the form (3.36), for instance, these operators are factor- 
izable, i.e. has zero invariants I3 — I31 —0, iff 

an = J aoxdx + fiiy), aio = j aoodx + f2{y) 

with two arbitrary functions on y, fi{y) and /2(y). Another interesting fact is 
that if coefficient an = an{y) is function of one variable y, then aoo is general 
invariant and there definitely should be some nice geometrical interpretation here, 
etc. 



3.5 Summary 

At the end of this Chapter we would like to notice following very interesting fact 
- beginning with operator of order 4, maximal number of general invariants is 
bigger then number of coefficients of a given operator, 

'" + ')'" + ^)<n.Vn>3. 
2 

It means that general invariants are dependent on each other and it will be a chal- 
lenging task to extract the subset of independent general invariants, i.e. basis in 
the finite space of general invariants. 

As to semi-invariants, notice that an operator of arbitrary order n can always be 
rewritten in the form of factorization with reminder of the form 

lndl^ + ln,idl^-^ + ... + ln,k-i, k<n 

and exact expressions for all /; are provided by BK-factorization procedure. The 
same reasoning as above will show immediately that /„ is always general invari- 
ant, and each /„ yt-jo is iQ-th semi-invariant, i.e. it becomes invariant in case if 
ln,k-i = 0, Vi < I'o- 

In this paper, explicit formulae for I2, h, hi are given. Formulae for higher or- 
der operators can be obtained by pure algebraic procedure described in [67] but 
they are too tedious to be derived by hand, i.e. programm package for symbolical 
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computations is needed. 



Already in the case of three variables, the factorization problem of a corresponding 
operator and also constructing of its invariants becomes more complicated, even 
for constant coefficients. The reason of it is that in bivariate case we needed just 
to factorize leading term polynomial which is always possible over C. It is not 
the case for more then 2 independent variables where a counter-example is easily 
to find (see Ex.5), i.e. there exist some non-trivial conditions to be found for 
factorization of polynomials in more then two variables. 



3.6 Exercises for Chapter 6 

1. Using this particular solution u = —j^^^yji of the Liouville equation 

{\ogu)jcy + 2u = 
find general solution of this equation. 

Hint: Liouville equation is invariant under the following change of variables: 

2. Prove that matrix (3.21) is degenerate. 

3. Let Laplace invariants are equal, i.e. d = b. Prove that initial operator 

dxdy + adx + bdy + c 

is equivalent to the operator 

dxdy + c. 



4. Let in Lemma 3.8 function w{x, y) is chosen as 

w{x,y)^X,{x)Yi{y)+X2{x)Y2{y). 

Prove that d2 = 0. 



3.6. EXERCISES FOR CHAPTER 6 

5. Check that 

x^+y^ -3xyz = {x + y + z){x^ + y^+z^-xy-xz 
and prove that +y^ + is not divisible by a Unear polynomial 

ax + l5y + Yz + S 
for any complex coefficients a, j8, 7, 5. 
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Chapter 4 

Commutativity of linear operators 



4.1 Introduction 



In this Chapter, notion of commutativity for two linear operators will be studied 
which leads us to construction of LAX pairs consisting of a given linear operator 
and a new one commuting with a given operator. This technique is used afterwards 
for construction of a commuting operator for a linearized operator with coef- 
ficients depending on the solution of nonlinear equation. Existence of LAX 
pair for nonlinear differential equation is closely connected with symmetries of 
this equation and can be regarded as a definition of its integrability. Connections 
of this new definition of integrability with definitions given in Chapter 1 will be 
demonstrated as well as some examples. 



4.2 Operators of one variable 

4.2.1 Commutativity of formal series 

Commutativity of differential polynomials is a classical problem, highly non- 
trivial because polynomials constitute a ring. On the other hand, formal series 
constitute a field and their commutativity problem can be solved easily and nec- 
essary and sufficient conditions can be written out explicitly. Therefore, our first 
step here is to study commutativity properties for Lorant series which demands to 
generalize for the case n e Z formulae (introduced at the beginning of Chapter 1) 
for differential operator 
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^" • « = E ( ^ ) Z)^(fl)D"-^ (4.1) 



and binomial coefficients 

n\ n{n— l)...{n — k+ I) 



k J l-2---k \ 



with ( " 1=1. (4.2) 



A pseudo-differential operator (PSDO) is introduced [58] as a following formal 
series 

n 

A= ^kD''^a„D" + ---+ao + a-iD-'^+a-2D~^ + --- ■ (4.3) 

k=—°° 

while the product of PSDOs defined as consequent monomial multiplication using 
(4.1) and standard formula DW^ = D^^^ where j and k are not necessary positive 
integers, with the same formula (4.2) for the superposition with n e Z. 

For example, using for simplicity notation 

D{a) = ax, D^{a)=axx, 

we can compute 

o a = £ ( J D^{a)D-^-^ = aD'^ - la^D'^ + 3axxD~^ + . . . . 



It is important to understand clearly that notation Z)^(a) means action of operator 
on a in case A: > and it is just a formal description in case k<Q. 

It follows from definition that a product of two PSDOs of order m and n has order 
m + n and its leading term is product of two leading terms correspondingly. Asso- 
ciativity of so defined multiplication can be derived directly using the properties 
of binomial coefficients. 
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Definition 4.1 An element A given by (4.3) of associative algebra A e 
is called pseudo differential operator (PSDO) of order n with leading term 

Notation A+ is used for Taylor part of A, 

n 

A+ = £ aid" = a„D" + • • • + ao (4.4) 

A:=0 

and notation A_ is used for Lorant part of A, 

-1 

A_= £ aA:Z)^ = a_iD-i+fl_2D-2 + .... (4.5) 

That produces presentation of algebra ^S) as direct sum of subalgebras 

^^ = ^+{9)^^-9 (4.6) 

Obviously, in case when a PSDO under consideration coincide with its Taylor part 
A+, i.e. a formal series is in fact a finite sum, multiplication defined for PSDOs 
gives the same result as multiplication of usual differential operators. 



Lemma 4.2 Let A be a PSDO (4.3) of order n with 1 as leading coefficient. Then 
following keeps true: 

1. for arbitrary n e Z there exists unique formal series L such that AL = ZA = 1 ; 

2. for positive n > there exists unique formal series M of the first order, with 
leading coefficient 1 and such that M" = A and M"A = AM". 

< Starting with 

5 = D + Z7o + ^-ioD-i + Z7_2oZ)-2 + ... (4.7) 
one finds by induction that 

5" = + nh^U'-^ + Z7i,„D"-2 + Z72,„^""^ + . . . 
where the new coefficients 

hj^n = nb-j + f[bo, . . .,b-j+i]. 



98 



CHAPTER 4. COMMUTATIVITY OF LINEAR OPERATORS 



Hereafter /[...] means that the function / is a differential polynomial in its argu- 
ments. For example 

&i,„ = nb.i +f[bo] ^nb-i + "^"~^\ ^o,x + ^o)- 
The system of equations 

an-i=nbo, a~^ = bi^„, a„_3 = &2,n, • • • 
has the triangular form and therefore uniquely solvable. I 

Definition 4.3 Two PSDOs L and M obtained in Lemma 4.2 are called inverse 
and n-th root PSDO correspondingly. They are denoted as and aV«. 

Example 4.4 As an illustration let us compute powers and A^ of the first 
order PSDO of the general form 

A^D + ao+D^^ai +D^^a2+D^^a3 + ... (4.8) 

We find 

A^ = D^ + 2aoD + 2ai + al + flo,x + fl2l£>~^ + fl22£>~^ + • • • , 
A^=D^ + 3aoD^ + (3ai + 3al + 3ao,x)^ + 030 + a3iD~^ + + ... 

where 

a2i = 2a2 + ai^x + 2aoa\ , 
(230 = 3(22 + 3(2 i^;c + 6aoai +aQ + 3aoao^x + <^o,xxx- 

Formulae become longer with each step but in principle the procedure goes smoothly 
and can be continued up to any power. Representing n-th power of original PSDO 
A as A" = A^, -h A" , we can, using formulae above, compute Taylor's part for 
n = 1,2: 

A\_=D + ao, A% = D'^ + 2aoD + 2ai+al + ao^x + '2ai. 



For instance, the first three coefficients of the series A = y/Z for Schrodinger op- 
erator L = D^ — u (introduced in Chapter 2) have form 

2ao = 0, 2a\ = —u, 021 = 2a2 — —Ux = 0. 
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Proposition 4.5 Let A be a PSDO of the of the form (4.3) with 1 as leading coef- 
ficient. Then arbitrary PSDO B commutes with a given A iff fi can be represented 
as formal series 

5 = £j8,aV«, (4.9) 
where all are constants, Pk £ C. 



Definition 4.6 Pseudo-differential operator C = AB — BA constructed from two 
PSDOs A and B is called their commutator and it is denoted asC = [A,B]. 

Proposition 4.7 Let A is a PSDO of order n and 5 is a PSDO of order m, then 
their commutator has order k <m + n, 

[A,B]= CkD\ 

k<m+n 

and its leading term has form 

In particular, by a„ = 1, n 7^ and [A, 5] =0 this formula implies that the series 
B has constant leading coefficient bm^x — 0- 

4.2.2 Commutativity of differential operators 

Ordinary differential operators 

A = a„D" + an- H h a\D + 

can obviously be regarded as a particular case of PSDOs and the operation of 
multiplication defined in previous section, can be understood in usual sense, as 
composition of ordinary differential operators. Also result of Prop. 4.7 keep true, 
i.e. if A and B are differential operators of orders m and n respectively, then their 
commutator C is a differential operator of order <m + n— I. 

On the other hand, in comparison with PSDO, differential operators have some 
specific properties which we are going to present below, beginning with following 
definition. 
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Definition 4.8 Centralizer C(A) of an ordinary differential operator is a set dif- 
ferential operators {B} which commute with A, i.e. [A,B] = AB — BA = 0. Cen- 
tralizer C(A) is called trivial if it contains only polynomials in A with constant 
coefficients. 



In order to construct classes of equivalence on the set differential operators, let us 
introduce two transformations of differential operators as follows: 

• change of independent variable x — > x 

d d dx , , , , , ^ 

a— = —^ — = a{x) (4.10) 
ax ax ax 

• similarity transformation 

A^A = f-^oAof (4.11) 
where a = a{x) and / = f{x) are appropriate chosen smooth functions. 



It is easy to see that transformations (4. 10) and (4. 1 1) allows bring the differential 
operator A of order n > to the canonical form 

A^D" + an-2D"~^ + --- + ao (4.12) 

in which the first two coefficients are fixed by 1 and 0, respectively. Particularly, 

1 /n 

the change of variable (4.10) with a = leads to the operator with unitary 
leading coefficient. An example of the differential operator in the canonical form 
is given by Schrodinger operator L = — u. Any differential operator of second 
order is equivalent Schrodinger operator L = Z)^ — a up to transformations (4. 10) 
and (4.1 1) as has been indicated above. 



Example 4.9 Let us describe centralizer for first order differential operator A = 
a{x)Dx + b{x). After change of a variable (4.10) it takes form 

^ C4' 

A=D + b,D=- 
ax 

and the second transformation another change of variables 

A = a~^Aa, with D\oga = b 
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gives finally 

A = D 

where all "tildes" and "hats" over A are omitted for simplicity of notations. It can 
be shown that both transformation (4.10) and (4.1 1) does not change structure of 
centralizer and Proposition 4.5 gives its description as 

{Z)0, D,d\d\...}. 

Thus, centralizer constructed in this example, is trivial. 



Theorem 4.10 Let C(A) is centralizer of a differential operator A of order 
n > 0, then any two elements of the centralizer, 8^,82 EC{A) commute with each 
other: 

[5i,52]=0. 

Proof. As in examples above, an appropriate change of independent variable pro- 
duces constant leading coefficient of A and without loss of generality let us put 1 
as leading coefficient. Then statement of the theorem is direct corollary of Propo- 
sition 4.5. 1 

It follows from Theorem 4. 10 that 

[Ai,A2]=0, [A2,A3]=0 ^ Ai,A3eC(A2) ^ [Ai,A3]=0. 

Therefore the binary relation [A, 5] =0 is an equivalence and one can replace 
the operator A in C(A) by any nontrivial (of non zero order) element lying in 
centralizer. The minimal order n > of nontrivial elements of C(A) gives some 
indication what is a structure of the centralizer. Thus, for n = 1 a centralizer 
is always trivial as it was demonstrated in Example 4.9. But in the case with 
n = 2 it is not true. In order to construct full description of nontrivial centralizers 
in the case n = 2 we need to study first a problem about commutative pairs of 
operators of the second and the third orders. In this case we need just compute the 
commutator C = [A, 5] and write down the condition when it vanish. 

Example 4.11 Let A = - m and B = D^ + aD + b. Then 

[A, B] = axD^ + a\D + ao, 
a2 = 2% + 3m;c, ai = axe + 3mxc + 2^x5 ao = "xtc + + ^xc- 
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Excluding coefficients of the third order operator one gets equation for the poten- 
tial u 

Uxxx — 6uUx = £Ux (4.13) 

where e e C is a constant of integration. Any solution of Eq. (4.13) gives rise a 
commuting pair of operators. Particularly, in the case 

A = D^+x'~^e, B = D^+x~^aiD+x~^a2 

the condition of commutativity [A,B] =0 reduces to linear algebra as follows 

3e = 2ai, ai + a2 = 0, 12e + aie + 6a = =^ e = -2, ai = -3, a2 = 3. 



Thus, we have now the operator A = — 2x^^ of the second order such that C(A) 
contains the third order operator B. Since all polynomials in A gives operators of 
even orders, the third order operator B can not be written as polynomial in A and 
thus we have at last example of nontrivial centralizer C(A). One can verify di- 
rectly that A^ = B^ in this case (Cf. Prop. 4.5). 

This way it was shown that any nontrivial centralizer C(A) of Schrodinger operator 
A — — u contains some odd order operator B. In the opposite case, when all 
orders of operators in centralizer are even one can easily prove that any operator 
B eC{A) of order 2n could be represented in the polynomial form 

B = CnA" + Cn-iA"-^ + ---+co, G C. 

Indeed, the leading coefficient should be constant due to Proposition 4.7. Denote 
this coefficient by c„ and consider the operator B = B — CnA^ e C{A) . By condition 
its order is even 2n, n<n again and we can proceed further in order to diminish 
it. 

Analogous reasoning in the case of nontrivial centralizer C{A) of Schrodinger 
operator leads to following algebraic formula for any operator ^2 G C(A) : 

bI = P2{A)Bx + Q2{A) (4.14) 

where B\ e C(A) has minimal odd order and PiiA)^ Qii^) polynomials in A 
with constant coefficients. Particularly, apply this formula to B\ we get 



5? = Pi(A)5i + (2i(A),5 = 5i-^Pi(A) ^ B^ = Q{A). 
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It is easy to see that new operator B e C(A) has the same minimal odd order 
m = 2n + l as Bi and can be used as generator as well as 5 1 in the general formula 
(4.14) and, hence degree of polynomial Q{A) in the last equation is the same m = 
2n+ \ . In the next Section we describe how to reconstruct centralizer using this 
polynomial Q{A) which defines in certain sense multiplication in the commutative 
ring C(A) with two generators A and B related by the equation 

B^^Q{A). (4.15) 

For this construction we need to introduce a very important notion of common 
eigenfunctions for two commuting differential operators L and M of orders n and 
m, respectively. Let us notice first that if 

Lx]/ = Xx]/, and [L,M] = 0, 
then for some function x(f such that iff — Mxj/ we have 

L\j/ = Ai/A. 

It means that for any fixed X the n— dimensional subspace Ker(L — A) constituted 
by solutions of equation Li// = A l/A for eigenfunctions of the operator L, is invariant 
under the action of operator M. 

Definition 4.12. Let L and M are commuting differential operators and M is the 
restriction of the operator M on an invariant subspace Ker(L — A ) . The eigenvector 
Mxj/ — iixj/ is called common eigenf unction of two differential operators L and 
M. 

It follows from Definition 4.12 that this common eigenfunction \j/ solves two 
equations: 

Lxj/^Xxi/, My/ ^ (4.16) 

In particular, for Schrodinger operator L = A = D^ — u and M = B given by (4. 1 5), 
the last formula above implies ju(/l) = ^/Q{X) and therefore Eq. (4.16) takes the 
form 

\j/xx={X+u)\j/, By= ^/Q{X)\j/. 

It could be shown that this over determined system of two equation for one func- 
tion y/ (common eigenfunction) can be reduced to the first order linear ODE (due 
the fact that greatest common divisor of orders of operators L and M is equal 1.) 
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Example 4.13 The differential operators from Ex. 4. 1 1 are related by the equa- 
tion = and equations (4.16) for common eigenfunction give 

2 3 3 

After differentiation with respect to x the first equation and comparison it with the 
second one, we obtain the first order ODE: 

Wx X^k^ + l , 1 

X— = — kx-\ 

y/ x^k^ - 1 kx-1' 



Remark. In many cases it is convenient to rewrite the scalar differential equa- 
tion Li// = A i/A for eigenfunctions as the systems of first order equations for the 
vector function y '• 

xirx = V{x,X)xir (4.17) 

For the general differential operator L of order n it is n— th order system with 
\(f= {\i/^,...,\i/"y. For instance, for Schrodinger operator we have 

and therefore the 2 x 2 matrix V for Schrodinger operator is 

It is possible to proceed with the same reasoning as above using matrix represen- 
tation of operators and generalize the results obtained for Schrodinger operator on 
the operators of higher order but this task lies beyond the scope of this text. 



4.3 Operators of two variables 

Above we regarded operators of one variable that implies a differentiation Dx with 
respect to single independent variable x. Now we try to insert additional differ- 
entiations D;^, compatible with Dx, and to introduce a corresponding set of inde- 
pendent variables ti^, k= 1, 2, 3, ... . First, we regard this problem algebraically 
and consider action of these additional differentiations A^, upon elements of the 
ring of formal series in powers of D. It immediately produces a nonlinear system 
of PDEs for coefficients of the formal series. In particular these PDEs include in 
itself famous Kadomtzev-Petviashvili equation (KP): 



^tx — {^xxx 6mMx)x ~l~ 3m^. 



(4.19) 
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Further we consider a broad class of evolutionary equations 

Ut=f{x,U,Ux,Uxx,Uxxx---) (4.20) 

with one space variable x. We show how to build up constructively a Lax pair for 
equations (5.2) possessing infinite number of higher symmetries. It is remarkable 
that in this case basic theorem is formulated again in terms of formal series in 
powers of D and highlights notions of algebraic theory of the first sections from a 
new point of view. 

4.3.1 KP-model 

In order to introduce an action additional differentiations acting upon elements 
of the set of PSDO (4.3) it is sufficient to define it action on the coefficients of 
formal series A = Y.o-kD'^- We postulate the formula as follows 

A, (A) = fl„,f,Z)" + • • • + ao,f, + a-x,t,D-^ + . . . [B^, A] (4.21) 

where the formal series 5^, k= 1, 2, 3, ... will be defined exactly below and 
n e Z denotes an order of A. 

Firstly, we remark that the basic differentiation D in the set of formal series could 

def 

be written down in this form and corresponds to the choice B\ = D. Namely, use 
the definition of multiplication of the formal series A = and D we have 

[D,A] =DA-AD = Y,Do a^d" - a^d"^^ = £afc,x£>^ (4.22) 

Thus, we verified Eq.(5.3) for A: = 1 with B\=D and = D^. 

Secondly, we remind that the Leibnitz rule (general feature of any differentiation) 
is satisfied for the operation of commutation with a formal series B since 

[5,AiA2] =5AiA2-AiBA2+Ai5A2-AiA25= [5,Ai]A2+Ai[5,A2]. 

Thus, left and hand sides of Eq.(5.3) obey Leibnitz rule and we wish now to 
compare orders of the two formal series in the formula under discussion. Due 

Proposition 3.7 the order of commutator [5, A] is equal « + m— 1 if the formal 
series B of order m is chosen independently from the series A of order n. On the 
other hand the order of A^ (A) has to be less or equal of the order n of A. In the case 
m = 1 all agreed since n+m — l = n and that opportunity used in Eq.(5.4). For m > 
1 in order to equalize the orders we must choose at least leading coefficients of Bj^ 
in the strong accordance with leading coefficients of series A. It is the source of 
nonlinearity imprinted in Eq.(5.3). Take in account Proposition 3.5 we formulate 
below very nice algebraic solution of this problem of equalizing orders in left and 
right sides of Eq.(5.3). Some times next proposition is considered as a model of 
modem integrability theory. 
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Lemma 4.1. (Sato ??) Let formal series B of the first order is normalized as 
follows 

B = D + biD^+b2D^ + b3D^ + ... (4.23) 

and the differentiations (5.3) of the coefficients of this series defined by formu- 
lae 

Dk =^Dt^{B) = Z7i,,,Z)-i +b2,t,D-^ + ■■■ = (4.24) 

Then the differentiation Di=Dx and DnD^ (B) = DmDn (B) for any n, m> 1 , and 
moreover 

Dn{B^)-Dm{B\) = [BI,B1], DniB^l) -Dm{Bi) = [B1,B'\_]- (4.25) 

M Obviously B+—D and, therefore, Di (B) = Bx- In order to prove (5.7) we start 
with the identity 

= = [Bl,B"l] + [Bl^B"!^] - [B^l^Bl] + [Bl^Bl]. 

On the other hand 

[D„{Bl)^Dm{B"+' -B1) = [BlX-] + [B"l,B"+] -DmiBl), 

Dn{Bl)=Dn{B^+'-B^) = [Bl,B^] + [Bl,B"l] -D„{Bl). 

and (5.7) follows from it. It is easy to see now the commutativity conditions 
DnDm{B) = DmDn{B) are corollary of Eq. (5.7). I 

Remark 4.2 In this chapter we do not consider additional differentiations Df^ for 
formal series of general form. Nevertheless, the equations Eq.(5.6) and Leibnitz 
rule imply that for powers 5" of the normalized series (5.5) we have 

A,(fi") = [<,5«]. (4.26) 

This generalization Eq.(5.8) and Proposition 3.4 allows to include in domain of 
differentiations formal series (4.3) of any order n but yet normalized (i.e. lead- 
ing and next one coefficients equal 1 and 0, respectively). In order to get off these 
restrictions one should to develop algebraic variant of transformations (4.10 and 
(4.1 1 used in previous Section yet it is not easy at all. We are going to discuss this 
normalization problem in next section from quite different point of view. 

The nonlinear action on elements of of additional differentiations Df^, com- 
patible with Dx defined by Lemma 4.1 is not unique. In particular, one can start 
with the formal series 



A = D + aQ + a\D ^+a2D ^ + 03!) ^ + . 



(4.27) 
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and use another way (Cf. (4.4), (4.5)) to separate polynomial and Lorant parts of 
A" : 

A"=Al+A^, where =^ D" + an,iZ)"-^ + a„,2^""^ + • • • + an,n-iD 

and n = 1,2,3,.... The proof of next proposition repeats the proof of Lemma 4. 1 
and will be omitted. 

Lemma 4.3 Let formal series A of the first order defined by (5.9) and the differ- 
entiations (5.3) of the coefficients of this series defined by formulae 

Dk=^Dt,{A) = ao,t,+ai,t,D-^+a2,t,D-^ + ---=[Al,A]. (4.28) 
Then Di = Dx and DnDm{A) = DmD„{A) for any n, m>l. 

One could notice that as well as Lemma 4.1 corresponds to representation 

as the direct sum of the polynomial and not polynomial parts (see Definition 3.1) 

Lemma 4.3 corresponds to another way do define this direct sum (Cf. (4.6)). 

We are going now to consider the action of differentiations D^, on elements 
in more details and write down an explicitly some nonlinear equations for coeffi- 
cients of the formal series (5.9) and (5.5). In order to stress that we take next: 

Definition 4.3 We call by mKP-hierarcJiy and KP-hierarchy the infinite sys- 
tem of PDEs for the coefficients and of the formal series(5.9) and (5.5) 
defined by Eq.(5.10) and Eq.(5.6), respectively. 

Below we provide basic, in certain sense, examples of equations of mKP and 
KP hierarchies. It will be equations for uq and bi , respectively, with derivatives 
on the first three independent variables h,t2,t2. In the case of mKP-hierarchy it 
called modified Kadomtzev-Petviashvili equation (mKP shortly) and is similar to 
KP-equation (5.1) mentioned in the introduction. 

Example 4.5 Historically, for KP-equation (5.1): 

Uxt = {Uxxx - ^UUx)x + 'iUyy 

the auxiliary linear problem arised in the form as follows 



¥y = ¥xx-u\i/, e¥t = ¥xxx + a\i/x + b\i/ 



(4.29) 
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In order to find the compatibility conditions which gives KP-equation we have to 
equate cross derivatives {\i/y)t = {¥t)y Denote xj/k = D^^v^q find after x— differentiation 
Eqs.(5.11) 

e{Wy)t = V^5 + {aW\ + b\if)xx -eutW- «( V^3 + axj/i + bij/), 

e{\lft)y = V5- (My)xxx + %V^i +a(v/3 -MV^i -MxV) + byXif + b{\ir2 ^ UX^f) . 

Equate now coefficients by ^fk in above formulae one gets for A; = 2, 1, 0, respec- 
tively 

lux + 2>Ux = 0, 2bx = ay- 3uxx -axx = ay + Uxx, 
by + eut = bxx + 
Rewrite the last equation in terms of u one obtains KP. 

In terms of Definition 4.3 it corresponds to differentiations (5.3) with ^ < 3. We 
have for 1 Dt^ : 

mKP: A%=D^ + 2aoD, A^ = + 3aoD'^ + 3{ai +ao^x + al)D. 
KP: Bl=D^ + 2bi, Bl = D^ + 3biD + 3{bi^ + b2) 

Therefore it coincides in KP-case with (5.11) up to notations. Moreover, Lemma 
4.1 equations 

D2{B) = [Bl,B], D3{B) = [Bl,B] 

yield 

bi,y = bi^xx + 2b2,y, b2,y = b2,xx + 2b3^x + 2bi^xbi. 
bi,t = bi^xxx + ^b2,xx + 3&3,x + bibi^x 
and, thus, exclude b2 and b^ we find: 

1 3 

[bi,t - ^bi^xxx - ^bib\^x]x = -^bi^yy 

Compare these two view points on KP one could say that these are equivalent and 
give the same nonlinear PDE. 

Analogously, using Lemma 4.2 equations one gets mKP 

%x = Dxiqxxx - 2ql) + 3qyy + 6fc^y, uq = qx- (4.30) 

In algebraic approach to integrability the notion of commutativity appears prefer- 
able in comparison with the notion of symmetry. Thus it seems awkward to re- 
formulate Lemmas 4.1-2 in terms of symmetries of the infinite system of PDEs 
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constituting hierarchies. Yet consider particular equations of hierarchies like (5.1) 
and (5.12) the notion of symmetry becomes more important and useful. In the 
next example we demonstrate reductions Eq.(5. 1) related with stationary points of 
very simple symmetries transformations. 

Example 4.6 Since the equations contain independent variables only inexplic- 
itly they are keep its form invariant under shift In the case bi not depends in y KP 
equation reduce and gives rise KdV: 

bl,y = ^ 

Bussinesque?? 

4.3.2 Conservation laws 

Interplay of symmetries and conservation laws for ODEs has been discussed in 
Introduction. This interconnection for PDEs is quite different but very interest- 
ing as well and we present here a terse introduction to pure algebraic aspects of 
the theory oriented on applications in next chapters. Hereafter in this section we 
restrict ourself by the case of KP-hierarchy and will denote by ^ the set of poly- 
nomial functions on coefficients bj of the formal series (5.5) and its derivatives 
with respect to / = 1,2, 

Definition 4.7 Corollary of equations of KP-hierarchy written in the form of 
divergence i.e. 

EA,(i8,) = (4.31) 
k 

is called conservation law of KP-hierarchy if sum in (5.13) is finite and the the 

variables j5k e ^. 

In order to build up conservation laws of KP-hierarchy one can use lemma below 
which based on next definition. 

Definition 4.8. For the formal series A e in powers of D 

n 

A= akD^^anD''-\ hao + a_iZ)"^ +fl_2^"^ + . . . . 

the coefficient a_i is called residue and is denoted as res (A). 
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Lemma 4.9 For all m, n e Z 

ros{[aD'",bD"])=D,a^^„ 

where „ is a polynomial function on a, b and it's derivatives. 

Proof. The residue vanish if the powers m, n obey the condition mn > 0. For 
instance in the case n — the commutator aD'^b — baD^ is a differential operator 
if m > and PSDO of order m— 1<— 2ifm<0 (see Prop. 3.7). Obviously the 
coefficient by is zero in both cases. Obviously as well that the residue vanish 
if m + n < 0. 

Let now m, n have different signs and m + n = k>Q. Then in virtue of (4.1) 

res([aD'«,Z,D«]) = ( ) {aD''^\b) + {-\fD'^+\a)b) 

since 

m + n = k =^ m{m — \)---{m — k) = ±n{n — \)---{n — k). 

Standard "integration by parts" proves now the lemma. Particularly for ^ = 0, 1 
we have, respectively 

aD{b) + D{a)b = D{ab) aD^{b) - D^{a)b ^ D{aD{b) - D{a)b). 

I 

Use (5.8) we find 

A^(res5") = res Df^B" = res[5* 

and Lemma 4.9 yields 

Corollary 4.10. KP-hierarchy possess infinite sequence of conservation laws of 
the special form as follows 

Dt,{pn)^D,{Gk,n). p„ = res(5"), ak,n e ^ (4.32) 

where ^ denotes the set of polynomial functions on coefficients bj of the formal 
series (5.5). 

Example 4.11. 

Pi=^i, P2 = 2b2-bi^, P2 = 2b2-bi^, p3 = ^b3 + ^bj-b2^-bi^ 
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Theorem 4.12 ([?]). Let in the case (5.5) 

oo oo 

A = D + Y^biA-', A"^=iA"^)++Y^bkmA-K 

i=l k=l 

Then for all /, n = 1 , 2, . . . 

Dnbi = Dbin, {D = Di). (4.33) 

Proof. Lemma 4. 1 yields 

oo oo 

D„{Ai) = D{An) + [Ai,A„], Ai = J^M"', A„=Y, bumA'^ 

i=l k=\ 

We have now 

CO oo CO oo 

Dn{Ai) = Y.^n{bi)A-' + Y.^i[A-\An] D{An) = D(Z7fe„)A-^+ £ bkn[A'\A,] 
i=l (=1 k=l k=l 



and 



[Ai,A„] + bkn[A-\Ai] - £Z7,-[A-',A„] = 0. 

k=l i=l 



In order to highlight the notion of conservation laws let us consider 7— periodic 
in X solutions of KP-hierarchy. Then 



rxo+T rxo+T 

Dt^ip) = D^iOk) A, / pdx= Ok^^dx = Q. 



Therefore, the integral R{t) = J^^^ pdx conserves it's value as times tk change 

and the integral R = /^^ pdx rather than density p has "physical sense." That 
is the reason to call two conservation laws in the differential form A^(py) = 
Dx{Oj^k)i j = 1,2 with distinct densities equivalent pi = p2 if 



fXo+T rxo+T 

/ pidx— / P2dx. (4.34) 

Jxo Jxq 



In general not necessary periodic case, we take 
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Definition 4.11 Two conservation laws of KP-hierarchy D;^ (pj) = Dx{<Jj^k) i j = 
1,2 to be said equivalent iff Pi — P2 e where latter is subspace of ^ consti- 
tuted by ImD. 

Obviously, equivalence relation introduced in Def.4. 1 1 implies the equality (5.16). 
Moreover, as we will see in the next section this relation p = could be verified 
in fully algorithmic way. 

4.4 Integrability problem 

In this section we discuss a general definition of integrability for evolutionary 
PDEs (5.2) with one space variable x using the theory of formal series in powers 
of D = Dx (sections 2 and 6). Follow traditions of differential algebra [63] we 
introduce the notation 



for the infinite sequence of derivatives with respect to x of the basic variable u. 
In these notations the general evolutionary PDE is rewritten as follows 



The highest derivatives m„ defines the order n of Eq. (5.17) and we consider 
usually equations of order n>2. The dependence on x of the right hand side in 
Eq.(5.17) became very important in the next Chapter in relation with godograph 
type transformations. 

4.4.1 Basic definitions. 

Consider the set ^ of smooth functions of variables a:,m,mi,M2,--- we will use 
the equation Eq. (5.17) in order do define action of the "additional" differentia- 
tion Dt on ^ (the first one is D = Dx). Namely we take as definitions of these 
differentiations the formulae 



M l-> Ml !->■ M2 '->^ •••'->■ Mfe •• • 



Ut= f{x,U,Ui,U2,...,Un). 



(4.35) 



Dx'. g{x,U,Ul,U2,...,Um)^ gx + g*{ui), Ul=D{u), 
A: g{x,U,Ui,U2,...,Um) ^ g4f)^ f = Dt{u) 



(4.36) 
(4.37) 



In fore cited book it called differential indeterminate. 
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where gx = dx{g) is partial derivative of the function g and g* denotes differential 
operator as follows 

^ OUk OU OU\ OUm 

It is easy to see that operators (5.18),(5.19) act on the "generators" of the set 
in a natural way, as partial differentiations should do: 

Dx{x) = \, Dx{u) = ui, Dt{x) = 0, Dt{u) = f{x,U,Ui,U2,...,Un) 

and 

DxDt{u)=Dx{f), D,Dx{u)=Dt{ui)=Dx{f) ^ DxDt{u) = DtDx{u). 

It is important to notice that action of m— th order differential operator g^ on func- 
tions V can be defined also as 

,.(v) = !?«!^|,.„. (4.38) 

In particular as corollary of Eq. (5.20) we find that 

{fg)* = fg* + gf*- 

One can prove now that Djc, Dt and A = DxDf —DfDx as well, obey Leibnitz rule 
and that A = on the set ^ . 



Definition 4.13. The operators (5.18) and (5.19) on the set are called total 
differentiation with respect to x and evolutionary differentiation related with 

/ e respectively. The Gato derivative (5.20) of a function g e we call 
linearization of the function g. 

Summing up we have 



Proposition 4.14. For any smooth function f & ^ corresponding evolutionary 
differentiation (5.19) commutate with (5.18). 

This statement looks trivial but if, as in previous Section, we try to find other 
evolutionary differentiations 

Dt^{u)=fk{x,u,ui,U2,...,u„^) (4.39) 
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which commutate each other we shall come to nontrivial obstacles. Situation here 
in certain sense corresponds to correlation of the Proposition 3.5 which is very 
general and abstract and more interesting problem about commutating differen- 
tial operators considered in Section 4. The role of Proposition 3.5 plays now 
Lemma 4.1 and commutativity conditions — have to play the role of 

commuting differential operators. This link became more apparent and instruc- 
tive if we bring up commutativity conditions to the level of a tangent space of 
using (5.20). Namely instead of writing down the cross derivatives equality 
D^Dt^iu) — Df^Dt^ (m) for the two evolutionary equations 

Dt^{u) = fl{x,U,Ul,U2, . . . ,Urn), Dt^iu) = f2{x,U,Ul,U2, . . . ,Un2) (4.40) 

we will consider analogous condition Dt^Dt^iv) — Dt2D^{v) for linearized ones 
A;t(^) ~ ifk)*{v)i k= 1,2. Translate this condition into condition of commuta- 
tivity of corresponding operators 

[Ai-(/i)*,A2-(/2)*]=0 (4.41) 

we come to next ^ 

Definition 4.15. Let us say that the two evolutionary equations (??) satisfy the 
consistency condition and, respectively, compatibility condition if the cross 
derivatives equation Dt^Dt^iu) — Df^Dt^ (u) is fulfilled or, and in the second case, 
if differential operators equality (5.23) is valid i.e. 

Ai((/2)*) - A2((/i)*) = [(/i)*, ifiU (4.42) 

We remind (Cf. (5.3)) that the action of any evolutionary differentiations Dt,^ on 
polynomials A = Y,cijD^ ■, cij & ^ is defined by Dt,^{A) = Y,Dt^{'^j)DK 

In order to compare two almost equivalent Definitions 4.15 we give simple exam- 
ple. 

Example 4.16. Let us consider the pair evolutionary equations as follows 

Ai (w) = fix-, U,Ui,U2,...,Un), (u) = Ml . 

Then, apply (5.18) and (5.19) and Proposition 4.14 we get 

DtMu)=Dt,{ui)=D,{f)=f, + Mui), DtMu)=Dt2{f)=Mui)- 
^It can be proved that cross derivatives equation Dt^Df^ {u) = Df^Dt^ (m) implies (5.23). 
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Thus the consistency condition is satisfied iff fx = and / = f{u, ui,U2, ■ ■ ■ ,Un). 
On the other hand the compatibility condition (5.24) yields 

fu,x = fuux = ■■■ = fu„,x = ^ f = f{u,Ui,U2,...,Un)+X{x) 

where is arbitrary function of the independent variable x. 

4.4.2 Main theorem 

We wish now to address more practical question. Let the evolutionary equation 
(5.17) of the order n is given i.e. A(m) = /• What conditions on its right hand 
side f or, more exactly, on the corresponding differential operator : 

= ^ + + + (4.43) 

OU OU\ OUn 

are needed in order to build up a hierarchy of evolutionary differentiations (w) = 
fk compatible with fore given. Take in account that in the compatibility condition 

A((A)*) -DM),) = [/., im (4.44) 

the order of the differential operator {fk)* could be arbitrary large whilst the 
order of Dty.{{f),) is fixed by n we see that leading terms in Dt{{fk),) and in 
the commutator in the right hand side should coincides. It gives rise the definition 
which reproduces the main formula (5.3) of previous section in a new environment 

Definition 4.16 We say that evolutionary equation 

Mf =/(x,M, Ml, M2, •••,««) 

of order n posses formal symmetry if there is the formal series of order m^O 

A = aniD"' + Um-iD"^^^ H h ao + U-xD^ + a-2D^ + . . . 

with coefficients aj e ^ such that 

A(A) = [/.,A]. (4.45) 



Apply this definition to Burgers and Korteweg de Vries equations which we often 
use as illustrative examples, we find 
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Example 4.17 For Burgers equation 

Ut = U2-\-2uu\ — f , f*=D +2uD + 2ui 
the formal symmetry is PSDO as follows A = D + u + u\D~ ^ . Indeed 

[/^,A] =/ + Z)(/)Z)-i and A ^D + m + miD"! ^ = "f + "1,?^"^ 

and thus A (A) = [/*,A]. 

Similarly, in the case of Korteweg de Vries equation (KdV) 

Ut = duui — U'i=f, /* = —D + duD + 6mi 
the second order PSDO A = D^ —Au — 2u\D~^ is the formal symmetry since 
[/^,A] = -4/-2Z)(/)Z)-\ and A^D^-Au-2uiD-^^Dt{A)^-Aut-2ui^tD~^ 
and thus A (A) = [/*,A]. 

Next statement can be considered as milestone of integrability theory of evolu- 
tionary equations with one space variable. 

Theorem 4.18 If there are series of evolutionary differentiations Dtf,{u) = fk of 
orders ^ o° hy k ^ oo satisfying compatibility conditions (5.26) for a given 
evolutionary equation A(") = / of order n>2. Then this equation Df{u) = f 
possess the formal symmetry A of the first order: 

A = aiD + ao + a-iD^^ + ... (4.46) 

Construction of the solution A of the basic equation Dt{A) = [/*,A] needs some 
technicalities related with the limiting procedure in Eq.(5.26) when k —> oo. We 
skip that formal part of proof of Theorem 4.18 (see Appendix??) and go instead 
directly to applications. 

Firstly, we explain how to get formal symmetry of the first order starting with any 
solution of Eq.(5.27) and describe it's general solution. Obviously, 

aiAi + a2^2, ,«yec and A1A2 
satisfy Eq.(5.27) if Ai and A2 do. Apply results of Section 3 we get. 
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Lemma 4.19 Let formal series in D of order m i-e- B = Y^bjD^ bj e ^ is 
solution of Eq.(5.27) and let m 7^ 0, n > 2. Then there exist the first order formal 
symmetry A (see (5.28)) with the leading coefficient 

/ „ N i ^ def of 

and any solution B of Eq.(5.27) of order m can be represented as formal series in 
powers of A with constant coefficients a/ e C: 

Proof. Slight generalization of Lemma 3.2 yields the formal series Bi with coeffi- 
cients from 'W such that 5^ = B. Let us consider the formal series R = Df{Bi) — 
Since 

Dt{B) = A(5i)Br' +5iA(5i)5r^ + " " " +5r'A(5i), 

and 

= [f,^A]A'"-^+A[f,,A]A"'-^ + ---+A'^-%A] 

we have 

5f - = ^""i +AiM'"-2 + • • • + A'"-ii? = 0. 

In this sum the all m series in have the same leading term and, hence, this 
leading term is zero. It is possible only if ^ = and therefore we proved that Bi 
is a solution. 

In order to find the formula for leading coefficients of solutions Eq.(5.27) one 
should find the leading term in —Bt. Since n > 2 it is contained in 

= Cn+m-lD"'^'"~^ + ■■■, Cn+m-l = nfnD{bm) - mbmD{fn) 

where /„ and bm are leading coefficients in and B (Cf Proposition 3.7 ). Van- 
ishing leading coefficient gives the equation 

nDlogbfn = mDlogfn <^ nlogbm — mlog/„ = const 

which implies for m—l the formula (5.29) for the leading coefficient of the first 
order formal symmetry. I 

Most known applications of Theorem 4.18 related with next proposition quite 
similar to Corollary 4.10: 
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Corollary 4.20. In conditions of Theorem 4. 1 8 the evolutionary equation (5.17) 
possess infinite series of conservation laws as follows: 

A(resA'')=resA(A")=res[/*,A"] elmD, n = -1, 1,2, . . . . (4.48) 

Hereafter ImD denotes the subset of functions g & which can be rewritten as 
g = Dx{h),h e In other words Corollary 4.20 provides series of conservation 
laws 

A(Pn) = Dx{an), pn: CT^ G ^, p\ =^res£/V 
of the equation (5.17)) under cosideration. 

Remark 4.20'. Series conservation laws (5.30)) should be completed by the 
conservation law with the dencity po = 

Namely, rewrite equation (5.27) for m— th order formal symmetry 

A = amJy^ + am- iD""' hao + a_iD"^+a_2D"^ + ... 

as follows 

A(A)A-i = /, - A/,A-i = [A-\Af,] 
we find res (A;A~^) e ImD due Lemma 4.9 and that 

res(AfA-i) = res [{am,,D'" + am-i,,D"'^^){am^D-"' + aD-'"-^)] = 
{am,t{oc-m^2r)^ ^• 

It remains to notice that 

am,x ^m— 1 

a = . 

dm 

Example 4.21. For n = — 1 we have due (5.29) 

resA-^ = res^^£)-i + ...^ =(/«)"»• 

Thus p_i = (y^^ " and by Corollary 4.20 the evolutionary equation (5.17) 
should always satisfy the condition 

A(P-i)=a(^^^ "elmD (4.49) 
if it is compatible with hierarchy of evolutionary differentiations. 



4.4. INTEGRABILITY PROBLEM 



119 



Definition 4.22. The conservation law (i.e. differential corollary of (5.17)) 

Dt{p)=D,{a), p,a^^ (4.50) 
is called trivial if p e ImD and there is /? e such that 

p=D,{h) Dt{p)=Dt{D^{h))^D^{D,{h))eJmD 

for any evolutionary differentiation Dt. The functions p and a in Eq.(5.33) are 
called density and flux, respectively, and two conservation laws with the densities 
Pi and P2 are considered equivalent if Pi — P2 e ImD. 



Come back to Example 4.21 we see now there is an alternative or (P-i) G ImD or 
A(P-i) G ImD. In both cases this condition is non trivial one and can be verified 
using classical statement by (Helmgolz??) proven in P.Olver book [59]. 



Theorem 4.23. The function h e belongs to the subspace ImD iff E{h) = 
where Euler operator E is defined as follows 

- def d d T d ^ d 

£4E^_-D— +d2— -d3— + ... . (4.51) 

ou ou\ au2 du-i 



Equations similar to Eq.(5.27) play major role in modem integrability theory and 
are called commonly by Lax pair. The first test which have to justify that that is 
not trivial one ^ could be formulated similarly to the definition 5— integrability be- 
low. In other words proper Lax pair must produce series of non trivial conservation 
laws. 



Definition 4.24. The evolutionary equation (5.17)) with formal symmetry (5 . 27)) 
is called 5— integrable if the series (5.30)) contains nontrivial conservation laws 
of arbitrary higher order and C— integrable in the opposite case. 

The famous Korteweg de Vries (see (2.163) in Section?) gives us an example 
of 5— integrable equation. As matter of fact the formal symmetry for KdV indi- 
cated in Example 4.17 related very closely with Lemma 4.6 ( see Eq.(2.89) and 
concluding remark in Section?). KdV hierarchy will be discussed in next section. 

Although Definition 4.24 looks too complicated to be really useful it is not the 
case. In the next chapter we apply this definition to the second order evolutionary 



^There is known a few cases of faked Lax pairs. 
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equations (5.17)) and verify that first three of conditions (5.30) fully describe all 
second order evolutionary equations possessing hierarchy of higher order symme- 
tries. Burgers equation from Example 4.17 is, in certain sense, typical representa- 
tive in this class and a class of C— integrable equations in general. 



Symmetries and cosymmetries 

Concluding this chapter we discuss duality between symmetries and conserva- 
tion laws starting with "zero" dimensional evolutionary equations i.e. dynamical 
systems 

^ =/(«), u={u\u\...yf (4.52) 
at 

where ()^ denotes hereafter matrix transposition. In this case the conservation 
law with the density p = p{u) is just a constant of motion 

at au^ ou" 

Denote by p the column- vector constituted by the partial derivatives ^ we notice 
that 

A+/Jp = 0. 



4.5 Summary 

The primary goal of this chapter is to establish a workable criterion that can be 
readily checked to determine whether a given evolutionary equation possess inte- 
grability features like symmetries, local conservation laws and Lax pair. 



4.6 Exercises for Chapter 4 

1. Prove the set of DO of the form L = D" + Y.^k^^^D"^'^ closed under multi- 
plication. Define the transformation which kill the coefficient Ci in L. 

2. Find the fist three coefficients in the series A = (D — f)^^. Prove that all 
coefficients are "bona defined" i.e. are differential polynomials in f. Hint: Use 
the formula D-f= q)Dq)~^,f = Dlogq). 
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3. Verify that 

L = D^-u =^ (l'^^ ^D^ -^{uD + Du). 

4. Prove that the change of independent variables yields 

dx = gdy ^ Dy = gDx and U = g^U + {Dy,g{y)} 
where a " Schwarzian derivative" {D,a} is defined as follows 

{D^,a{x)) = -. — -^ — • (4.53) 

One can verify readily that Liouville transformation is invertible and 

lJ = g\lJ-{D,,h{x))) ^ U = h\U-{Dy,g{y)}), g{y)h{x) = \. 

5. For evolution equation 

n-t = Un + F{u,ui, . . . ,ui^), k<n, {n>2) 

possessing a formal symmetry prove that = du^F is a density of a local conser- 
vation law. 

6. Check out that 

L = D -u, Lf = [(L2)+,L] 4^ 4ut + Uxxx^6uUx. 

7. Find non-trivial densities p = p(m,mi) for Burgers equation from Example 
2.2. The picture would be incomplete without next 

8. Prove that condition (g) = g^ (/) coincides with cross derivatives equation 
Ai (") = Ai^fi (w) if DtiU = f and Dt2U = g. 

9. Prove that the bracket {f , f} =^ f} (f) - fl {f ) (Cf Definition 2.8) satis- 
fies the Jacobi identity. 

10. Find the constant e such that Dt^u = etUx+ 1 is the symmetry of equations 
D^u = f from Examples 2.2, 2.3 in the sense that Dt^Dt^u = Dt^Dt^u. 
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Chapter 5 

Integrability of non-linear PDEs 



5.1 Introduction 

In this Chapter we are going to introduce some definition of integrability for non- 
linear differential equations in partial derivatives. (Cf. Intro, in previous Ch.3) It 
links together the notion of higher symmetries of PDEs and the theory of commu- 
tativity from previous chapter. While in Chapter 3 we have deal with a differen- 
tiation Dx with respect to single independent variable x, now we try to insert ad- 
ditional differentiations A^., compatible with and to introduce a corresponding 
set of independent variables ^yt, ^ = 1, 2, 3, . . . . In the first section we address 
this problem algebraically and consider action these additional differentiations 
Df^, upon elements of the ring of formal series in powers of D studied in Sec- 
tion 3.2. It immediately produces a nonlinear system of PDEs for coefficients of 
the formal series under consideration. In particular these PDEs include in itself 
famous Kadomtzev-Petviashvili equation (KP): 



In the next, oriented on PDEs theory, section we consider a broad class of evolu- 
tionary equations 



with one space variable x. We will show how constructively build up a Lax pair for 
equations (5.2) possessing infinite number of higher symmetries. It is remarkable 
that in this case basic theorem is formulated again in terms of formal series in 
powers of D and highlights notions of algebraic theory of the first section from 
the new point of view. 




6UUx)x + 'iUyy. 



(5.1) 



— f{^^ Mxc) Uxxx •••) 



(5.2) 
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5.2 KP-model 

In order to introduce an action additional differentiations acting upon elements 
of the set of PSDO (4.3) it is sufficient to define it action on the coefficients of 
formal series A — Y.<^kD^- We postulate the formula as follows 

A, (A) = an,t,iy + • • • + flo,^, + a-i,tP' ^ + . . . =^ [BkA (5-3) 

where the formal series 5;^, k= 1, 2, 3, ... will be defined exactly below and 
n e Z denotes an order of A. 

Firstly, we remark that the basic differentiation D in the set of formal series could 

def 

be written down in this form and corresponds to the choice B\ = D. Namely, use 
the definition of multiplication of the formal series A = Y^^^kP^ and D we have 

[D,A] =DA-AD = Y,Do a^D^ - a^^D^+i = £ ak,^D'' (5.4) 

Thus, we verified Eq.(5.3) for k = I with B[ = D and A, = D^. 

Secondly, we remind that the Leibnitz rule (general feature of any differentiation) 
is satisfied for the operation of commutation with a formal series B since 

[B,AiA2] =BAiA2-AiBA2+AiBA2-AiA2B=[B,Ai]A2+Ai[B,A2]. 

Thus, left and hand sides of Eq.(5.3) obey Leibnitz rule and we wish now to 
compare orders of the two formal series in the formula under discussion. Due 
Proposition 3.7 the order of commutator [5, A] is equal n + m— I if the formal 
series B of order m is chosen independently from the series A of order n. On the 
other hand the order of Dt^ (A) has to be less or equal of the order n of A . In the case 
m = 1 all agreed since n + m~l=n and that opportunity used in Eq.(5.4). For m > 
1 in order to equalize the orders we must choose at least leading coefficients of Bj^ 
in the strong accordance with leading coefficients of series A. It is the source of 
nonlinearity imprinted in Eq.(5.3). Take in account Proposition 3.5 we formulate 
below very nice algebraic solution of this problem of equalizing orders in left and 
right sides of Eq.(5.3). Some times next proposition is considered as a model of 
modem integrability theory. 

Lemma 4.1. (Sato ??) Let formal series B of the first order is normalized as 
follows 

B = D + biD-^+b2D-^ + b3D-^ + ... (5.5) 

and the differentiations (5.3) of the coefficients of this series defined by formu- 
lae 

Dk =^Dt^{B) = bij^D-'+b2fyD-^ + ■■■ = [bX,B] (5.6) 
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Then the differentiation D\=Dx and DnDm{B) = DmDn{B) for any n, m> 1 , and 
moreover 

DniB'l) -DUBl) = [BI,B1], Dn{B"l)-D^{B"l) = [BlX+i (5.7) 

< Obviously B^=D and, therefore, Di (B) = 5;^. In order to prove (5.7) we start 
with the identity 

= = [Bl,B"^] + [5« ,5"!] - [B"l,B'L] + [B'L,B'!^]. 

On the other hand 

=Z)«(5«+i -B"_) = ] + [B1,BI] -DUBI), 

D„{B"l) = Dn{B'"+^ -B1) = [Bl,B'!^] + [B^^B"^] -Dn{B"l). 

and (5.7) follows from it. It is easy to see now the commutativity conditions 
DnDm{B)=DmDn{B) are corollary of Eq. (5.7). I 

Remark 4.2 In this chapter we do not consider additional differentiations Dt^^ for 
formal series of general form. Nevertheless, the equations Eq.(5.6) and Leibnitz 
rule imply that for powers fi" of the normalized series (5.5) we have 

A,(5'') = [4,5«]. (5.8) 

This generalization Eq.(5.8) and Proposition 3.4 allows to include in domain of 
differentiations formal series (4.3) of any order n but yet normalized (i.e. lead- 
ing and next one coefficients equal 1 and 0, respectively). In order to get off these 
restrictions one should to develop algebraic variant of transformations (4.10 and 
(4. 1 1 used in previous Section yet it is not easy at all. We are going to discuss this 
normalization problem in next section from quite different point of view. 

The nonlinear action on elements of of additional differentiations Z)f^, com- 
patible with Dx defined by Lemma 4. 1 is not unique. In particular, one can start 
with the formal series 

A = D + aQ + aiD^^ + a2D^^ + a^D^^ + . . . (5.9) 

and use another way (Cf. (4.4), (4.5)) to separate polynomial and Lorant parts of 
A": 

A" = A« + , where A« =^ D" + fl„,iD"- ^ + fl„,2£>""^ + • • • + fl„,„-iD 

and n=l,2,3,.... The proof of next proposition repeats the proof of Lemma 4. 1 
and will be omitted. 
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Lemma 4.3 Let fomal series A of the first order defined by (5.9) and the differ- 
entiations (5.3) of the coefficients of this series defined by formulae 

Dk=^Dt^{A) = aofy + aifyD-'+a2fyD-^ + ---^ K,A]. (5.10) 
Then Di = Dx and DnD„t{A) = DmDn{A) for any n, m>l. 

One could notice that as well as Lemma 4.1 corresponds to representation 

as the direct sum of the polynomial and not polynomial parts (see Definition 3.1) 

Lemma 4.3 corresponds to another way do define this direct sum (Cf. (4.6)). 

We are going now to consider the action of differentiations Z);^, on elements 
in more details and write down an explicitly some nonlinear equations for coeffi- 
cients of the formal series (5.9) and (5.5). In order to stress that we take next: 

Definition 4.3 We call by mKP-hierarchy and KP-hierarchy the infinite sys- 
tem of PDEs for the coefficients ai^ and bi^ of the formal series(5.9) and (5.5) 
defined by Eq.(5.10) and Eq.(5.6), respectively. 

Below we provide basic, in certain sense, examples of equations of mKP and 
KP hierarchies. It will be equations for qq and b\ , respectively, with derivatives 

on the first three independent variables ?i,?2,?3- In the case of mKP-hierarchy it 
called modified Kadomtzev-Petviashvili equation (mKP shortly) and is similar to 
KP-equation (5.1) mentioned in the introduction. 

Example 4.5 Historically, for KP-equation (5.1): 

the auxiliary linear problem arised in the form as follows 

\l/y = Vxx-u\l/, m = Vxxx + a\l/x + b\l/ (5.11) 

In order to find the compatibility conditions which gives KP-equation we have to 
equate cross derivatives {y/y)t — {Vt)y- Denote y/k — D'^xi/'wc find after x— differentiation 
Eqs.(5.11) 

e{Yy)t = Y5 + {aYi+bY)xx-£iitV-u{il/3+aYi+b\i/), 
e{\i/t)y = V^5 - iu\l/)xxx + ay\l/i +a{\i/3 - uxj/i - Ux\l/) + byXj/ + b{\i/2 - uxj/). 
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Equate now coefficients by Xj/^ in above formulae one gets for = 2, 1 , 0, respec- 
tively 

2ax + 3ux = 0, 2bx = ay- 3uxx -axx = ay + axe, 
by + euf = bxx + Uxxx + (^Ux- 
Rewrite the last equation in terms of u one obtains KP. 

In terms of Definition 4.3 it corresponds to differentiations (5.3) with A; < 3. We 
have for Dj^ , Df^ : 

mKP: Al = D^ + 2aoD, AI = + 3aoD^ + 3(ai+ao^ + al)D. 
KP: Bl = D^ + 2bi, Bl = D^ + 3biD + 3ibi,x + b2) 

Therefore it coincides in KP-case with (5.11) up to notations. Moreover, Lemma 
4.1 equations 

D2{b) = [bI,b], D3{b) = [bI,b] 

yield 

bi,y = bi^xx + '2b2,y, b2,y ^ b2,xx + 2b3^x + 2bi^xbi- 

bi,t = bi^xxx + 3b2,xx + 3b3^x + bibi^x 
and, thus, exclude b2 and b^ we find: 

1 3 

[bi,t - ^bi^xxx - 3bibi^x]x = 4^i,xy 

Compare these two view points on KP one could say that these are equivalent and 
give the same nonlinear PDE. 

Analogously, using Lemma 4.2 equations one gets mKP 

4qtx = Dx{qxxx-'2ql) + 3qyy + 6qxxqy, ao = qx. (5.12) 



In algebraic approach to integrability the notion of commutativity appears prefer- 
able in comparison with the notion of symmetry. Thus it seems awkward to re- 
formulate Lemmas 4.1-2 in terms of symmetries of the infinite system of PDEs 
constituting hierarchies. Yet consider particular equations of hierarchies like (5.1) 
and (5.12) the notion of symmetry becomes more important and useful. In the 
next example we demonstrate reductions Eq.(5. 1) related with stationary points of 
very simple symmetries transformations. 
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Example 4.6 Since the equations contain independent variables only inexplic- 
itly they are keep its fonn invariant under shift In the case bi not depends in y KP 
equation reduce and gives rise KdV: 

bl,y = ^ 

Bussinesque?? 

5.2.1 Conservation laws 

Interplay of symmetries and conservation laws for ODEs has been discussed in 
Introduction. This interconnection for PDEs is quite different but very interest- 
ing as well and we present here a terse introduction to pure algebraic aspects of 
the theory oriented on applications in next chapters. Hereafter in this section we 
restrict ourself by the case of KP-hierarchy and will denote by ^ the set of poly- 
nomial functions on coefficients bj of the formal series (5.5) and its derivatives 
with respect to i— 1,2, 

Definition 4.7 Corollary of equations of KP-hierarchy written in the form of 
divergence i.e. 

£A,(i8,) = (5.13) 
k 

is called conservation law of KP-hierarchy if sum in (5.13) is finite and the the 
variables fi^e 

In order to build up conservation laws of KP-hierarchy one can use lemma below 
which based on next definition. 

Definition 4.8. For the formal series A e in powers of D 

n 

A= akO'' = UnD" -\ \-ao + a-iD~^ +a-2D~^-\- ... . 

k=—°° 

the coefficient fl_i is called residue and is denoted as res (A). 



Lenuna 4.9 For all m,nez 

Tcs{[air,bD''])=D,am,n 
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where am,n is a polynomial function on a, b and it's x— derivatives. 

Proof. The residue vanish if the powers m, n obey the condition mn >0. For 
instance in the case n = the commutator aD'^b — baD"^ is a differential operator 
if m > and PSDO of order m— 1<— 2ifm<0 (see Prop. 3.7). Obviously the 
coefficient by is zero in both cases. Obviously as well that the residue vanish 
if m + n < 0. 

Let now m, n have different signs and m + n — k>0. Then in virtue of (4. 1) 

res([aD'",Z7Z)"]) = ^ {aD^+\b) + {-lfD''+\a)b) 

since 

m + n = k =^ m{m— I) ■ ■ ■ (m — k) = ±n{n — I) ■ ■ ■ (n — k). 

Standard "integration by parts" proves now the lemma. Particularly for ^ = 0, 1 
we have, respectively 

aD{b) +D{a)b = D{ab) aD^{b) -D^{a)b = D{aD{b) - D{a)b). 

I 

Use (5.8) we find 

A,(res5") = res = res[5^,5"] 

and Lemma 4.9 yields 

Corollary 4.10. KP-hierarchy possess infinite sequence of conservation laws of 
the special form as follows 

def 

D,^{Pn)^D,{ak,n), p„ = resl^"), C7fe,„e^ (5.14) 

where ^ denotes the set of polynomial functions on coefficients bj of the formal 
series (5.5). 

Example 4.11. 

Pl=^l, p2 = 2Z72-^l,x, P2 = 2Z72 - P3 = 3Z73+3Z7f-Z72,x-^l,xc 

A,(pl)=£>x(pn+l) 
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Theorem 4.12 ([?]). Let in the case (5.5) 

oo oo 

A^D + Y^biA', A'"^{A'-)++Y^bkmA-K 

i=l k=\ 

Then for all /, n = 1,2,... 

Dnbi = Dbin, {D = Di). (5.15) 

Proof. Lemma 4. 1 yields 

oo oo 

Dn{Ai)=D{A„) + [Ai,A„], Ai = Y,biA-', A„ = bkmA''' ■ 

i=l k=l 

We have now 

oo oo oo oo 

Dn{Ai) = YDn{bi)A-' + Ybi[A-\An] D{An)=YD{bkn)A-''+Ybkn[A-\Ai] 

i=l i=l k=l k=l 

and 

oo oo 

[Ai,A„] + Y bkn[A-\Ai] -Ybi[A-\An] = 0. 

k=l i=l 



In order to highlight the notion of conservation laws let us consider L— periodic in 
X solutions of KP-hierarchy. Then 

rxQ+L rxo+L 

Dt^{p) = Dx{Gk) Dt^ / pdx= Ok^dx = 0. 

Jxn Jxn 



Therefore, the integral R{i) ^= j^^^^ pdx conserves it's value as times change 



and the integral R = j^^^^^ pdx rather than density p has "physical sense." That 



7X0 



is the reason to call two conservation laws in the differential form Z);^(py) = 
Dx{(yj,k): j = 1,2 with distinct densities equivalent p\ = P2 if 

rXQ+L rxo+L 

/ pidx= / P2dx. (5.16) 

Jxn Jxn 



IXQ JX{i 

In general not necessary periodic case, we take 



Definition 4.11 Two conservation laws of KP-hierarchy Z)f^(p^) —Dx{Oj^k}, j = 
1 , 2 to be said equivalent iff pi — P2 e where latter is subspace of ^ consti- 
tuted by ImD. 

Obviously, equivalence relation introduced in Def.4. 1 1 implies the equality (5.16). 
Moreover, as we will see in the next section this relation p = could be verified 
in fully algorithmic way. 
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5.3 Integrability problem 

In this section we discuss a general definition of integrability for evolutionary 
PDEs (5.2) with one space variable x using the theory of formal series in powers 
of D = Dx (sections 2 and 6). Follow traditions of differential algebra [63] we 
introduce the notation 

M !->■ Ml !->■ M2 !->■. ..l->- MA: • 

for the infinite sequence of derivatives with respect to x of the basic variable u. ^ 
In these notations the general evolutionary PDE is rewritten as follows 



Mf = /(x,M,Ml,M2,...,M„). (5.17) 

The highest derivatives m„ defines the order n of Eq. (5.17) and we consider 
usually equations of order n>2. The dependence on x of the right hand side in 
Eq.(5.17) became very important in the next Chapter in relation with godograph 
type transformations. 



5.3.1 Basic definitions. 

Consider the set ^ of smooth functions of variables x,u,ui,U2, ... we will use 
the equation Eq. (5.17) in order do define action of the "additional" differentia- 
tion Df on (the first one is D = Dx). Namely we take as definitions of these 
differentiations the formulae 

Dx: g{x,u,ui,U2,...,u,n)^ gx + g*iui), ui=D{u), (5.18) 
Df-. g{x,U,Ui,U2,...,Um)i-)- g^{f), f = Dt{u) (5.19) 

where gx = dx{g) is partial derivative of the function g and denotes differential 
operator as follows 

g/Ml^^D^=^ + ^D + ... + 4^D-. 

^ OUk OU OU\ OUm 

It is easy to see that operators (5.18),(5.19) act on the "generators" of the set 
in a natural way, as partial differentiations should do: 

Dx{x) = \, Dx{u) = ui, Dt{x) = 0, Dt{u) = f{x,U,Ul,U2,...,Un) 



'in fore cited book it called differential indeterminate. 
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and 

D^D,{u)^D^{f), D,D^{u)^D,{ui)^D^{f) ^ D^Dt{u) ^ D,D^{u). 

It is important to notice that action of m— th order differential operator on func- 
tions V can be defined also as 

dg[u + ev] 

8*{v)^ 3 le=o- (5-20) 

ae 

In particular as corollary of Eq. (5.20) we find that 

{f8)* = f8*+8f*- 

One can prove now that Dx^ Df and A = DxDf —DfDx as well, obey Leibnitz rule 
and that A = on the set . 



Definition 4.13. The operators (5.18) and (5.19) on the set are called total 
differentiation with respect to x and evolutionary differentiation related with 
f e respectively. The Gato derivative (5.20) of a function g e ^ we call 
linearization of the function g. 

Summing up we have 



Proposition 4.14. For any smooth function / e ^ corresponding evolutionary 
differentiation (5.19) commutate with (5.18). 

This statement looks trivial but if, as in previous Section, we try to find other 
evolutionary differentiations 

Ai(w) = A(-«,W,Wl,W2,---,Mni) (5.21) 

which commutate each other we shall come to nontrivial obstacles. Situation here 
in certain sense corresponds to correlation of the Proposition 3.5 which is very 
general and abstract and more interesting problem about commutating differen- 
tial operators considered in Section 4. The role of Proposition 3.5 plays now 
Lemma 4.1 and commutativity conditions [D^^D^J = have to play the role of 
commuting differential operators. This link became more apparent and instruc- 
tive if we bring up commutativity conditions to the level of a tangent space of 
^ using (5.20). Namely instead of writing down the cross derivatives equality 
Dt^Dt^{u) = Dt^Dty (m) for the two evolutionary equations 



Ai(w) =/l(j^,W,Wl,W2,---,"ni), Dt2{u)=f2{x,u,ui,U2,...,u„2) (5.22) 
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we will consider analogous condition Dt^Dt2{v) = Dt2Dt^{v) for linearized ones 
^hi^) = ifk)*iv), k — 1,2. Translate this condition into condition of commuta- 
tivity of corresponding operators 

[Ai-(/i)*,A2-(/2)*]=0 (5.23) 

we come to next ^ 

Definition 4.15. Let us say that the two evolutionary equations (??) satisfy the 
consistency condition and, respectively, compatibility condition if the cross 
derivatives equation Dt^Dt^iu) — Dt^Dt^ (m) is fulfilled or, and in the second case, 
if differential operators equality (5.23) is valid i.e. 

Ai((/2)*) -Dt.am = [(/i)*, (/2)*] (5.24) 

We remind (Cf. (5.3)) that the action of any evolutionary differentiations on 
polynomials A = '^ajD^, aj e is defined by Dt^{A) = '^Dt^{aj)D-' . 

In order to compare two almost equivalent Definitions 4. 15 we give simple exam- 
ple. 

Example 4.16. Let us consider the pair evolutionary equations as follows 

Dfiiu) = f{x,U,Ui,U2,...,Un), Dt^iu) = Mi- 

Then, apply (5.18) and (5.19) and Proposition 4.14 we get 

Dt,Dt2iu)=Dt,iui)=D,if) = f, + Mui), Dt2Dt,iu)=Dt2if)=Mui)- 

Thus the consistency condition is satisfied iff fx — and / = /(m, mi, M2, . . . , m«). 
On the other hand the compatibility condition (5.24) yields 

/«,x = /«i,x = ••• = /«„,x = <^ / = /(m,mi,M2,...,m„)+X(x) 

where X (x) is arbitrary function of the independent variable x. 



^It can be proved that cross derivatives equation Df^Df^ (m) = Df^Df^ (m) implies (5.23). 
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5.3.2 Main theorem 

We wish now to address more practical question. Let the evolutionary equation 
(5.17) of the order n is given i.e. A(w) = /• What conditions on its right hand 
side / e ^ or, more exactly, on the corresponding differential operator : 

/, = ^ + ILd +... + (5.25) 

OU OU\ OUn 

are needed in order to build up a hierarchy of evolutionary differentiations (m) = 
/yt compatible with fore given. Take in account that in the compatibility condition 

A((A)*) - A,((/)*) = [/*, (A)*] (5.26) 

the order % of the differential operator {fk)* could be arbitrary large whilst the 
order of Dti^{{f)^:) is fixed by n we see that leading terms in Dt{{fk)^) and in 
the commutator in the right hand side should coincides. It gives rise the definition 
which reproduces the main formula (5.3) of previous section in a new environment 



Definition 4.16 We say that evolutionary equation 

Ut= f{x,U,Ui,U2,...,Un) 

of order n posses canonical Lax pair if there is the formal series of order m^O 

A = a^D"' + am- iD""' ^ H h ao + a_ iD~ ^ + a_2^~^ + • • • 

with coefficients aj e ^ such that 

A(A) = [/*,A]. (5.27) 

Some times solutions A of non-zero order of this equation are called formal sym- 
metries. 

Apply this definition to Burgers and Korteweg de Vries equations which we often 
use as illustrative examples, we find 

Example 4.17 For Burgers equation 

Ut = U2 + 2uui = f, f*=D +2uD + 2ui 
the formal symmetry is PSDO as follows A = D + u + uiD~^. Indeed 
[/^,A] =/ + D(/)D~^ and A = D + u + uiD'^ ^ Dt{A) = Ut + ui^tD~^ 
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and thus A (A) = [/*,A]. 

Similarly, in the case of Korteweg de Vries equation (KdV) 

Ut = 6uu\ — U2=f, /* = —D + 6uD + 6ui 



the second order PSDO A = D — Au — 2u\D is the formal symmetry since 

[/^,A] = -4f-2D{f)D-\ and A = -Au-luiD'^ =^ Dt{A) = -Auf-lui^tD'^ 
and thus A(A) = [UA]- 

Next statement can be considered as milestone of integrability theory of evolu- 
tionary equations with one space variable. 

Theorem 4.18 If there are series of evolutionary differentiations A^(m) = fk of 
orders ^ °°hy k °° satisfying compatibility conditions (5.26) for a given 
evolutionary equation Dt{u) — f of order n>2. Then this equation Dt{u) = f 
possess the formal symmetry A of the first order: 



Construction of the solution A of the basic equation Dt{A) = [/*,A] needs some 
technicalities related with the limiting procedure in Eq.(5.26) when ^ — )■ oo. We 
skip that formal part of proof of Theorem 4.18 (see Appendix??) and go instead 
directly to applications. 

Firstly, we explain how to get formal symmetry of the first order starting with any 
solution of Eq.(5.27) and describe it's general solution. Obviously, 



satisfy Eq.(5.27) if Ai and A2 do. Apply results of Section 3 we get. 

Lemma 4.19 Let formal series in D of order m i-e- B — J^^jDJ , bj e is 
solution of Eq.(5.27) and let m 7^ 0, n > 2. Then there exist the first order formal 
symmetry A (see (5.28)) with the leading coefficient 



A — aiD + ao + U-iD +... 



(5.28) 



aiAi + a2A2, ,ajec and A1A2 



«1 = ifn) " , ft 



def df 



(5.29) 
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and any solution B of Eq.(5.27) of order m can be represented as formal series in 
powers of A with constant coefficients a/ e C: 

k<m 

Proof. Slight generalization of Lemma 3.2 yields the formal series with coeffi- 
cients from such that B'^ = B. Let us consider the formal series R = Dt{B\) — 
Since 

Dt{B) = Dt{Bi)B'^-' +BiDt{Bi)B'^-^ + ■■■ +B'^-'Dt{Bi), 

and 

[f,,B] = [f,,A]A'--'+A[f,,A]A'--^ + ---+A'--'[f,,A] 

we have 

Bt - [f^,B] = RA""-^ +ARA'"-^ + ■■■ -{-A'^-^R = 0. 

In this sum the all m series in have the same leading term and, hence, this 
leading term is zero. It is possible only if R = and therefore we proved that Bi 
is a solution. 

In order to find the formula for leading coefficients of solutions Eq.(5.27) one 
should find the leading term in —Bf. Since n > 2 it is contained in 

= C„+«_iZ)"+'""^ + . . . , Cn+m-l = nfnD{bm) - mbmD{fn) 

where /„ and bm are leading coefficients in and B (Cf Proposition 3.7 ). Van- 
ishing leading coefficient gives the equation 

nDlogbm = mDlogfn <^ nlogbm — mlog/n = const 

which implies for m = 1 the formula (5.29) for the leading coefficient of the first 
order formal symmetry. I 

Most known applications of Theorem 4.18 related with next proposition quite 
similar to Corollary 4. 10: 

Corollary 4.20. In conditions of Theorem 4. 1 8 the evolutionary equation (5.17) 
possess infinite series of conservation laws as follows 

A(resA'')=resA(A") = res[/*,A"] e ImD, n = -1, 1,2, . . . . 

Hereafter ImD denotes the subset of functions ^ G which can be rewritten as 
g — Dx{h),h e In other words Corollary 4.20 provides series of conservation 
laws : 

Hpf 

A(pn)=£>x(o-„), p„,a„e^, p„ = res(A«). (5.30) 
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of the equation (5.17)) under consideration. This series should be completed by 
the conservation law with the density 

Po= • (5.31) 

Indeed, rewrite equation (5.27) for m— th order formal symmetry 

A = UmD"" + am- iD'"" ^ H h ao + a_ iD ^ + a-2D^^ + ■■■ 

as follows 

A(A)A-i = /, -A/,A-i = [A-\Af,] 
we find res(AfA^^) G ImD due Lemma 4.9 and that 

res(A,A-i) = res [{am,tD"' + am-i,tD"'-^){a-^D-"' + aD-""-^)] = 
am,t{0C-m^2r)^ • 

It remains to notice that 

^mjc ^m— 1 

a = . 

Example 4.21. For n = — 1 we have due (5.29) 

resA~^ =res {^D-^ + ..^ = {/„)'". 

Thus p 1 = (^^j " and by Corollary 4.20 the evolutionary equation (5.17) 
should always satisfy the condition 



A(p-l)=A(^^J elmD (5.32) 
if it is compatible with hierarchy of evolutionary differentiations. 



Definition 4.22. The series of conservation laws (5.30)), (5.31)) we call inte- 
grability conditions. 
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Definition 4.23. The conservation law is differential corollary of (5.17)) written 
in a form of divergency 

Dt{p)=D,{a), p,ae^ (5.33) 

where functions p and o are called density and flux, respectively. It is called 
trivial if p e ImZ) and there is h e 'W such that p = Dx{h). In this case 

A(P) = Dt{D,{h)) = D,{Dt{h)) e ImD 

for any evolutionary differentiation Z)j. Two conservation laws Eq.(5.33)with the 
densities pi and p2 are considered equivalent if pi — p2 G ImD. 

Come back to Example 4.21 we see now there is an alternative or (P-i) G ImZ) or 
Dt{p-i) G ImD. In both cases this condition is non trivial one and can be verified 
using classical statement by (Helmgolz??) proven in P.Olver book [59]. 



Theorem 4.24. The function h&'^ belongs to the subspace ImD iff E(h) =0 
where Euler operator E is defined as follows 

fief d d o d id 

£^1^ -D— +d2— -d3— + .... (5.34) 

OU OU\ OU2 OU2 



Equations similar to Eq.(5.27) play major role in modem integrability theory and 
are called commonly by Lax pair. The first test which have to justify that that is 
not trivial one ^ could be formulated similarly to the definition 5— integrability be- 
low. In other words proper Lax pair must produce series of non trivial conservation 
laws. 



Definition 4.25. The evolutionary equation (5.17)) with formal symmetry (5.27)) 
is called 5— integrable if the series (5.30)) contains nontrivial conservation laws 
of arbitrary higher order and C— integrable in the opposite case. 

The famous Korteweg de Vries (see (2.163) in Section?) gives us an example 
of 5— integrable equation. As matter of fact the formal symmetry for KdV indi- 
cated in Example 4.17 related very closely with Lemma 4.6 ( see Eq.(2.89) and 
concluding remark in Section?). KdV hierarchy will be discussed in next section. 



'There is known a few cases of faked Lax pairs. 
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Although Definition 4.25 looks too complicated to be useful it is not the case. In 
the next chapter we apply this definition to the second order evolutionary equa- 
tions (5.17)) and verify that first three of conditions (5.30) fully describe all sec- 
ond order evolutionary equations possessing hierarchy of higher order symme- 
tries. Burgers equation from Example 4.17 is, in certain sense, typical representa- 
tive in this class and a class of C— integrable equations in general. 



Symmetries and cosymmetries 

Conclude this chapter we discuss duality between symmetries and conservation 
laws starting with "zero" dimensional evolutionary equations i.e. dynamical sys- 
tems 

^ = /(m), u={u\u\...yf (5.35) 

where ()^ denotes hereafter matrix conjugation. In this case the conservation law 
with the density p = p (m) is just a constant of motion 

Infinitesimal version of conservation law is obtained by replacement uhy u + ev 
in above equation (Cf Definition 4.13): 



Proposition 4.25. Let p = p{u) is the density of conservation law for Eq.(5.35)) 
and p denotes column- vector constituted by the partial derivatives ^ . Then 

A+/Jp = 0, {f*)ij=^ (5.36) 
Proof. The dynamical system Eq.(5. 35)) yields 

^-/*v (5.37) 
dt 

for infinitesimal dynamical variables v. The same change m by m-|- ev gives in the 
infinitesimal approximation 

p{u + ev) = p{u) + ep^v. 

Substitute this expression into conservation law with the density p we get 

d ^ _ dp _ ^ dv 

<p,^;>=<^,i;> + <p > (5.38) 

dt dt dt 
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where bracket denote the scalar product. Take into account Eq.(5.37)) one gets 
Eq.(5.36)). I 

In the case of evolutionary equation (5.17) the linearization (5.20) gives for in- 
finitesimal variable v PDE as follows 



vt = f*v ={fu + fu,D + ... )v. (5.39) 
Similarly to zero-dimensional case Definition 4.22 of conservation laws implies 

Proposition 4.26. Let p = p (x, m, mi, M2, . . . ) is the density of conservation law 
(5.33)) for evolutionary equation (5.17) and 

p = p„ -Z)(p„J +D\pu,) (5.40) 

Then 

Dtp+f^p = (5.41) 



where /J is the differential operator 
formally conjugated to the operator 

The proof of Proposition 4.27 (see [59]) is analogous to the zero-dimensional case 
above. The function p defined by Eq.(5.40) is called variational derivative of the 
function p. One has to notice that in the L— periodic case for the variation of the 
functional 

def Fo+i^ 
Rlu] = / p(x,u,u\,U2,- ■ ■)dx 



def f^o+J^ 
l[u] = / p{x,u,ui,U2,- ■ ■)d^ 

Jxn 



JXO 

we have in the infinitesimal approximation 

/fXo+L fXo+L 
P:t:VdX = £ 
Xo Jxo 

The last formula in particular plays the role of the scalar product (5.38). 



R[u + ev]-R[u] = e p^vdx = el pvdx. (5.42) 

Jxn Jxn 



Compare Propositions 4.26 and 4.27 one can notice that in both cases p e KerZ)f + /, 
while symmetries are in Ker(Z)f — /^). 



For future references we formulate 
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Definition 4.28. The order of nontrivial conservation law (5.33)) with the den- 
sity p is defined as order of variational derivative (5.40). 

We remind that density p is called trivial if p e ImD and in virtue of Theorem 
4.23 variational derivative p = in this case. One can verify that the order n of 
variational derivative p always is even number and if n = 2m > then 

= a{x,U,Ui, . . . ,Um)u2m + ■ ■ ■ , (-1)'"^-^ ^0. 

For n — one can assume without loss of generality that p = u. The latter case take 
place iff the right hand side / of evolutionary equation (5.17) satisfies a condition 
/ e ImD. 



5.4 Summary 

The primary goal of this chapter is to establish a workable criterion that can be 
readily checked to determine whether a given evolutionary equation possess in- 
tegrability features like symmetries, local conservation laws and Lax pairs. A 
systematic linearization of all equations is one of basic tools in this chapter. 



5.5 Exercises for Chapter 4 

1. For evolution equation 

Ut = u„ + F(u,ui,...,Uk), 

possessing a formal symmetry prove that 
vation law. 

2. Check out that 

L = D^-u, Lt = [{Li)+,L] 4ut + Uxxx = 6uUx. 

3. Find non-trivial densities p = p(m,mi) for Burgers equation from Example 
2.2. The picture would be incomplete without next 



k <n, {n>2) 
— du^F is a density of a local conser- 
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6. Prove that condition = g*{f) coincides with cross derivatives equation 
Dt^Dt^{u) = Dt^Dt^ (m) if Dt^u = f and Dt^u = g. 

7. Prove that the bracket {f , f} =^ /J (f) - fl {f ) (Cf Definition 2.8) satis- 
fies the Jacobi identity. 

8. Find the constant e such that = etUx + 1 is the symmetry of equations 
D^u = f from Examples 2.2, 2.3 in the sense that D^Dt^u = Dt^D^u. 



Chapter 6 

Burgers-type equations 



6.1 Introduction 

We are going to prove follow by [64] that any evolutionary equation of the second 
order 

Ut = F{x,u,ui,U2), Uj = D^{u) (6.1) 

which possess a canonical Lax pair can be reduced by invertible transformation to 
one of the following three canonical forms 

Ut = D{ui+u^+f{x)) (6.2) 
Ut = Z)(^-2x) (6.3) 

Ut = + eixu + e2u) (6.4) 

or to the linear equation ut = U2 + f{x)u with the arbitrary function f{x). How to 
do it constructively??? 

(6.2)-(6.4) 

6.1.1 Rough classification 

The special feature of evolutionary equations of the second order which simplifies 
a general structure of integrability conditions is defined by 

Lemma 1 A nontrivial conservation law (5.33)) for evolutionary equations of 
the second order (6.1) has order n = 2m less or equal 2. 
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Proof. One can use the formulae from previous Chapter 

P ^ Pu-D{pui)+D^{pu2) =P ^ a{x,U,Ui, . . . ,Um)u2m + ■ ■ ■ 

and the equation (5.41) for variational derivatives 

Dfp+Fjp = 0. 

Separate in this equation leading terms which contain m„+2 one finds 

o def dF 

air {F)+D\aF2Un) + ---^0, F2 = ^ (6.5) 

and if the order n = 2m>2 that results in 2aF2Un+2 = a = 0. 1 

Thus the order of conservation law for evolutionary equation (6.1) can be equal 
or 2. If order n = 2m = 2 Eq.(6.5) in the proof of Lemma gives 

2F2 + U2F22 = ^ F = {fu2 + g)~^+h (6.6) 

where functions /, g and h may depend x, w, wi . 

In the zero order case with p = p{x,u), p = Pu'/'Owe have 

DtP = PuF elmD =^ F = f{x,u,ui)u2 + g{x,u,ui). (6.7) 

— 1/2 

Since p_i = F2 should be density of conservation law of order less or equal 
2 (see Lemma 1) this function on w mi, U2 should be linear in U2 and we have 
three possibilities: 

— 1/2 

F2 ' = Da{x,u,ui), (6.8) 
F^^^^ = Da{x,u,ui) + P{x,u), j8„y^0 (6.9) 

— 1/2 

F2 = Da{x,u,ui) + ^{x,u,ui), j8„j„i7^0. (6.10) 

The first case corresponds to trivial conservation law and it has order and 2 in 
the cases (6.9), (6.10), respectively. 

Rough classification of Eqs.(6.1) uses only the first integrability condition Z)f(p_i) e 
ImD. 

Svinolupov's Theorem 2. Up to invertible change of variables 

d u 

x= q>{x, u, ui), u = \i/{x, u, ui) , ui = \i/i{x, U, Ul) = — (6.11) 
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the general nonlinear evolutionary equation (6.1) of the second order can be re- 
duced to one of two quasi-linear forms 

Ut = U2 + f{x,u,ui)j (6.12) 
ut = D{^+f{x,u)). (6.13) 

if it satisfies the first integrability condition. 



Outline of the proof. Since the first mtegrability condition (i.e. local conservation 

-1/2 

law with the density F2 ) is invariant under invertible transformations (6. 1 1) it 
is sufficient to find change of variables such that 



F2^^^dx = dx, or F2^^^dx = udx. (6.14) 

The first case corresponds to trivial conservation law with the density (6.8) and 
the second one to cases (6.9) and (6.10). 

In order to find the change of variables in the case (6.9) it is sufficient to use point 
transformation: 

x=q>{x, u), u = \i/{x, u) =^ udx= \\f{<Px-\-u\<Pu)dx. 

Indeed, compare formulae (6.9) and (6.7) we see that a{x, u,u\) = a{x, u) in (6.9) 
and we have 

\if{(px^-u\(pu)=Da^-^ = uiau + ax + ^. 
Thus, functions 9 and X]/ should satisfy two equations as follows 

OCx + P = ¥(px, OCu = ¥% — — = — • 

«« % 

That yields the first order PDEs for the functions (p and \j/ : 

def 

L<p = 0, L\i/ = Pu¥, L = {ax + P)du-0Cudx. 

In order to find exact expression in quadratures of the function i// in terms of 9 
one can substitute y = 4>(9,m) into the equation Li// = /3„V^. That results in the 
ODE equation ( + j8 )4>„ = j8„^ in which x— variable plays the role of parameter. 
Jacobian of the transformation x = 9, m = ^> is ^x'^u 7^ 0. 

In the case (6.8) the transformation (6.11) which resolve the first of equations 
(6.14) is called contact one. Contact transformations are fully defined (see next 
section) by the first two equations in (6.1 1) 



X= q>{x, U, Ml), U= \l/{x, U, Ml) 



(6.15) 
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which have to satisfy one and same equation 

+ = (6.16) 

with arbitrary coefficient a = a{x,u,u\) depending on three variables x, m, mi. In 
our case Eq.(6.14) gives,obviously, (p = a and Eq. (6.16) for ^ can be rewritten 
as follows: 

.dv/ da , , 

In the case (6.10) the corresponding contact transformation (6.11) is defined by 
the second of equations (6.14): 

ax + uiau + u2aui+ P = \i/{x,u,ui)[(px + ui% + u2%i]. 

That implies 

V(pui = aui, ¥(Pu = 0Cu + pui 
since due Eq.(6.16) we have 

(PxYui - (Pui ¥x = "1 (<P«i Vu - (puVui)- 



du\ du du\ du\ du 



6.1.2 Canonical densities 
Lemma 2 Canonical densities are: 

p_i = F-''^, (6.17) 

Po = Fl/F2-(7_lF2"'/^ (6.18) 
Pi = F-''\]^F,+ ^-ao + \al,-y^^\lF,-{D{F2)f)), (6.19) 

and 

dF dF dF 

^2 = ^, Fi = — , Fq=^ 

OU2 ou\ du 

and C7_i , c7o have to be recursively found. 
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Remark . The integrability conditions 



Dt{pk)^D,{Ok), -1,0,1. 



(6.20) 



are based upon consistence conditions for 

Ut = F, Ux = G, or more precisely Vj = (v) , = G* (v) 
where = F2D^ + F\D + Fq and G* defined analogously: 



It has been proved by [64] that evolutionary equations of the second order (6.1) 
which possess the symmetry of order m > 3 satisfies integrability conditions (6.20). 

Definition 



6.2 Linear case 

In the general theory one can suppose without loss of generality the leading coef- 
ficient Mo = 1 up to change of independent variables introduced by 

Definition 3.4. Two potentials U{x,X) and U{y,X) are said to be related by 

Liouville transformation if there is function g — g{y) such that 



d d d . , 

-7-^* "I" G*F* = -^F^ -\-F^G^ =^ —L = G*J 



dx = gdy ^ Dy = gD, and ?7 = g^U + {Dy,g{y) } 



where a " Schwarzian derivative" {D, a} is defined as follows 



{Dx,a{x)} 



def 3a|_ 

4a2 2a' 



(6.21) 



One can verify readily that Liouville transformation is invertible and 



U^g\U-{D„h{x)}) ^ U^h\U-{Dy,g{y)]), g{y)h{x) 



1. 
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Let us consider change of variables x, m, mi, . . . defined by the formula (6.22): 

X — ^{x, U, Ui, . . . ,Um), U = \j/{x, U, Ui, . . . ,Un) (6.22) 

and the next condition of invariance of the vector-field D 

Do^=aD, (6.23) 
where 4> notates the transformation defined by (6.22) and 




Recall that D — Dx is the total differentiation with respect to x of functions of 
variables x, u, ui = Ux, 

Applies 6.23) to Jc, U one finds, respectively, 

D{q)) = a, D{\\f) = aD{u) = au\. 

Therefore, 

ux=D{\if)/D{<p). 

Express u\ as a function of x, m, mi,... and apply (6.23) to u\. It gives M2 = 
D{u\)/D{(p). Thus (6.23) is equivalent to formulae 

Ui+i=D{ui)/D{((>), / = 0,1,.... 

which together with (6.22) define the transformation totally. 

<I> : (x,M, ai,M2, . . .) (x, a,Mi,M2, • • •)• (6.24) 

The order of the transformation (6.24) is grater from m" n'm formulae (6.22). 

It is easy to show that the condition (6.23) is equivalent of invariance conditions of 
= : G)i = dui — Ui-^idx contact differential forms. Reformulation on multi-dim 
case is also easy. 

For zero order transformations {point transformations) 

X = q>{x, m), m = \lf{x, m), 

jacobian of (6.24) d{x,u)/ d{x,u) ^ 0. Less obviously, that Invertability condi- 
tions for a general (6.22) demand that order of 4> does not exceed first one and 
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moreover that transformation should be closed in the space of variables x, w, ui. 
Jacobian 

d {x, u,ui)/d {x, u, u\ ) 

of the transformation 

X = <p{x, U, Ml), U = \lf{x, U, Ml), Ml = Xi^J "l)) (6.25) 

where x — /D{(p) " (p^x — Wui can be rewritten as follows 

¥u-X(pui^O- (6.26) 

Defined by (6.25), (6.26) transformations (6.24) are called contact. Legandre 
transformation, provides an example of the contact transformation. 

6.3 Contact Transformations 

The action of the differentiation Dx defined by the infinite dynamical system 

Dxix) = l, Dxn = ui, DxUk = Uk+i. (6.27) 
The corresponding vector field in scalar case is 

d ^ d 

ox Q OUk 

Contact transformation is described as invertible transformation 

X = q){x,u,ui, . . .), u= \i/{x,u,ui, . . .), Ml = VAi(x, M, Ml, ...),.. . (6.28) 
which leaves the dynamical system (6.27) invariant. This conditions means that 
¥k+i^Dx{\irk)/Dx{(p), k = 0,1,2,... 

or 

D;cA(x,M,Mi,...) =Dx{(p)DjcA{x,u,ui,...). (6.29) 
We have the classical result 
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Proposition. Invertible transformations (6.28), (6.29) are 

point transformations x= ^{x,u),u= u) and contact one with 

¥ui (mi<P« + (px) = (pui {uiYu + ¥x), %¥ui - ¥u(pui 0. 

Composition of two contact transformations is contact again. 

In order to define a contact transformation one can fix the function a — a{x^ w, u\) 
and choose the functions 9, l//^ as two distinct solutions of the first order PDE 



0JC + Ml 0M = a{x, U, Ml)0„i . 

The usual example of the legandre transformation with q) = ui,\j/ = xui — u cor- 
responds to a = 0. On the other hand the infinitesimal description of the one pa- 
rameter group of contact transformations is equivalent to more common theory of 
the first order PDE 

Ut =A{x,u,Ux). 

One-parameter group of CTs (contact transf.) are generated by the first order PDE 

Ut ^ F{x,u,Ux), (6.30) 
and corresponding ODEs (bi-characteristic system of eqs.): 
dx du „ du\ ^ 

^ = -7, ■^ = ^-')'"i' -^=Dx{F)-ru2 = F, + uiFu (6.31) 

where 7= dF /du\. 

Let us consider the vector- field 

d d 
Dt = F^+D,{F)^ + ... 
du oUx 

related with (6.30). It is easy to see that the CT 

t = t, x=^{x,u,ui), u = \(f{x, u, ui), 

where ^(x, m, ui) is any first integral of (6.31) reduces F to the form F = 
a{x,u)ux + b{x,u). Related to the latter quasi-linear case the dynamical system 
(Cf. (6.31)) simplifies: 

dx r \ du . , 

— = — fl(x,M), — = b{x,u). 
dt ^ ' dt ^ ' 
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If q){x, u), 9„ 7^ is first integral of these then the point transformation 

t = t, x=^{x,u), m = i/a(jc, m), 

reduces F to F = b\\ru — a\\rx = c(x, u) . In other words, a composition above trans- 
formations (CT and point one) reduces PDE (6.30) to ODE Ut — c{x^ u). 

Contact Transformations, end 

Backlund problem. Invertability conditions for (6.24) in the case of two j8m = 
(m, v) or more fields. 

Hodograph and legandre 

{X = M, ( X = Ml, 

U = X \ U = XUl — u. 

are involutions. 

For differential substitutions (6.22) of the first order x = D{\lf)/D{(p) generally 
depends of U2 but in the scalar case it cancelled for invertable transformations. In 
the case jSw = (m,v) invertability does not guarantee cancellation j8m2 in jSmj = 
D{\i/) /D{(f>). For instance it is the case for the transformation 

_ Vx _ Vx 

x—x^ u= — , v — —v-\-u — 

which is involution since Vx = uux- Notice, that h{fiu) = uvx/ux keeps invariant 
under this transformation. 

Generalization of criterion of inversability of DS of the first order, indicated above, 
on the case j8m = (m,v) is interesting open problem. For degenerated DS there 
is invariant of the first order like example above. For non-degenerated DS we 
assume inversability and the following condition: 

order of composition of DS is sum of individual ones. 

In this case DS defines shift transformation n ^ n+l, and corresponds with the 
chain of equations like follows 

= a(M„)(v„+i - Vn), Vnx = b{Vn){Un-l - Un). 

6.4 Scalar differential substitutions 

Let us say that there exist scalar DS w = P{x, u, Ux) from 

ut = f{x,u,Ux,Uxx,---), (6.32) 
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into 

Wt = g{x,w,Wx,Wxx,---), (6.33) 

ifP„^y^Oand 

{Pu^D, + Pu){f)^g{x,P,D,{P),Dl{P),...). (6.34) 

Differential substitutions (DSs) defined by module of point and contact ones. We 
say that DS ai from (6.32) into (6.33) is equivalent DS a2, if ai = j8 o a2 o 7, 
where j8 and yinversable transformations eqs. (6.32) and (6.33), respectively. 

Let us relate with P{x, u, Ux) the hyperbolic equation 

p 

Uxy + b{x,u,Ux)uy = 0, b=^, (6.35) 

and its linearization 

[DxDy + bDy + Uy{bu + ^W.^);,)]/ = {Dy + Uybu,) O {D^ + b)f = 0, (6.36) 

where / = /(x, w, Ux-,... and Dx-, Dy denote differentiations / with respect x and y 
on solutions of the equation Uxy = —Uyb. Recall that solutions of (6.36) are sym- 
metries of the hyperbolic equation. 

It is easy to see that DyP = and that (6.36) follows from (6.34). Thus if 

w = P{x,u,Ux). (6.37) 

is DS from (6.32) then / is symmetry of corresponding hyperbolic equation. In- 
verse is also true (see [?]). 



Definition. Substitution (6.37) called regular if the corresponding hyperbolic 
equation has as x— invariant as well as DyP{x, m, Ux) = j— invariant DxQ{y., m, Uy,Uyy,...) 
0. 



Theorem [?]. DS (6.37) is regular iff 

Ux = oc{x,w)u^ + P{x,w)u + y{x,w). 

Examples: 

P = Ux, P = Ux + e", P = Ux-u^, P = Ux + e" + e' 
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It can be proven that any DS of the form 

W ^ \}/{x,U,Ux), X^(p{x,U,Ux) 

is equivalent with (6.37). 

Differential substitutions. 

The problem of classifications all evolutionary equations (5.2) of orders 2 5 

with formal symmetries considered by many publications. Progress here depends 
mainly from a development an appropriate theory of the differential substitutions. 

We say that two evolutionary equations 

Mf =/(a:,M,Mi,M2,...), Vt=g{x,V,Vl,V2,...) 

are related by the differential substitution 

u = h{x,v,vi,. . . . 
if in result of total linearization 

(pt=M<P), (p = h*i¥), ¥t = 8*{¥)- 
one arrives to the equations 

Dt-g*=h~^{Dt-f^)h^. ^ {h^)t + Kg*= f^K. (6.38) 

which could be rewritten in (x,m) as well in {x^v) variables. That implies, in 
particular, that the two equaations connected by differential substitution have the 
one and same order. 

In the simplest case 

Ut=Dg{u,Ul,...), Vt = g{vi,V2,...), u = D{v) 
having in mind the general formula (6.38) one proves 

2.9 Proposition Let A be a formal symmetry of the evolutionary equation ut = 
Dg. Then the similarity transformation 

A = DAD-^ ^ A = D-^AD 

defines the formal symmetry for the equation Vt = g.v = u\. Vice versa if A is the 
formal symmetry for latter equation and coefficients of that series can be rewrit- 
ten in terms of u — v\ (i.e. invariant with respect shift in v) then the similarity 
transformation yields formal symmetry for ut = Dg. 
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Proof. Using the definition (5.20) we find that 

f = D{g) ^ f,=Dog, ^ Dt-f,=D{Dt-g,)D-\ 

Differentiating A = D~^AD we find that [Z)j —g*,A] —0 and vice versa. 
I 

This proposition allows one, in particular, to find the formal symmetry indicated 
in Example 2.2. Namely, the Burgers equation has a form Ut — D{u\ + u^) and in 
result of integration one obtains the equation 

which equivalent to linear one 

Vt = Vxx + vl, W = e", Wt = Wxx- 

Starting with the formal symmetry A=D for the last equation we find in notations 
of Proposition 2.9 

A = e~''De'' ^D + vi ^D + u, A^DAD^ ^D + u + uiD~'^. 

Closing this chapter we provide the list the most popular third order evolutionary 
equations of the form 

Ut = Ut, +F{x,U,Ul,U2) 

which posess the formal symmetry. 

Ut = U3 + P{u)ui, P"' = 0, (6.39) 
Ut = M3-^M^ + (ae^" + j8e"2")Mi, (6.40) 

Ut = U3-:^{ul + Q{u)), Q^ = 0. (6.41) 



6.5 Exercises for Chapter 5 
1. 



2. 
3. 
4. 
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